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Abstract 

The non-invariant energy functional of the preceding paper is improved in order 
to obtain its gauge-invariant form by strictly taking into account the non-Abelian 
character of Relativistic Schrodinger Theory (RST). As an application of the results, 
the dichotomy of positronium with respect to singlet and triplet states is discussed 
(ortho- and para-positronium). The degeneracy of the ortho- and para-states oc- 
curs in RST if (i) the magnetic interactions are neglected (as in the conventional 
theory) and (ii) the anisotropy of the electric interaction potential is disregarded. 
In view of such a very crude approximation procedure, the non-relativistic positro- 
nium spectrum in RST agrees amazingly well with the conventional predictions. 
Indeed using here a very simple trial function with only two variational parameters 
yields deviations from the conventional predictions of some 10% (or less) up to large 
principal quantum numbers (n ~ 100). This hints at the possibility that the exact 
RST calculations could yield predictions for the bound systems which come very 
close to (or even coincide with) those of the conventional quantum theory. Such 
a numerical correctness of the RST predictions supports a fluid-dynamic picture 
of quantum matter where any single particle appears as the carrier of a wave-like 
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structure (such a structure, however, is generally believed to be associated exclu- 
sively with statistical ensembles but not with single particles). Here, more detailed 
predictions of the positronium spectrum would be desirable but this will necessitate 
to first elaborate more exact solutions of the RST eigenvalue problem, especially 
with regard to the anisotropic and magnetic effects. 



PACS Numbers: 03.65.Pm - Relativistic Wave Equations; 03.65.Ge - 
Solutions of Wave Equations: Bound States; 03.65.Sq - Semiclassical 
Theories and Applications; 03.75.B - Matter Waves 
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I Introduction and Survey of Results 



The present paper continues the investigation of the fluid- dynamic aspects of quantum 
matter which, however, in the literature is mostly described in terms of probabilistic point- 
particle concepts. The frequently quoted paradigm of "wave-particle duality" seems to 
suggest a certain equivalence of the wave and the particle picture, but in practice the 
particle picture appears preferably in combination with the probabilistic interpretation 
and thus has overrun the fluid-dynamic approach. But here one should not overlook the 
fact that the probabilistic point-particle picture has its interpretative difficulties, too. 

Indeed, even one century after the first tentative steps [1] towards a quantum theory 
of microscopic matter there is still considerable uncertainty about the true nature of the 
quantum phenomena. Despite the existence of a powerful mathematical apparatus for the 
description of the quantum world the controversy persists about the right understanding 
of what is really going on during the course of a quantum process. Even those workers in 
this field, who are able to handle excellently and successfully the mathematical apparatus, 
apparently like to conceive the quantum process as a "mystery" from the conceptional 
viewpoint [21 E]- Thus the proponents of the probabilistic point-particle picture have to 
concede (more or less involuntarily) that the core of the philosophical confusion about 
the quantum behavior of matter refers to just that notorious wave-particle duality which 
Feynmann seemingly had in mind when he said: " ...which has in it the heart of quantum 
mechanics. In practice, it contains the only mystery [1]". 

Concerning now this mysterious wave-particle duality, there seems to have evolved a 
certain wide-spread consensus which may perhaps be best expressed as follows: "It is 
frequently said or implied that the wave-particle duality of matter embodies the notion 
that a particle-the electron, for example-propagates like a wave, but registers at a detector 
like a particle. Here one must again exercise care in expression so that what is already 
intrinsically difficult to understand is not made more so by semantic confusion. The 
manifestations of wave-like behavior are statistical in nature and always emerge from the 
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collective outcome of many electron events. In the present experiment nothing wave-like 
is discernible in the arrival of single electrons at the observation plane. It is only after 
the arrival of perhaps tens of thousands of electrons that a pattern interpretable as wave- 
like interference emerges. " [5]. According to this view, the wave-particle duality refers to 
statistical ensembles but not to individual "particles" . As a consequence of this general 
belief, the wave concepts in quantum theory must not be applied to single particles: " Al- 
though one can in principal measure the mass, charge or energy of a single electron (held, 
for example, in an electromagnetic trap), one cannot measure its de Broglie wavelength 
except by a diffraction or interference experiment employing many such electrons" . [5]. 

But in contrast to this generally accepted picture of the wave-particle duality, it seems 
to us that those one-particle self-interference phenomena [21 E] do distinctly hint at the fact 
that already any single particle is equipped with some kind of wave-like structure which 
becomes manifest in some special situations (in other situations the wave-like structure 
may remain hidden). Naturally, such kind of philosophical basis must motivate one to 
look for a fluid-dynamic description of quantum particles which is able to compete in 
many areas with the conventional probabilistic paradigm. Here the conventional prob- 
abilistic theory of quantum matter and the present fluid-dynamic approach should not 
be understood as mutually falsificating logical systems but should rather be conceived in 
the sense of Bohr [6] as complementary ways to look upon the quantum domain. Indeed, 
the very successful density functional theory [HI [9] demonstrates that a fluid-dynamic ap- 
proach to the quantum phenomena is logically possible and at least as useful in atomic 
and molecular physics as the conventional probabilistic approach. 

An alternative theoretical framework of the fluid-dynamic type has recently been pro- 
posed in form of the Relativistic Schrodinger Theory (RST) (see the list of references at 
the end of the paper). This (rigorously relativistic) theory is also based upon the con- 
cept of wave functions, albeit not in the probabilistic sense of the conventional theory, 
but rather in the fluid-dynamic sense where the wave functions serve to build up the 
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physical densities of an N-particle system (e.g. the charge and energy-momentum density, 
etc.). But clearly, the mathematical formalism of such a fluid-dynamic approach must 
necessarily differ considerably from the probabilistic calculus of the conventional theory 
(i.e. Whitney sums of fibre bundles vs. tensor products of Hilbert spaces); however, de- 
spite this mathematical difference the physical phenomena must be describable equally 
well in the fluid-dynamic picture. For instance, the outcome for the energy of a bound 
N-particle system should be the same, no matter whether it is calculated with reference 
to the probabilistic or to the fluid-dynamic terminology. In the first case, one has to look 
for the eigenvalues of the N-particle Hamiltonian H(l . . . N); and in the second case one 
has to look for the values of the total energy functional Em upon the solutions of the 
RST eigenvalue problem, which itself consists of the coupled mass eigenvalue equations 
and the gauge field equations (see below). Obviously, there arises in both cases a very 
similar problem, namely (i) either to define the Hamiltonian H(l . . . N) for a relativistic 
N-particle system [TQl E] or (ii) to find the right relativistic energy functional Em for 
such a system (which is the main subject of the present paper). The latter object E[t] 
does then not only specify the energy carried by the solutions of the RST eigenvalue 
problem but simultaneous it provides also the basis of the principle of minimal energy, 
i.e. the variational principle 

SE [T] = , (1.1) 

whose extremal equations are nothing else than just those coupled field equations of the 
RST eigenvalue problem. 

The goal of the present paper is now to set up the RST energy functional -Ejt] for 
bound two-particle systems (N = 2); the case of general particle number N might then 
be discussed along quite similar lines of reasoning. Here, an important point refers to 
the question whether the considered N-particle system consists of identical or different 
particles. From the principal point of view, the interaction of the particles occurs always 
via the mechanism of minimal coupling; i.e. the N-particle wave function \l/ enters the 
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RST field equations in form of its gauge-covariant derivative T> 

V^ = d^ + A^ . (1.2) 

The wave function \1/ itself has the status of a section of a spinor bundle; and the nature 
of the bundle connection A^ ( "gauge potential") determines the type of interactions be- 
tween the particles. For instance, when the particles are subjected to the electromagnetic 
interactions one splits up the structure algebra U(N) (as the Lie-Algebra valued range of 
the connection A^) into its maximal Abelian part a(N) = u(l) © u(l) © ... © u(l) and 
the complement b(N) 

u(N) = a{N) © b{N) . (1.3) 

The RST philosophy says now that the electromagnetic interactions between the par- 
ticles are always to be described by the projection of the connection A^ to the Abelian 
subalgebra a(N) (i.e. A^ — > «4 M | ). Furthermore if all the particles are different, the 
projection of A^ to the complement b is required to vanish: A^\ b = 0. However, if the 
particles are identical the restriction of A^ to the complement b is non-trivial and then 
describes the additional exchange interactions between the identical particles. Thus the 
gauge interactions of non-identical particles are of Abelian character whereas the identi- 
cal particles do feel an additional force (exchange force) which requires to work with a 
non-Abelian gauge field theory. In the latter case, the theory becomes also non-linear 
which of course renders more difficult the construction of an adequate energy functional. 

In view of such a logical constellation, it seemed meaningful to first set up an energy 
functional Em for a system of different particles where the associated gauge field theory 
is Abelian (A^ = Ay\ ). The corresponding results were applied in their non-relativistic 
form to the positronium spectrum and yielded physically reasonable results [121 H3] , albeit 
the magnetic interactions have been neglected. Furthermore, in order to have at least an 
approximative energy functional for identical particles (A^\ b ^ 0) the gauge field theory 
was linearized [T3]. But clearly, such a linearization procedure must necessarily fail to be 
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rigorously gauge-invariant since for this property the non-linear terms of the non-Abelian 
gauge theory are indispensable (see below). Therefore the wanted gauge-invariant energy 
functional Em for a three-particle system with two identical particles (e.g. electrons) and 
one different particle (e.g. positron or proton) is now constructed in the present paper. 
As a demonstration for the obtained results one considers again the case of positronium, 
but now with emphasis on the well-known ortho/para degeneracy and its breaking [TO] . 
Indeed, this is a crucial test situation for the constructed RST functional E^ because 
here one can explicitly display the conditions (or approximations, resp.) by which there 
first arises the expected degeneracy which then afterwards is broken by the magnetic (and 
anisotropic) effects. 

Unfortunately, exact solutions of the RST eigenvalue problem are not yet available 
(neither analytically nor numerically) so that the subsequent investigation of the or- 
tho/para degeneracy must necessarily be of a more qualitative character. But despite 
the use of very crude approximation techniques, the obtained positronium spectrum is 
amazingly close to the conventional predictions and the observations (see fig. 4). This 
should provide sufficient motivation for an extensive elaboration of RST in order to test 
its predictive power and usefulness to higher accuracy in a wider range of physical appli- 
cations. 

The scope of the present paper is restricted to the construction of a gauge-invariant 
energy functional {Em) and is organized as follows: 

A. Gauge Structure 

The intention of setting up a gauge-invariant energy functional Em requires to first 
define the general gauge structure of the theory. To this end, the basic features of RST 
are briefly sketched in Sect. II. The emphasis lies here on the fact that the very field 
equations of RST do automatically imply certain local conservation laws, e.g. that of total 
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charge 



v% = o , 



(1.4) 



or that of total energy-momentum 



™ (T)rp 
V ± i 



. 



(1.5) 



Indeed, these conservation laws play an important role for the present problems. 
Namely they both concern the question whether or not any entangled particle can be 
thought to carry its own individual properties, e.g. charge and energy-momentum, or 
whether perhaps only the whole system of entangled particles (quite similarly as one non- 
entangled particle) can be considered to be the carrier of well-defined electric charge and 
energy-momentum? (The latter question will be answered in the positive). Here, the 
general belief says now that those well-defined and observable properties must be gauge- 
invariant, whatever the corresponding gauge group should look like. Therefore the gauge 
structure of RST in general, and especially also for the bound stationary systems, must 
first be clarified in detail. Of course, the wanted energy functional Em must be required 
to be invariant under the action of the ultimately established gauge group. 

In Sect. Ill it is made clear that the proper gauge group of electromagnetically in- 
teracting particles is the magnetic subgroup of U(l) x U(l) x U(l) x ... x U(l); this 
means that the electric potentials ^Aq{t) are left invariant and only the magnetic poten- 
tials A a {f) become transformed, namely in the usual inhomogeneous way (see equations 
fllll.20ap - fllll.20cp below). An important point refers here to the exchange potential 
which mediates the exchange forces among the identical particles; indeed, this potential 
acts like a tensor object and thus transforms homogeneously as is the case with the wave 
functions, too (see equations (1III.6P and (1III.9P below). 

B. RST Entanglement 

Quite generally speaking, the phenomenon of entanglement shows up in the occurrence 
of the so-called exchange effects. From the physical viewpoint, the main difference be- 
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tween the conventional theory and the present RST refers to the status of these exchange 
interactions: The conventional theory incorporates this type of interactions among iden- 
tical particles via the postulate that the N-particle wave functions should be symmetric 
for bosons and antisymmetric for fermions (or one postulates the anticommutation rela- 
tions for fermionic field operators and the commutation relations for bosonic operators, 
resp.). In the first line, such commutation postulates are of kinematical character and 
their dynamical consequences (e.g. exchange effects) are then studied afterwards, mostly 
in some approximative way. However, in RST the exchange effects are treated on the 
same footing as the electromagnetic interactions! This is briefly sketched in Sect. III. 

Such an ontological similarity of the electromagnetic and exchange forces implies that 
the latter are to be described by some four-potential as the exchange counterpart of 
the electromagnetic potentials A^, and both types of four-potentials together combine 
to the bundle connection A^, see equation (111.21) below. But clearly, such an equivalent 
treatment of the electromagnetic and exchange forces entails the common emergence of 
electromagnetic and exchange fields at all the relevant places of theory; for instance, the 
source and curl equations of the Maxwellian type for the electric (E) and magnetic (H) 
fields must be complemented by a source equation for the electric exchange field X and a 
curl equation for the magnetic exchange field Y, see equations (lIII.31ap - (lIII.31bj) below. 

But despite this formal similarity of both kinds of forces there are also important 
differences, especially concerning the gauge-fixing conditions. For the inhomogeneously 
transforming three- vector potentials A(r) it is still adequate to resort to the usual Coulomb 
gauge (V* A = 0); but for the homogeneously transforming tensor objects the analo- 
gous relation must be gauge- invariant, see equation ( 1III.33I) below, in order to not spoil 
the gauge covariance of the RST eigenvalue problem. 

In principle, the latter problem consists of the coupled set of Dirac equations for the 
matter fields, see equations ( 1II.6al) - (1II.6cl) below, and of the (non-Abelian) Maxwell equa- 
tions (see (111.81) below). But if one splits up these non-Abelian gauge field equations 
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into its electromagnetic and exchange parts, one observes the occurrence of an entangle- 
ment current which is built up by the complex-valued exchange potential like a 
Klein-Gordon type of current, see equation (IIH.43P below. 

The phenomenon of entanglement for two identical particles shows up also through 
the circumstance that the entanglement current Z M modifies the original electromagnetic 
currents of the particles in such a way as if some portion (5 e ) of charge (e), i.e. the 
entanglement charge, would have been exchanged between both identical particles, see 
equations ( 1III.66l) -( 1III.67bl) below. However this "virtual" charge exchange among the 
identical particles occurs in such a way that the total charge (e s ) of all the entangled 
particles together is unmodified and equals the product of the elementary charge times 
the number of identical particles, see equation (1III.63P below. In this way, charge conser- 
vation holds for the system of entangled particles as a whole, albeit not for the individual 
entangled particles themselves! 

C. Exact Energy Functional 

It is well known that the non-linear character of the non-Abelian gauge field theories 
is an implication of just that non-Abelian structure of the gauge group; and therefore 
the transition to the linearized theory must necessarily spoil the desired construction of 
a gauge-invariant energy functional, see ref. [13] . Nevertheless, the preliminary results 
of the linearized theory are very helpful for constructing now in Sect. IV the desired 
gauge-invariant form of the exact energy functional. Namely, that truly invariant form 
is essentially the same as for the approximative result of the linearized theory; however 
(i) the ordinary derivatives of the RST fields are now replaced by the gauge-covariant 
derivatives and (ii) there do emerge additional exchange terms (which vanish again for 
non-identical particles because these are unable to feel the exchange forces). 

But apart from this, the general structure of the gauge-invariant functional E^] (see 
equations ( 1IV.20I) and ( 1IV.23I) below) resembles very much its non-gauge-invariant pre- 
decessor [13J : there are the contributions with an immediate physical meaning (denoted 
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by , cf. (II V. 230 ) and these physical contributions are to be complemented by various 
constraints, e.g. the normalization constraints for the wave functions. The other con- 
straints do refer to the Poisson and exchange identities which, roughly speaking, ensure 
the equality of the gauge field energies and their "mass equivalents", i.e. the interaction 
energies of the coupled matter and gauge fields. This is immediately clear for non-identical 
particles where all the exchange terms are zero as a consequence of the vanishing of the 
exchange potential B^. But the physical meaning of the relationship between the gauge 
field energies and their mass equivalents becomes more intricate for identical particles 
where some of the constraint terms (i.e. those of the type Nq) are subjected to the varia- 
tional process while others (i.e. those of the type uq) are to be treated as true constants, 
see equations fjlVTl) . flTOl) . flIV.111) and flIV.141) below. Thus it seems that the true 
constants (of the type uq) are to be conceived as a kind of quantum numbers for the 
exchange interactions. 

However, under strict regard of all these constraints, the variational process (SE^] = 0) 
yields again both the gauge- covariant mass eigenvalue equations for the matter fields and 
the gauge field equations. This pleasant result does refer to both the relativistic situation 
and its non-relativistic approximation, where fortunately the gauge covariance survives 
the passage to the non-relativistic limit. A manifestly gauge-invariant form of the energy 
functional Em is presented in Appendix A. But observe here again, that the notion 
of "gauge-invariance" does refer here exclusively to the gauge transformations of the 
magnetic type! 

D. Ortho- and Para-Positronium 

As a plausibility test for the constructed energy functional E^ one will select some 
simple but sufficiently non-trivial system (here: positronium). As is well known, this 
simple two-particle system does occur in two forms: ortho- and para-positronium; see, 
e.g., ref.s [HH [TTj . Since these two forms have slightly different binding energies, the 
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corresponding ortho/para level splitting appears to provide an ideal test situation for the 
present energy functional which, if physically relevant, should be able to account for 
just this kind of splitting. 

However, a closer inspection of this ortho/para dichotomy in Sect. V reveals its 
complicated appearance within the framework of RST. Here, it turns out that the en- 
ergy difference between the ortho- and para-configurations is not only a consequence of 
the magnetic interactions between both positronium constituents but is also caused by 
their different anisotropic charge distributions. This difference of the charge distributions 
is itself an implication of the different angular momenta of the ortho- and para-states: 
whereas the para-states must always have vanishing total angular momentum, cf. flV.46p . 
the ortho-states may carry an integer number of angular momentum quanta, cf . (1VI.119|) . 
Here it seems plausible that the charge distribution of the para-states is concentrated 
in the immediate vicinity of the symmetry axis, which then implies lower angular mo- 
mentum; on the other hand, the charge distribution of the ortho-states splits into two 
subforms (fig. 5) which are both zero on the symmetry axis and adopt their maximal 
value aside of this axis. Clearly such a larger distance of the rotating matter from the 
axis of rotation will entail also a larger angular momentum. 

But the crucial point is now that the difference of the anisotropic charge distributions 
entails a corresponding difference of the electrostatic interaction energies and these are nor- 
mally thousand times greater than the experimentally verified difference of the magnetic 
interaction energies (see the discussion of this point in connection with the ortho/para 
degeneracy ( IV. 13j) below) . This hierarchy problem, concerning the strength of the electric 
and magnetic interaction, forces one to reconsider the ortho/para degeneracy and to think 
about some approach which at least approximately attributes the same electrostatic inter- 
action energy to the different charge distributions of the ortho- and para-states: this is the 
spherically symmetric approximation where the electric interaction potential ( 6 / p l4o(r, •&) 
is replaced by a spherically symmetric potential ' 6 / p l4o(r) : ( 6 / p l4o(r, d) =^ ' b / p L4o(r) (the 
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wave functions if> a (f) are still admitted to be non-spherically symmetric). And indeed, 
this assumption of an SO (3) symmetric interaction potential ' b / p l4o(r) lets coincide the 
eigenvalue problems for the ortho- and para-states to one and the same form, cf. the para- 
system f lV1.23al) - m.23bp . (IVO8|) - flVT^0l) and its ortho-counterpart flV1.115al) - (IVI.115bp . 

The common non-relativistic spectrum of both eigenvalue problems is sketched in 
fig. 2. For a rough quantitative estimate of the energy levels a two-parameter variational 
solution is chosen (see equation ( 1 VI. 5511 ). This choice predicts the RST energy levels with 
deviations not more than some 10% in comparison to the conventional predictions, see 
fig.s 4a-4b and the table on p. 11071 

A difference between both level systems refers to the quantum numbers i-p and Iq for 
orbital angular momentum: while for para-positronium the associated quantum number 
is £j> = 0,1,2,3,4.. ., the range of the corresponding ortho-number £q is to = 2, 4, 6, 8, . . . 
(see appendix E). It is true, the range of the ortho- number £q is half the range of the para- 
number i-p, but for any possible value of the ortho-number Iq there exist two different 
charge distributions, see fig. 5, which are supposed to cause different binding energies E^. 
However, a discussion of the corresponding breaking of the ortho/para degeneracy must 
be deferred until more exact solutions of the RST eigenvalue problems can be elaborated, 
especially with regard to the anisotropic and magnetic effects (for a first estimate of 
the magnitude of the latter effects see ref. [15]). Concerning the relative magnitude 
of both effects, the anisotropic electric correction is found to be of order 1 [eV] (i.e. 
roughly 10% of the binding energy [12]), whereas the magnetic corrections amount to 
roughly 10" 3 [eV] [E]. 
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II Relativistic Schrodinger Theory 

All those important symmetry principles such as Lorentz and gauge invariance, being 
considered as indispensable for any successful field theory of quantum matter, are encoded 
manifestly in RST through the general form of the basic field equations. The subsequent 
brief account of those field equations will be taken as the point of departure for the 
construction of the desired gauge- invariant energy functional E^. (For a more detailed 
discussion of the RST field equations, see the preceding papers p~2], HH] . 

1. Matter Fields 

First, the key equation for the N-particle wave function \I/ is chosen to be the Dirac 
equation 

ihcT^V^ = Mc 2 ^ . (Ill) 

Here, the derivative (T>) of the wave function \1/ has already been specified by equation (11.21) 
where the gauge potential (bundle connection) A^ takes its values in the Lie algebra U(N) 
of the unitary group U(N), provided all the particles are identical. When some (N — N', 
say) of the iV particles are non-identical (e.g. differing by their masses or charges), then the 
structure group U(N) of the fibre bundles becomes reduced to some subgroup U(N,N'). 
In order to present a concrete example, consider a three-particle system (N = 3) where 
two particles are identical (N f = 2) and one particle is different (think, e.g., of a positron 
or proton and two electrons). For such a situation, the original nine-dimensional structure 
group t/(3) for three identical particles is reduced to the product group {7(1) x {7(2) (for 
the reduction process, see ref. [H]). Accordingly, the gauge potential A^ for such a three- 
particle system can be decomposed into a five- dimensional Lie algebra basis {r a , x, X.} 

5 3 
a=l a=l 

where the electromagnetic generators r a (a = 1,2,3) do commute and are anti-Hermitian 
(f a = —r a ); and the exchange generators x,X are taken to be (pseudo-) Hermitian con- 
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jugates. The anti-Hermiticy of the gauge potential Ay 



-Ay) requires then that the 



electromagnetic 'potentials A° 



1,2,3) are real- valued whereas the exchange poten- 



tial is admitted to be complex-valued. 

Returning to the iV-particle Dirac equation (IH.ip there remain two elements to be 
explained, i.e. the total velocity operator F M and the mass operator M.. The first one of 
these objects is the direct sum of the ordinary Dirac matrices 7^, e.g. for the mentioned 
three-particle system 

Ty = (- 1 y)®ly ®ly. (II.3) 

Furthermore, the mass operator M. can be adopted to be diagonal in a suitably cho- 
sen U(3) gauge 

/ M P 
M M e 
M e 

where M p (M e ) is the rest mass of the positively (negatively) charged particles. But 
contrary to the conventional {probabilistic) theory, the iV-particle wave function \& in 
RST (as a fluid- dynamic theory) is not built up by suitable tensor products of single- 
particle wave functions but rather appears as the Whitney sum of single-particle bundle 
sections ip a , e.g. for the considered three-particle system 



V 



\ 



7 



(II.4) 



(II.5) 



As a consequence, the matter equation (III. lft may be resolved in component form and 
then yields the following coupled system of ordinary Dirac equations (for the considered 
three-particle system) 



ihcrfDpfa = -M p c 2 V>x 
ihcrf^D^ = M e c 2 V>2 
iturfD^ = M e c 2 ip 3 



(II.6a) 
(IL6b) 
(II.6c) 
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where the gauge-covariant derivatives of the wave function components ip a (a = 1,2,3) 
are given by 

D^ x = d^ x - % (A\ + A\) fa (II.7a) 
D^ 2 = d^ 2 - i (A\ + A\) fa - iB^fa (II.7b) 
D^fa = d^fa - i (A\ + A\) fa - iB^fa ■ (IL7c) 

Observe here that the exchange potential establishes a direct coupling of the wave 
functions fa and ^3 of the identical particles (111.7bl) - (jll.7cp . whereas the non-identical 



particle flll.Taj) does couple only indirectly (i.e. via the electromagnetic potentials A 



2 



and A 3 ^) to both identical particles. This effect of RST entanglement of the identical 
particles must of course be discussed in more detail below. 

2. Gauge Fields 

The role of the Dirac equation (111.11) for the matter field \I/ is played for the gauge 
field An by the (non-Abelian) Maxwell equations 

V^ u = -Ama s J v (II.8) 
e 2 

« s = ~r 

nc 

Here the bundle curvature (field strength) is defined as usual in terms of the connec- 



tion An as 



J> = V^A- V v A lt + [A^Au] . (II.9) 



The curvature may be decomposed for the presently considered three-particle system quite 
similarly to the connection A^ ( 111.21) as 

5 3 

^> = Yl = Yl F<1 ^ Ta + G ^ - G»vX (H.10) 

a=l a=l 
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so that the curvature components F" 1 ^, G^ appear in terms of the connection compo- 
nents A a /1 , in the following way 

fV = V^- v ^ ( ILlla ) 

F\ v = V^A\ - V V A\ + i \bJ* v - B u B^j (Il.llb) 
F\ u = V^A 3 U - V V A\ - i \bJb v - B V BA (II. 11c) 

G>, = V„B„ - V,5 M + % (A\ - A\) B v - i (A\ - A\) B„ . (II. lid) 



Similarly, the Maxwell equations ( 111.81) read in component form 



^F\ u 




(II. 12a) 


B^G^ 


= i7ra s j 2 v 


(II.12b) 


B^G^ 


= Aira e j 3 v 


(II.12c) 


F\ v ) B» 


* 

= 47ra s h v . 


(II.12d) 



Here the current operator f ill. 8f) has been decomposed according to 

^ = *Z^V« = 1 fe A r « + - Kx) , (H.13) 

a=l \a=l J 

where the required charge conservation law (11.41) emerges in operator form as 



T>^ = . (11.14) 

Indeed this is an automatique implication of the Maxwell equations ( 111.81) since the cur- 
vature T^v must obey the bundle identity 

2WJ> = . (11.15) 

Surely, the emergence of the charge conservation ( III. 141) as an immediate consequence 
of the field equations represents a pleasant feature of RST and therefore deserves some 
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closer inspection. First observe here that the current components (Maxwell currents) 
j a „ = {j a fl ,h^,h^ l },a = 1...5,a = 1...3, play the role of the sources for the gauge 
fields A a p but cannot be immediately identified with the well-known Dirac currents k afl 



Ku =F Ipal^a ■ (11.16) 



The precise relationship is 



J% = feiAi = $i7/iV>i (II. 17a) 

j'V = -Aty, = -^27/^^2 (II. 17b) 

J'V = = "^37^3 (II. 17c) 

j\ = h, = ^37^2 (H.17d) 

J'V = -fyi =4= -^27^3 • (H-17e) 

With respect to these identifications, the operator form of the charge conservation law ( 111.1411 
transcribes to the corresponding component form in the following way: 

= (II.18a) 

VAfy - i (b*K - b4A = (II. 18b) 



V^ 3M + i IB"^ - B»hA = (IL18c) 
- i {A\ - A\) W + ib» (k 2lI - k 3ll ) = . (II.18d) 

A further pleasant feature of the RST dynamics concerns the automatic conservation 
of energy-momentum (II. 5\i . Indeed, this works as follows: defining the energy-momentum 
density ^T^ u carried by the Dirac matter field \I/ through 

{D %» = ^ [W M (Z>„tt) - ipjf) + VT V (V^) - (Vfi) V v *] (11.19) 

and, analogously, the energy-momentum density concentrated in the gauge field 

through [15] 

(G) T„ = ^- K aP (f% x F? u x - \F\ X F^ (11.20) 
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then the total density 

(T) T„ = V% u + ^T, u (11.21) 

actually obeys the required (local) law (1151) of energy momentum conservation, provided 
that only both sub-dynamic equations for the matter subsystem ( lll.lj) and for the gauge 



field subsystem ( 111.81) are satisfied! Clearly with this splitting ( III. 211) of the total den- 



sity ^'Tjy,, the total energy E T [T3] also splits up into two constituents, i.e. 

E T = E D + E G , (11.22) 

with the obvious definitions of the partial energies (for stationary bound systems) 

E D = J d 3 f^T 00 (r) (II.23a) 

Eq — J d 3 f^T 00 (r) . (II.23b) 

This aims now just at the central point of the present paper: namely the question, to 
what extent can both the stationary matter equations (III.6ap - (III.6cl) and the (non-Abelian) 
Maxwell equations ( 1II.12al) -( lII.12dl) be conceived as the extremal equations due to the to- 
tal field energy E^ (1II.22j) . to be understood as an energy functional Eml Or in other 
words, do the stationary solutions of the coupled matter and gauge field systems (111.6ap - 



f lll.6cl) plus (1II.12al) -( lII.12dl) actually extremalize the total energy functional E^ (1II.22I) ? 
For non-identical particles (~~* = 0), this question has been answered in the positive 
so that it is not necessary to modify the form of that functional Em; but for identi- 
cal particles (~» ^ 0) it could be shown (albeit only in the linear approximation) 
that the orig inal definition (ITL2^ -( tTL23bl) of the energy functional E^] must be slightly 
modified [13] . But of course, one wishes to have the question clarified for the original non- 
Abelian (and therefore non-linear) theory and this is just what we intend to do by means 
of the subsequent elaboration. Observe here the fact that the complication, being caused 
by the presence of some identical particles, is due to the non-vanishing exchange poten- 
tial B^ which renders the theory non-Abelian and therefore non-linear! However, this 
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specific feature of the theory is indispensable for the gauge-covariant description of the 
exchange interactions which are to be conceived as the RST equivalent of the conventional 
phenomenon of entanglement. 
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Ill Gauge Invariance and Entanglement 



It should be obvious that there must exist an intimate relationship between the struc- 
ture group of the theory and the exchange effects which can occur in such a theory of 
the RST type. The point here is that the interaction potential A^ takes its values in the 
Lie algebra of the structure group U(N, N'), cf. ( 111.21) for the presently considered three- 
particle system, where this structure group results by reduction of the original structure 
group U(N) for N Actively identical particles. This process of reduction my be attributed 
to the fact that only the permutation of identical particles leaves invariant the physical 
properties of the N particle system where (in the true spirit of the gauge theories) the 
discrete set of permutations or relabelings, resp., within the subset of identical parti- 
cles is generalized to a continuous group U(N'). Therefore, avoiding the permutations 
of non-identical particles, the reduced structure group of the bundle arrangement must 
contain U(N') as a factor group. If we further assume that all the residual N — N' 
non-identical particles do feel exclusively the ordinary (i.e. Abelian) electromagnetic in- 
teractions, we come to the conclusion that the reduced structure group should be the 
product U(N,N') = U(l) x U(l) x . . . x U(l) x U(N'), i.e. for our presently considered 
three-particle system U(l) x £7(2). Indeed, the permutative decoupling of the first (posi- 
tively charged) particle a = 1 from the other two (negatively charged) particles a = 2, 3 
become manifest in both relationships ( 1II.llaj) -( 1II.lldl) between the curvature and the 
connection components and becomes manifest also in the component form of the Maxwell 
equations flII.12aD - flII.12dl) : moreover an analogous observation does apply also to the 
source equations flII.18al) -( lII.18dl) . Notice however, that the total fields of the identical 
particles obey an ordinary Abelian structure! 

In order to realize this more clearly, add up both equations (IH.llbl) and (III.llcl) in 
order to find 

F\ u = V,A S U - V U A% (III.l) 
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where the total fields due to both identical particles are defined through 

A% = A\ + A\ (III.2a) 
F% V = F\ V + F\ V . (III.2b) 

Indeed, these relationships for the total fields obviously avoid any reference to the exchange 
potential B^. A similar conclusion does hold also for both the Maxwell equations flII.12bj) - 
( III. 1 2c[) whose sum appears as 

= 4vr« s = 4ira s f v , (III.3) 

or similarly for the sum of equations (1II.18bj) plus (III. 18c[) 



(k 2 ^ + k 3ll ) = Wk sft = (III.4) 

which is now a proper conservation equation for the total Dirac current k Sfl of identical 
particles. Thus the subsystem of identical particles acts in an Abelian way like one (albeit 
compound) particle relative to the outside world. 

1. Gauge Structure 

The reduction of the original structure group U(N) to U (N, N') does not mean that 
we adopt the latter structure group as our proper gauge group because the general RST 
philosophy says that any particle of the system must occupy a well-defined one-particle 
state. This requirement excludes the state mixing of (both the identical and non-identical) 
particles through the gauge action of a general U(N, N') or U(N') element. This becomes 
especially clear by considering the stationary bound systems which are to be described by 
a set of one-particle wave functions i/) a (r, t) appearing as products of a time-independent 
factor if} a (f) (a — 1, . . . N) times the usual exponential time factor, e.g. for the presently 
considered three-particle system (a = 1,2,3) 

Va(r,t) = exp f ^—tj ■ r(> a (f) . (III.5) 
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Indeed it is easy to see that an 17(2) rotation of the wave functions ip a (f, t) of both identical 
particles (a = 2, 3) with different mass eigenvalues M a would in general generate a mixed 
time dependence violating the manifest stationarity. 

From this reason, we dismiss the structure group U(N, N') as the gauge group and in 
place of it take the simple product 1/(1) x 17(1) x U(l) . . . U (1) as our proper gauge group. 
Thus a change of gauge acts upon the one-particle wave functions in the following form: 

^ a (x) tf a (x) = e-^W ■ V«(a:) (III.6) 

where the group parameters a a (a = 1,2,... N) are real- valued functions over space-time. 
As is well-known, the gauge potentials A a ^ do not transform homogeneously as do the 
wave functions ip a (1III.6P but rather obey some inhomogeneous transformation law, e.g. 



for our three-particle system: 



(\ => (\ = % - d,a, (III.7a) 
=>■ ( n k; = C% - d,a 2 (III. 7b) 

{ll X m K = (ll X - d,a 3 , (III.7c) 



where the modified electromagnetic potentials are defined through 



A\ + A\ (III.8a) 
A\ + A\ (III.8b) 
A\ + A\ . (111.8c) 



However in contrast to these electromagnetic potentials, the exchange potential un- 
dergoes a homogeneous transformation law, i.e. 

S^B; = e-^-^ ■ . (III.9) 

For instance, by means of this gauge structure it becomes possible to gauge off the time 
dependence of the stationary wave functions ip a (r, t) flIII.5p . namely by simply putting for 
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the gauge parameters a a (t) 

a a = -^t (111.10) 



which then leaves the spatial components of the modified four potentials ( 1III.7al) -( lIII.7cf) 
unchanged but recasts the time components to 

<% =► <% = <% + ^ (Hl.lla) 

("k =► v% = ( n k + (ni.iib) 

("% ^ ("X = (H % + ^ • (ULllc) 

We will not make use of such a purely temporal gauge transformation for the following dis- 
cussions but one point deserves here some attention. Subtracting the equations 1 lcH 
minus flHI.llbj) yields a special combination 21m of potentials 

2l M = {U % - {1 % = % ~ — (HI.12) 

(sto = ( % - (3 k ) 

where the exchange radius <2m is evidently defined by 

«m = . H M , ■ (HI.13) 
(M 2 - M 3 )c 

Indeed, this potential 21m (IIII.12I) will frequently occur in connection with the wanted 
energy functional Em. The reason here is that, for the stationary bound systems, the RST 
field equations must be time- independent, and this time independence may be conceived 
to be the result of a temporal gauge transformation fllll. 10p - (IIII. 1 lc|) . 

Finally, observe also that the time dependence of the exchange potential can easily 
be determined from its gauge behavior: assuming that the new exchange potential B 1 
(1III.9P is time-independent, i.e. B' = B^ir), it follows immediately from the equations 



flIII.9l) -f lTlI.10l) that the original time-dependent B^r, t) must look as follows 

B^t) = exp(-z— ) ■ B,(f) . (111.14) 
au 
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This is the temporal countercurrent of the exchange current h^f, t) ( 1II.17dl) -( lII.17el) whose 
time behavior for the stationary states is given by 

h^f, t) = exp ( i— ) • h^r) . (111.15) 



On the other hand, the Dirac currents k apu (1II.16P will be found to be time-independent 
for the stationary states, i.e. 

£V = { (a) Mr);-4(r)) , (HI.16) 



and since these are the sources of the electromagnetic potentials A a ll , cf. (1II.12al) - (1II.12cl) 



one will find the latter objects also being time-independent for the stationary states 

A% = {^A (f);-A a (r} . (in. 17) 

Clearly, the Dirac currents k a {r) 1 6H must have vanishing divergence for such a situ- 
ation 

V^ = V-4(r)=0, (111.18) 
and if this is combined with the former source equations fill. 18bj) - fill. 18cl) one concludes 



that the Lorentz scalar products B^h^ must be real numbers: 

B^ = B*h . (111.19) 

However, for the subsequent discussion we do not wish to modify the standard time 
behavior of the wave functions ip a (r,t) flIH.5p and of the exchange potentials B^ir^t) 



( lIII.14p . Therefore we restrict ourselves to the purely magnetic gauge transformations 
where all the gauge parameters a a are time-independent (a a = a a (f)). In this case, the 
magnetic three-vector potentials Ai,Aa,Aai f lHI.7aP -f lIH.7cl) transform as follows 

A x =>- A[ = Ai + Vai (III. 20a) 

A u A n = Aa + Va 2 (III.20b) 
-4*iii => An = An + V« 3 , (III.20c) 
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whereas the electric potentials ^Mo, ^ u Aq, ^ m Ao remain invariant. Thus the gauge invari- 
ance of the wanted energy functional will refer exclusively to those gauge transformations 
of the magnetic kind! 



2. Mass Eigenvalue Equations 

This restriction to the magnetic gauge transformations is also important for the gauge 
covariance of the mass eigenvalue equations which are nothing else than the stationary 
form of the original matter field equations 0II.6ap - flII.6cp . In order to set up this stationary 



form, one first splits up the Lorentz covariant derivatives flll.7ap - flll.7cp into their time and 



space components and combines this with the fact that the Dirac four-spinor fields tp a (r) 
can be conceived as the Whitney sum of two-component Pauli spinors ^<p±(r), i.e. 

Mr) = ( "Wr) © (111.21) 
'a = 1,2,3^ 

This then lets appear the stationary version of the coupled system of Dirac equations 
f lll.6ap - flll.6cp in the following gauge- covariant Pauli form 

±M P - Mi 



ia-fv + iA^j <»<p±(f) + ( % ■ (1 VtW = \ c ' (1 W*0 (III.22a) 
,cr • (V + zl n ) (2 V±(r) + (I % • (2 V T (r) + B • ( 3 V T (f) -B-a < 3 V±(r 

M ^ M ^.^r) (III.22b) 



M3±M -c^(r,. (111.22c) 



ft 

Observe here again that the gauge covariance of this coupled mass-eigenvalue sys- 
tem refers exclusively to the magnetic gauge transformations (lIII.20ap - flIII.20cp whereas 



the electric potentials QAq, ^ u Aq, ( 1u Aq remain invariant, in contrast to the case of the 
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temporal gauge transformations ( 1III.llal) -( 1III.llcl) . Notice also that for the mass eigen- 
value equations ( 1III.22al) -( 1III.22cl) we preferred the use of the original electric poten- 
tials (I k , (H % hi order that the mass eigenvalues M a (a = 1, 2, 3) become explicitly 
displayed. But observe that the latter do appear in such a specific way that they in 
principle could be absorbed also into the transformed potentials ( lIII.llal) -( lIII.llc l ) after 
performing the gauge transformation (1III.7al) - (1III.7cl) in combination with ( 1III.10I) . 

3. Gauge Field Equations 



Concerning the magnetic gauge covariance, the situation with the gauge field equations 
is somewhat different from the matter field equations. Namely, building up the Lorentz 
covariant objects F a ^ (or G^, resp.) by the usual electric and magnetic three-vector 
fields E a , H a (or X and Y, resp.) one finds the three- vector versions of the former four- 
vector relations (1II.llal) - (1II.lldl) to look as follows, namely for the electric fields 



e\ = -v (1 k 

E 2 = -V ( % - i 
E 3 = -V (3 k + i 



B B — B B 

' t * 
B B — BqB 



(III.23a) 
(III.23b) 

(III.23c) 



and similarly for the magnetic fields 



H x = V x A l 

* 

H 2 = V x A 2 - iB x B 
H 3 = V x A 3 + iB x B . 



(III.24a) 
(III.24b) 
(III.24c) 



Obviously these electromagnetic field strengths E a {r) and H a (f) are invariant under the 
magnetic gauge group U(l) x U(l) x U(l) and consequently one expects this invariance 
to hold also for the corresponding field equations (see below). In contrast to this, the 
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exchange field strengths X and Y of both the electric and magnetic type are sensitive 
against a magnetic gauge transformation 

X(r) = -V'B (r) - i Ql M B(r) (III.25a) 
Y(r) = V' x B(r) . (III.25b) 

Here 21m is the electric potential difference introduced by equation (1III.12I) and V is the 
gauge-covariant generalization of the ordinary gradient operator V, i.e. 

V' = V - i% (111.26) 
(2l = A 2 - 1 3 ) • 

Evidently, both exchange field strengths X and Y obey the homogeneous (tensorial) 
gauge transformation law 

X => X' = e ~ i{a2 - az) ■ X (III.27a) 

Y => Y' = e- i{a2 ~ a3) ■ Y , (III.27b) 

and consequently one expects gauge-covariant field equations for both objects X, Y in 
contrast to the situation with the gauge-invariant electromagnetic field strengths E a , H a . 
Let us remark here that the check of the claimed transformation laws ( 1III.27al) -( lIH.27bl) 



is a nice exercise where the vector potential difference 21 can be shown to change like 

a=>a' = 2t-V(a 2 -a 3 ) . (111.28) 

Of course, the field equations of those three-vector objects must again be deduced 
from their Lorentz-covariant predecessors (1II.12al) - (1II.12dl) . Here, the gauge invariance of 
the electromagnetic part is again manifest, i.e. one has for the electric part 

V • E x = 4vra s • (1) j (III. 29a) 

B • X - B • X ] = 47ra s • (2) j (III.29b) 

B-X-B'X) = 4vra s • (3) j , (III.29c) 
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and similarly for the magnetic part 

V x E x = 47ra s Ji (III.30a) 

B X - B X + B xY - B xY \ = Ana s j 2 (III.30b) 

/ * * * * \ 

VxJJ 3 + i 5 X - B X + B xY - B xY) = 4vra s j 3 . (III.30c) 



(For the link of the Maxwell currents {j a M } = {^'o(^0; — Ja(^*)} to the Dirac currents 
{ka[i} = {^ko (r) ; — k a (r)} see equations (1II.17al) - (1II.17el) ). However the corresponding 
equations for the exchange field strengths display the expected gauge covariance mani- 
festly through the emergence of that gauge- covariant derivative V 26[) 



V' • X + B = 4ira s h (III.31a) 
V' x Y - M M X + iB € - iij x B = 4vra s h . (III.31b) 

Here the three-vectors S and are defined through the differences of the electric and 
magnetic fields due to the identical particles (a = 2, 3), i.e. 

i = E 2 -E 3 (III.32a) 
f> = H 2 - H 3 . (III.32b) 



Summarizing, the matter equations flIII.22aP - flIII.22cP and the gauge field equations 



( lIII.29aP -f lIII.31bp represent the essential RST field equations for the stationary bound 
systems; and the task to be tackled now is to construct an energy functional Em in such 
a way that its extremal equations turn out to be identical to just this coupled set of matter 
and gauge field equations! However before this can be done there remains one point to 
be settled which refers to the different transformation behavior of the electromagnetic 
potentials A" 1 ^ and the exchange potential B^, cf. f lIH.7aP -f lIH.7cP vs. (1III.9p . 
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4- Linearization and Gauge Invariance 

Both the electromagnetic potentials A a ^ and the exchange potential appear as 
independent parts of the bundle connection ( 1II.2j) . but there arises an important 
difference when the structure group U(l) x £7(2) becomes reduced to the proper gauge 
group £7(1) x U{\) x U{1). Namely, through this reduction process both bundle con- 
stituents A a M and do acquire a rather different geometric meaning: On the one hand, 
the electromagnetic potentials A a ll preserve their property as connection components for 
the reduced bundle and therefore do transform inhomogeneously, cf. ( IIII.7al) - (lIII.7cl) . in 
the true spirit of the generalized equivalence principle for the gauge theories. But on the 
other hand, the exchange potential B^ acquires tensorial properties under the reduction 
of the structure group and therefore transforms homogeneously, cf. ( Till. 91) . Thus it is no 



longer possible to locally gauge off the exchange potential B^ as it is the case with the 
electromagnetic potentials A a ^. Clearly, such a change of the geometric meaning of B^ 
must have its consequences for the field equations to be obeyed by this object. 

Recall here that the field equations (1II.12al) -( lII.12cl) for the potentials A a ll are based 
upon the "curls" (III. 11 aft - fill, llcl) and such a type of equation is just sufficient for the 



connection components. Indeed, it is not viable to fix also the divergences of A a fl by 
means of some source equation since this would unduly restrict the set of possible gauge 
transformations of the reduced bundle connection. However, since B^ has lost its in- 
homogeneous transformation properties through the reduction of the structure group to 
the gauge group, this object becomes thereby fixed to a greater extent than is the case 
with the potentials A a ^. As a consequence there must exist, besides the "curl type" of 
field equation (III.12dj) . also a further equation for B^ which is of the "source type". We 
will satisfy now this requirement by selecting such an additional "source type" equation 
for B^, namely 

WB^ = , (111.33) 
which is nothing else than the vanishing of the invariant divergence (i.e. a kind of invariant 
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Lorentz gauge condition). The covariant derivative (D) is defined here as 



%B V = V^B V + i%B v (111.34) 
= A 2 M - A 3 ^ (111.35) 

so that the source relation (1III.33I) may be also rewritten as 

V'*B = -m M B (111.36) 

with the covariant gradient operator V being defined by equation (IHI.26I) and the po- 
tential difference 21m by flIII.12p . This is to be conceived not as an additional dynamical 
equation but rather as an auxilary constraint (in the sense of a gauge constraint) in order 
to ensure the existence and uniqueness of the wanted energy functional E^ (see below). 

The latter relation (1III.36H is well-suited in order to explicitly demonstrate how the 
gauge invariance of the theory may become spoiled through the process of linearization. 
Namely, the explicit form of ( 1III.36I) reads 

V'B-i%*B = -i( 2l - — ) B . (111.37) 



But linearizing the theory in the literally sense means that we omit the products of gauge 
potentials 21 • B and A B so that the residual divergence relation reads 

V'B = —B . (111.38) 

Indeed this truncated non-invariant relation was used for the linearized theory (see equa- 
tion (1III.18P of ref. [13]), but here it is revealed now as one of those places where the strict 



gauge invariance is spoiled through the linearization process. 

5. RST Entanglement 

Attention to the strange phenomenon of "entanglement" has been drawn for the first 
time by Schrodinger. He discussed this effect within the conceptual framework of the 
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conventional (i.e. probabilistic) quantum theory, where its observable consequences refer 
to the well-known exchange effects (e.g. shift of atomic energy levels by the "exchange 
energy" [10J). However, this can not mean that other (non-conventional) quantum ap- 
proaches could not take account of such exchange effects, too. Indeed, it has already 
been demonstrated in some preceding papers that those exchange effects can on principle 
also be described within the framework of RST. But whereas these demonstrations re- 
ferred mainly to concrete numerical examples, it seems now time to discuss the exchange 
phenomenon on a more general level. 

Here, the point of departure is the component form ( lII.12al) -( 1II.12dl) of the non-Abelian 
Maxwell equations (111.80 . It seems rather obvious that these equations suggest to intro- 
duce an entanglement vector G^, namely through 

G (t = -^-(B v G V(t -B v G vl ^\ . (111.39) 



By this arrangement, the Maxwell equations ( lII.12bl) - (1II.12cl) for both "entangled" parti- 
cles (a = 2, 3) adopt the following form 

V^F V + 4vra s G v = 47ra s j\ (III.40a) 
VF V - 4vr« s G v = 4vr« s j\ . (III.40b) 

This could be interpreted in such a way that the Maxwell currents j a (a = 2, 3) as the 
sources of the field strengths F a ^ u , become modified through some kind of "entangling cur- 
rent" (~ Gfj) and thus generate a modified electromagnetic potential A a M (a = 2, 3) as the 
solution of just those equations ( 1III.40al) - (1III.40bl) . This RST entanglement effect would 



then consist partly in the modification of the electromagnetic interactions (~ ^4%) of the 
entangled particles and partly in the emergence of a new additional kind of interaction 
being mediated directly through the exchange potential B^, see the covariant derivatives 
( 111.7bl) -( fll.7cj) . However observe here that the entanglement vector G^ formally drops out 



for generating the total field strength F s ^ 3 j) of the entangled particles so that the 
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influence of the entangled subsystem on the outside world (here the first particle, a = 1) 
works exclusively via the electromagnetic potential A s M ( 1III.2al) which coincides with the 
former potential ^'A^ (1III.8a[) . 

There arises now a certain ambiguity concerning the normalization of the wave function 
for an entangled subsystem; but this problem can be enlightened by a closer look at the 
properties of the entanglement vector ( 1III.39I) . For this purpose, introduce the skew- 
symmetric tensor ( 



-B VfJ ) through 



B 



Ana,, 



B^B U — B^By , 



(111.41) 



and then find by straightforward differentiation and observation of the vanishing diver- 
gence (1111.330 



G v 



B' 



KB, 



5' 



KB^ 



(111.42) 



From here it becomes obvious that the entanglement vector can be written as the 
sum of the divergence of a (four-dimensional) entanglement field strength B^ v and an 
entanglement current l v which is defined through 



47ra s 



B» CD u B, l ) - B^ 



KB, 



(111.43) 



Thus the splitting of G u (IIII.42I) reads 



G v = — (V^ + k 



(111.44) 



and if this is inserted in both equations ( lIII.40al) -( 1III.40bl) they do reappear in the form 

WF\„ = 4vm s (j\ + l u + VBj = 47ra s ■ ^f u (III.45a) 
V»F\ V = 4vr« s (j\ -l v - WB^) = 4vr« s • (III.45b) 

But here the entangled currents ^j 2 ^ and < -~- ) j 3 A1 have now vanishing divergence 

= = , (HI.46) 
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for by straightforward differentiation one finds 

VG, = -V% = i (b^ - ) (HI.47) 



which then just annihilates the sources of j 2 ^ (= — A^) and j 3 ^ (= —k 3fJi ), cf. (1II.18bj) 



f lII.18cl) . so that the source equations ( 1III.46I) can actually be true. The ordinary Maxwell 
equations flIII.45aP - flIII.45bp then say that the field strengths F 2 ^ u and F 3 ^ can be 



thought to be generated by those entangled and source-free currents ^j 2 ^ and < - _ - ) j 3 M . 
Such a situation would provide us with the possibility to postulate normalization condi- 
tions for the entangled wave functions ifraif) { a — 2, 3) of the following kind 

dS» = J dS" = -l. (111.48) 

(S) (S) 

Here, the integrals over some three-dimensional hypersurface (S) would then be indepen- 
dent of the choice of that hypersurface (S). However, we will not resort to this construction 
but rather will prefer the following normalization conditions: 

N D (a) = [ d 3 f {a) k (r) - 1 = (111.49) 



t=const. 



(a = 1,2,3) . 

The reason for this preference, where the independence of the 3-surface of integration is 
lost, will become clearer below in connection with the discussion of the gauge- invariant 
energy functional. 

Finally, it is also interesting to remark that the "external" field strength F s ^ flIII.2bj) 



of the entangled two-particle system may be complemented by an "internal" counter- 
part f\ u 

/V - F V - F \» ■ ( m -50 

The combination of both non-Abelian Maxwell equations f lIII.45al) -f lIII.45'b1) then yields 



V (/V - 2B„ V ) = 47ra s (f u - f u + 2l v ) . (111.51) 
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The comparison of this "internal" Maxwell equation to its "external" counterpart ( 1111.31 ) 
says that the internal gauge interactions do not obey the ordinary (Abelian) Maxwell 
scheme but take account of the non- Abelian character whose imprint on the internal gauge 
dynamics results in the additional objects B^ v and l v occurring in equation (lIII.5ip . 



6. Boundary Conditions and Entanglement Charge 

The preceding arguments were based upon a Lorentz-covariant presentation of the 
RST entanglement, but the wanted energy functional Em can naturally not fit into this 
scheme because the concept of energy is quite generally not a Lorentz-invariant object, 
even in the relativistic case. Therefore, in search of a gauge- invariant energy functional, 
it is more convenient to represent the preceding discussion also in three-vector notation. 
This will be helpful for building up the RST energy Et from its various constituents. 

First, decompose the tensor B^ v ( IIII.4ip into its "electric" and "magnetic" part ac- 
cording to 



W = i \^B B - B Bj (III.52a) 
V = i(^Bx B^j , (III.52b) 

so that the electric fields E a ( 1III.23bl) -( lIII.23cl) for the entangled particles appear now as 

E 2 = -V {2) A + W (III. 53a) 

E 3 = -V (3 k - W , (III.53b) 

and similarly for the magnetic fields H a ( lIII.24bl) -( 1III.24cT) 

H 2 = V x A 2 ~V (III. 54a) 

H 3 = V x A 3 + V . (IIL54b) 

Furthermore, it becomes obvious from equation ( 1III.34I) that, after bundle reduction, 
the homogeneously transforming exchange potential B^ acts like a charged matter field 
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with current ( 1III.43f) and thus feels the induced connection 2l M (= A 2 ^ — A 3 ^). The 
corresponding induced curvature follows from equations f lIII.53al) -( 1III.54bl) in terms of the 
original curvature components <£, 5} (1III.32a|) - (1III.32b[) and the entanglement fields W,V 
(llll.52al) - (IIII.52bp as 

- V2t = £-2W (111.55) 

and 

Vx2l = £) + 2V. (111.56) 

Conversely, the time component of the entanglement vector ( 1III.44I) reads in terms of 
the electric entanglement field strength W and density Iq 

V • w 

Go = k , 111.57 

47ra s 

and similarly for the space component 

— * —t 

where evidently V plays the part of a magnetic entanglement field strength and / is 
the entanglement three-current. Obviously, the latter two equations represent the three- 
vector reformulation of the Lorentz-covariant equation flIII.44p . In a similar way, the 
"internal" Maxwell equation (IIII.51I) may be understood as the field equation for the 
induced connection {2to, 21} = {SU*} and appears in three-vector notation as 

-A2lo + 4 (V ■ W\ = V • (i + 2W\ = 47ra s ( (2) j - (3) Jo + 2/ ) (III.59a) 
-A21 - 4 (V x V) = V x (ft - 2V \ = Ana s (j 2 - j 3 + 2l \ . (III.59b) 

This three-vector formalism for the induced "exchange bundle" is also well-suited in 
order to consider (albeit somewhat superficially) the question of boundary conditions. 
Surely, one will require that the electrostatic potentials ^A fjlll. 17|) will decay at infin- 
ity (r — > oo) like an ordinary Coulomb potential (a = 1, 2, 3) 

lim {a U (f) = e a ■ — (111.60) 

i — >oo r 
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where e a is a "pseudo-charge number". Of course for the non-identical particle (a = 1) 
one will put e% = 1 which is the actual charge number, because the non-identical particle 
is not subjected to the exchange forces; but for the two entangled particles one naturally 
expects e a 7^ — 1 (a = 2,3). Recall here the fact that the total potential A s M (IIII.2a,[) of 
both entangled particles feels exclusively the total current j*, cf. ( 1III.3I) 

V^FV = DA\ = -AA\ = 4Tra sJ s u (111.61) 

so that we can require at spatial infinity (r — ► oo) 

^A = ( 2 k + ( 3 k = -2 • ^ . (111.62) 

r 

This means that one always has 

e s = e 2 + e 3 = -2 (111.63) 

but not necessarily e 2 = e?, = — 1 since the exchange interaction can transfer entanglement 
charge S e from one identical particle to the other. 

In order to discuss this charge exchange a little bit more thoroughly, integrate over the 
equation ( 1III.59aj) and apply also Gaufi' integral theorem in order to find under observation 



of the normalization conditions ( 1III.49I) 



j> (i + 21^) • alS = 8na s J al 3 fl (r) , (111.64) 

(Soo) 

where is a closed two-surface at infinity (r — > oo) with surface element dS. Here the 
two-dimensional surface integral of the vector field W ( 1III.52al) can now be assumed to be 
zero because it is a nearby hypothesis that Bo should decay at spatial infinity (r — > oo) 
at least like r _1 (monopole) and B at least like r~ 2 (dipole), which yields W ~ r -3 . But 
for such a situation the equation (1III.63P then links the entanglement charge S e = ( e2 ~ e3 ^ 
to the integrated entanglement charge density lo(f) through 

€ ■ dS = 8na s 5 C = 8vra s / d 3 rl (r) , (111.65) 

(Soo) 
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i.e. the entanglement charge 

f d 3 rl (r) (111.66) 



is shifted from one identical particle to the other. 

Applying quite similar arguments as for the present equation (IIII.59a[) also to the 



former equations ( 1III.29bl) -( lIII.29cT) for the individual particles one concludes that the 
pseudo-charges e a ( 1111.60)) for the identical particles are given in terms of the entanglement 



charge S c ( \lll.66\i through 



e 2 = -1 + S e (III.67a) 
e 3 = -l - 5 e . (III.67b) 

Thus summarizing the results, one arrives at the conclusion that it is only the total charge 
number e s (IIH.63P (of a subsystem of identical particles) which is generally fixed and 
leads to a well-defined electrostatic potential ( s 'Aq ( 1111.62)) at spatial infinity, whereas the 
individual particles do carry charge numbers (here e 2 and 63) which are neither integer nor 
are generally fixed but depend upon the considered situation (i.e. the concrete quantum 
state). Clearly, the tacitly made presumption is here that the charge number e a of any 
particle (a) is defined by that asymptotic form ( a l4o(r) of its emitted potential (1111.60)) . 
Here, it should become also obvious that future work has to be spent upon a more thorough 
inspection of the boundary conditions of all the concerned gauge fields. 
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IV Gauge- Invariant Energy Functional 

It seems unlikely that the wave functions (or even densities) in an atom or molecule can 
be observed directly in order to test the physical truth of the RST solutions; and therefore 
one has to look for some other physical observable which is essentially determined by the 
RST solutions but on the other hand provides also a handle for the observations. Surely, 
such a quantity will be the total field energy being concentrated in any RST solution 
since this quantity determines uniquely the optical spectrum of the bound system under 
consideration. A nearby definition of the wanted RST energy Et refers of course to 
the total energy momentum density ^T^, cf. equation (III. 2 1 p . This proposal has been 
studied extensively in the preceding paper [T3], however only for the linearized theory 
which has been shown (by the preceding arguments) to miss the exact gauge invariance. 
Therefore we will readily set up once more the wanted energy functional Em , but this time 
with properly taking into account the non-Abelian (and therefore non-linear) character 
of RST. 

The non-relativistic limit Em of the relativistic Em is also worked out; and here 
it is important to preserve the gauge covariance just in this limit. It should also be a 
matter of course that the subsequent exact result for Em must be closely related to the 
previous approximative result of the linearized theory. Indeed, this relationship shows up 
in the formal identity of the present exact result ( 1IV.20I) with the analogous expression 
of the linearized theory, cf. equation (1IV. 14|) of ref. [TSj. Thus the present progress refers 



here mainly to the individual contributions to the energy functional which are now made 
rigorously gauge-invariant! This invariance becomes manifest through the rigorous use 
of the gauge- covariant derivative D and V, see equations 26[) and flIII.34p . Clearly, 



the expectation is now that such an energy functional, which meets with the symmetry 
requirements of all the modern gauge field theories, should also have something to do with 
physical reality. 
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1. Generalized Poisson Identities 



The Poisson and exchange identities are important constraints for the principle of 
minimal energy and therefore must be considered before the corresponding energy func- 
tional can be set up. The electric Poisson identity is obtained by multiplying through 
any of the three electric source equations (lIII.29al) -( lIII.29c l ) with the corresponding elec- 
trostatic potentials PMo(f), ( - II ko(r), ( m ko(r*) and integrating over whole three-space with 
finally adding up all three results. This procedure yields the electric Poisson identity in 
the following form 

ivf ee ng) (IV.l) 

where the left-hand side is defined in terms of the electric potentials and charge densi- 
ties ( a 'jo(r) through 



4*) = _^L_ J [(v(%) • (v< 2 k ) + (v(%) • (v( 3 k ) + (v< 3 k ) • (v^k 

- ^ J d 3 f [(% ■ « j + ("k ■ (*> jo + (ni ko . (3) Jo] -JlL-J d 3 r 2l (V • W) . 

(IV.2) 

And furthermore, the right-hand side of the identity ( 1IV.1I) is given in terms of the ex- 
change objects through 

»g ^-y / ^sto-'o- (iv.3) 

It should be obvious that here must exist also a magnetic analogy to that electric 
identity (lIV.ip . Indeed this magnetic Poisson identity may be deduced from the magnetic 
part (IIII.30ap - (lIII.30cP of the (non-Abelian) Maxwell equations in just the same manner 
and appears then as follows: 

Nt ] = ng° , (IV.4) 
where the left-hand side of this identity is defined quite analogously to its electric coun- 
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terpart flTOD through 

r (m) fa 



N, 



G 



Ana* 



V x AA • V x A 2 + V x A 2 ) • V x A, 



+ V x A 3 ) • VxAj 



he 

— I d 3 f[ 3 \ • Ai + j 2 • A u + j 3 • A m 



(IV.5) 



and the right-hand side appears again as the magnetic counterpart of Uq (1IV.3[) . i.e. 



n. 



(m) 
G 



he 



d 3 m- V x V 



he 



d 3 f 21-/ . 



(IV.6) 



47ra s y V / 2 

Now with a little bit intuition concerning the relativistic formalism, one easily realizes 
that the electromagnetic Poisson identity 



Ng-N^-ln^-ng*) =0 



(IV.7) 



is built up by certain Lorentz invariants; namely this identity reads explicitly 
J d 3 f (t l \ v ■ ( 2 k^ + ■ ( 3 k^ + ■ ^A" u + 8na s VSTB^) 
+4vra s j d 3 f (% • W + • ( 2 ) / + < m k M • ( 3 ) / - %F) = . 
Here it is made use of the source equation for the object ( 1III.44I) 



(IV.8) 



i'H 



(IV.9) 



whose space-time splitting is given by equations ( 1III.57I) -( 1III.58I) ; and furthermore the 
skew-symmetric tensors ^'A^ (= — ' a U Vfl ) are defined through (a = 1,2,3) 



(a) 4 — V7 4 a _ \7 4 a 



(IV. 10) 



But let us stress here the somewhat strange fact that it is only the left-hand side n'q 
(1IV.2P and (1IV.5[) of the electric and magnetic Poisson identities (lIV.ip and (1IV.4I) 
which play the role of a constraint for the subsequent RST principle of minimal energy. 
Clearly one expects that a similar constraint must be active also for the exchange subsys- 
tem. 
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2. Exchange Identities 



The procedure for deducing the exchange identities is quite similar to that for the pre- 
ceding electric and magnetic Poisson identities. More concretely, one multiplies through 

* 

the electric exchange equation (1III.31a|) with B (and its complex conjugate with B ) and 
then one arrives after some simple mathematical manipulations at the following electric 
exchange identity 



00 



(iv.il) 



Here the left-hand side is given by 



N, 



(h) ^ 



he 
Ana. 



d 3 r 



V'B< 



(v'B 



he it * * 

— I d 3 r ( B h + B h 



o 
he 



d 3 f 2lo [V'W 



(IV.12) 



and the right-hand side by 

(h) ^ he 



71, 



G 



d 3 r % [V-W] - 



he 



d 3 fW 2 



(IV.13) 



Quite analogously, the magnetic exchange identity is obtained from the magnetic ex- 
change equation (1III.31bl) by multiplying through with B (or B, resp.); and this then 
yields the following identity 



(g) 



(IV. 14) 



with the left-hand side Nq^ being defined by 



N, 



(g) 



G 



+ 



he 
he 

he 



d^f 



V x B ] • [V'xS 

he 



V'B) ■ ( V'B) -2& ( B'B 



d 3 fQln[V'W 



he 



Ana* 



d 3 fV 2 - 



Ana, 



d 3 fW 2 - 



d 3 f 21- ( V x V 

he 



d 3 r\B'h + B-h 



(IV. 15) 
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and similarly the right-hand side Uq* by 

,0s) ^ ^ c /^3-> T ?2 ^ 



= d i fV 2 + d 6 r%a(V'W) . (IV.16) 

4iia s J 8na s ' v 



The exchange identity fllV.lip of the electric type can again be combined with the ex- 



change identity of the magnetic type (1IV.14I) in order to yield the total exchange identity 



n g ~ n G ~ ( N g ~ n^) = . (IV.17) 
This may be written again in terms of Lorentz invariants as 

d 3 r [G^G^ - 47ra s WB^ ) + 4vra s / d 3 f + B^ + 8vra s B^B„ V ) = 



(IV.18) 

and is of course the exchange counterpart of the electromagnetic case (IIV.7[) - (IIV.8[) . Here 



we have made use also of the fact that the exchange field strength G^ v (IH.lldj) can be 
expressed in terms of the gauge-covariant derivative (1III.34I) as 



G>„ = 3^B U - B V B^ . (IV. 19) 

However observe again that for the subsequent energy functional Em one cannot use the 
Lorentz invariants (1IV.18P immediately but rather one has to resort to their non-invariant 
parts Nq and NqP separately; the complementary parts and Uq do work as simple 
constants which must be held fixed for the variational process! 

3. Relativistic Energy Functional 

The essential structure of the RST energy functional Em has already been discussed 
in the preceding paper [13], albeit only for the linearized theory. But this principal form 
of E[t] does also apply to the original non-Abelian (and therefore non-linear) theory, i.e. 
we put 

Em = 4]° + £ A D (a) • N D (a) + A«. (ivf - ng>) 

^ V ' (IV.20) 

+ AS - (jv£° ~ + )&• - ng») + A«. (jv£> - nf 
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The leading feature of this functional consists in its composition of a physical part 
(i.e. Errq) and of the constraints of electromagnetic and exchange type, besides the nor- 
malization conditions (1111.491) for the wave functions. The four gauge field constraints are 
based upon the demand that the objects i\^ e) flTOj) . fHV3D . Jv£ h) flTvU2l and JV^ 



(IIV.15I) must be kept fixed for the variational process, i.e. one interpretes the function- 



als rig (1IV-3I) . (1IV.6I) . Uq (1IV.13I) and Uq (1IV.16I) as constants not undergoing 
the variational manipulations 6n^ = Suq^ = Sn^ = 6riQ = 0). Especially, for 
a system of non-identical particles (such as the previously considered positronium [T3] ) 
all exchange fields do vanish and thus the real numbers 

^G ? ^G ' ^G ' ^G diJCG all zero 

per se. For such a situation the variational procedure does reproduce correctly the RST 
system of mass eigenvalue and gauge field equations, see ref. [13J. But the progress with 
the present generalized case (1IV.20I) is now that both the RST mass eigenvalue system 
(IIH.22al) - flIH.22cl) plus the gauge field system (1III.29al) - flIII.31bl) does reappear in a gauge- 
invariant form from the extremalization procedure with fixed numbers Uq , Uq , Uq , tiq 
(not to be interpreted as functionals of the gauge fields): 

5E m = (IV.21) 

(5ng = 5r4 m) = <5r4 h) = 5n^ = 0) . 

Thereby the Lagrangean multipliers turn out as 

A D (l) = M lC 2 , A D (2) = -M 2 c 2 , A d (3) = -M 3 c 2 (IV.22a) 

Aj. m) =-Ag } =2 (IV. 22b) 

\M=-\®=2. (IV.22c) 

Recall here also the former discussion below equation (1III.47P concerning the right choice 



of the normalization condition for the wave functions ip a (r). The reason for that preference 
of ( 1III.49I) over (1III.48I) becomes now clear because it is just this specific condition (1III.49I) 
which (by virtue of equations (1IV.22al) ) lets emerge the mass eigenvalues M a as desired 
in the mass eigenvalue equations flIII.22ap - flIII.22cp . 
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If the system of extremal equations due to ( 1IV.21I) could be solved exactly, one could 
substitute the exact solutions back into the energy functional E^ ( 1IV.20I) and would then 
find that the total RST energy E^, concentrated in that hypothetical exact solution, is 
determined by the physical part E^O alone. This physical part itself is the sum of all the 
partial energies of the coupled RST field system, i.e. 

= Z 2 (l) -M p c 2 + Z\ 2) -M e c 2 + Zf 3) -M e c 2 + 2T kin + e£ ] - ^ m) - + E { ^ • (IV.23) 

Here, the mass renormalization factors Z 2 ^ (a = 1, 2, 3) are defined in terms of the Dirac 
spinors ^) a {r) as follows 

Z 2 {a) = J d 3 r^ a (r)^ a (r) . (IV.24) 

Furthermore, the total kinetic energy T^ n is the sum of the three single-particle ener- 
gies T kin (a) 

T kin = T kin (l) + T kin (2) + T kin (3) , (IV.25) 
with the kinetic energy of the first particle (a = 1) being given by [13] 



he 

T kin (l) = i— I c? 3 r^i(r)7- W>i(r) (IV.26) 



2 

and similarly for both identical particles (a = 2, 3) 

T kin (a) = -i— I d 3 f^ a (r)7- V^ a (f) . (IV.27) 

The last contribution to the physical energy comes from the gauge field 

and subdivides into the electrostatic energy E^ 

E * = 4^T J rf3f v^ 1 * ^2 + ^2 * ^3 + ^3 * Ei) , (IV.28) 
magnetostatic energy E^ 1 

4 m) = / d3 ? [Hi 'H 2 + H 2 'H 3 + H 3 - H^j , (IV.29) 
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exchange energy of the electric type 



j#> = ^/ ,/V(a'.a-) . av.m 



and finally the exchange energy Eg* of the magnetic type 



E <g = Jc_ I d 3 f (Y'Y) . (IV.3i; 



An interesting point of this RST principle of minimal energy (IIV.2ip refers now to 



the fact that only the Lagrangean multipliers A D (a) (jlV.22aj) . being linked to the wave 
function normalization ( 1III.49I) . depend upon the special RST solution; but the four mul- 
tipliers Xq\ Xq] Aq , X^q referring to the Poisson and exchange constraints (1IV.1I) . (1IV.4I) . 
( II V. lip and (IIV.14I) are universal constants, cf. (IIV.22bj) -( HV.22cl) . Therefore it is possible 



to collect all four gauge field constraints into a single one, i.e. we may put for the RST 
functional E [T] ljIV20jl 

3 

E m = + Ad(o) ■ N D (a) + X G ■ (N G - n G ) (IV.32) 

0=1 

where for the single residual multiplier one has Aq = —2. The compact gauge field 
constraint is then given by the identity 

N G = n G , (IV.33) 

with Ag working as the proper constraint condition whereas n G has the status of a fixed 
number (5n G = 0) 

Ag = JV£ } - A^ - A^ h) + A^ g) (IV.34a) 
n G = n g } - 4 m) - ng 5 + ng } . (IV.34b) 

Finally it is also worth while to point out to the circumstance that the mass eigenval- 
ues M a do not only emerge in the mass eigenvalue equations ( lIII.22al) - (IIII.22cj) but also 
in the gauge field equations (IIII.29bl) - (IIII.31bl) . namely implicitly in form of the exchange 
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length om ( 1HI.13I) occurring in the transformed potential 21m ( 1HI.12I) and also in the 
exchange field strength X ( 1III.25af ). This entails certain complications when one wishes 
to exploit the principle of minimal energy (IIV.2 1[) for developing variational techniques 



in order to obtain approximate solutions for the RST eigenvalue problem (consisting of 
the matter equations (1III.22al) -( lIII.22cl) plus the gauge field equations (1III.29al) - (1III.31bl) ). 
It is true, the exchange length au is to be kept fixed during the variational procedure 
(IIV.2ip . but for solving the extremal equations due to the functional Em one has to ex- 



press <2m in terms of the mass junctionals, i.e. one has to render more precise the former 
definition of au ( 1IH.13I) to 

1 (M\2]-M\ 3 ])c 

— = V E w> (IV.35) 

(for a discussion of the mass functionals M[ a ] see ref. [13J). Clearly by this arrangement 
(IIV.35P the RST eigenvalue problem is revealed to be essentially an integro-differential 



system. However this complication is not present for systems of different particles where 
the exchange subsystem becomes trivial (i.e. = 0), see below for the treatment of 
positronium. See also Appendix A for the manifestly gauge-invariant form of the energy 
functional E^y 

4- Non-Relativistic Variational Principle 

Surely, it is very satisfying to have a rigorously relativistic energy functional because 
this admits to predict the frequencies of spectral lines as precisely as the observational 
data do require it in order to test the quality of the theory. However for the practical 
applications it will mostly not be possible to solve exactly the relativistic eigenvalue prob- 
lem so that one is forced to resort to approximative solutions. But when the solutions 
are not exact, it makes no sense to insert them into the exact (relativistic) energy func- 
tional; but rather it is sufficient (and simpler) to use a similar approximation of the energy 
functional, i.e. we have to look for its non-relativistic limit. In this sense we must now 
deduce (from the relativistic predecessors) the non-relativistic approximations of both the 
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corresponding RST field equations and the energy functional. Logically, both concepts 
must fit together in the sense that it does not matter whether one first establishes the 
non-relativistic functional QE[t])> sa y) an d then looks for the non-relativistic field equa- 
tions as the extremal equations of Em ; or whether one deduces the non-relativistic field 
equations directly from their relativistic predecessors ( lIII.22al) - (1III.22cl) . i.e. the following 
approximation scheme must be commutative: 



Relativistic 


non-relativistic 


Non- Relativistic 


Field Equations 


approximation 


Field Equations 



p 
•x: 

o 
o 

p 
ft 



p 
•x: 

o 
o 

p 
ft 



Relativistic Principle 
of Minimal Energy 


non-relativistic 


Non-Relativistic Principle 
of Minimal Energy 


approximation 



The point here is that we have to establish the non-relativistic approximations in such 
a way that the magnetic gauge covariance (or invariance, resp.) is properly respected; 
this is not the case in the preceding paper [T3] where the process of linearization did spoil 
the strict magnetic gauge invariance/covariance. 

Non-Relativistic Eigenvalue Equations 

First let us deduce the non-relativistic energy eigenvalue equations directly from 
their relativistic predecessors (1III.22a[) - (1III.22c[) . The method of non-relativistic approx- 
imation consists here in resolving approximately the equations for the "positive" Pauli 
spinors ^ip + {r) with respect to the "negative" Pauli spinors ^ip^(r) and substituting 
this back into the equation for the "negative" Pauli spinors. Thus in the first step, the 
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"negative" Pauli spinors do appear in the following form [13J: 



(1 V-(r) - R(i) (IV.36a) 
( 2 V-(r) - R (n) (2 V+(r1 + r ( ii) (3 V+(r) (IV.36b) 
( 3 V-(r) - ^(iii) (3 V+(0 + r(in) (2 V+W • (IV.36c) 

This demonstrates again the exchange coupling of both identical particles (a = 2, 3) 
as well as the corresponding decoupling of the first (non-identical) particle (a = 1). The 
resolvent operators, mediating the approximative transition to the non-relativistic regime, 
are given by 

«<">=-2^ J -(™') ' f M=2h S - 8 (IV - 37b) 
A "^ = -2K-c 3 -^ + lA ' a ) • ^^ikl 3 ' 6 - (IV ' 37C) 



But now the crucial point is here that those non-relativistic approximations ( 1IV.36al) 



( 1IV.36cl) must be covariant under the group of magnetic gauge transformations! These 
transformations are defined for the magnetic potentials A\, An, Am through the equations 
( 1III.20aft - 11III.20cl) . furthermore for the exchange potential {B^} = {B ,—B} through 
equation ( 1III.9I) . and finally for the Pauli spinors through 



(<V ± (f) = e~ taa ■ {a V±(r) • (IV.38) 

Clearly, the latter transformation law is entailed on the two-component Pauli spinors 
^ a V±(r) by the corresponding law (IIII.6P for the four-component Dirac spinors ^(r), 



namely just via the direct sum construction ( 1III.21I) . But the point with this group 
of magnetic gauge transformations is now that those approximative relations ( 1IV.36al) - 



(1IV.36cp must transform covariantly under the magnetic gauge transformations. Indeed, 
it is easy to see that this magnetic covariance requirement transcribes to the resolvent 
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e n 


• e~ 


-iai 


= Rm 




• e" 


-ia.2 


= R(ii) 


ioi3 . p' 


• e" 


-iaz 


— -R(in) 



operators in the following form: 

(IV.39a) 

, e ia2 ■ f In) • e~ ia3 = f ( n) (IV.39b) 
, e^-f{ m) -e- ia2 =f m . (IV.39c) 

But this required transformation behavior of the resolvent operators can easily be ver- 
ified by simply evoking the corresponding transformation laws for the three-vector po- 
tentials B, Ai, An, Am as given by equations ( 1III.9I) and ( lIII.20al) -( 1III.20cl) . Thus, the 
magnetic gauge covariance is preserved in the first step. 

The second step consists now in substituting those appropriate "negative" spinors 
( a V-(r) flIV.36al) -f lI\r36c1) back into the eigenvalue equations flHI.22al) -f l!TL22c1) for just 
the negative spinors. This procedure then leads us to the following Pauli-like equations 
for the "positive" spinors ^(f + (r), a= 1,2,3 

H m W<p+(r) = E m ■ W<p+(r) (IV.40a) 
H (n) ( 2 V+(r) + $(n) (3 V+W = #p(2) • (2 V+W (IV.40b) 
U m (3 V+(r1 + 6(iii) (2 V+(r) = £p(3) ■ (3 V+W • (IV.40c) 

Here, the Pauli energy eigenvalues -Ep( a ) are non-relativistic approximations of the former 
mass eigenvalues M a , i.e. 

£ P(1) - - (M p + MO c 2 (IV.41a) 
E P{2) ^ (M 2 - M e ) c 2 (IV.41b) 
E P{3) -(M 3 -M e )c 2 , (IV.41c) 
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furthermore the Pauli Hamiltonians (Tt) are given by 



T^(iii) 



h 2 



2M n 



h 2 



2M e 

h 2 

2M P 



V + iAi) -a> Hi 

V + iAu) -B'Hn 

2 



V + iA 



in 



o-H 



in 



+ he ( % 
- he ( n k 
he ( n % 



(IV.42a) 
(IV.42b) 
(IV.42c) 



and finally the exchange Hamiltonians (t)) are found to be of the following form: 

h 2 



-hcBn 



2M P 



6(in) = -hcB + 



/i 2 



2M P 



a • Y - i ( V • B ) - 2iB • V + zA m 
-t • Y - / ( V' • 5^1 - 2iB • ( V + zk*ii 



(IV.43a) 
(IV.43b) 



Observe here that the present result (lIV.40ap - (lIV.40cP for the non-relativistic eigen- 
value equations is the strictly gauge-covariant version of the former results (III. 74) and 
(III.76a)-(III.76b) of ref. [13]. Indeed, the present Pauli Hamiltonians 7t^(i,ii,iii) (IIV.42a|) - 
(IIV.42cl) are exactly the same as in the preceding paper [T3j, see equations (III. 75) and 
(III. 77) of that paper. The reason for this is that the former Pauli Hamiltonians of ref. p3] 
are already invariant with respect to the magnetic gauge transformations: 

(IV.44a) 
(IV.44b) 
(IV.44c) 

Therefore the progress made here refers only to the exchange Hamiltonians f)(n) and t)(iii) 
flIV.43al) -f lI\r43bl) . These objects appear now manifestly gauge-invariant because they 
contain the gauge-covariant gradient operator V 26[) rather than its counterpart V 
of ref. [13]. Thus the magnetic covariance of the exchange Hamiltonians emerges in the 
following form: 



iai nj 1 

e n 


e" 


-ia\ 


= H(i) 


e n (jj) 


e" 


-io<2 


= 'Hiii) 


ia 3 nj 1 

e n (I ) 


e" 


-ia 3 


= T^(iii) 



,J«2 U I 



'""f)(n)e 



(III) 



= f)(II) 
^ = 6(111) • 



(IV.45a) 
(IV.45b) 
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Summarizing, it has been demonstrated by explicit construction that it is possible to 
deduce the non-relativistic eigenvalue equations from their relativistic originals under the 
desired preservation of the magnetic gauge covariance. Now it should be a matter of course 
that one wishes to see also the corresponding gauge-invariant energy functional (E[t], say) 
whose extremal equations would just coincide with the present non-relativistic eigenvalue 
equations (1IV.40al) -( lIV.40cl) . Clearly, a nearby guess suggests here to deduce this invariant 
functional Em from its relativistic original Em (IIV.20I) . again under preservation of the 
magnetic gauge invariance. Moreover this non-relativistic functional E[t] is required to 
produce also, through its extremalization with respect to the gauge fields, just the non- 
relativistic version of the original gauge field equations (1III.29a[) - (1III.31b[) . 

N on- Relativistic Gauge Field Equations 

However, resorting to such a non-relativistic limit of the gauge field equations would 
imply that one neglects (at least in lowest-order approximation) the fields of magnetic 
type altogether! The reason is that the magnetic fields couple to the particle velocities 
which in the non-relativistic domain are negligibly small. On the other hand, the present 
gauge invariance refers exclusively to the magnetic gauge transformations; and therefore 
the neglection of all the magnetic-type fields would eliminate also the magnetic gauge 
structure. Therefore in order to preserve this gauge structure for the non-relativistic limit 
of the gauge field subsystem, one should neglect only certain products of the magnetic- 
type fields provided this neglection does not spoil the magnetic gauge structure. However 
such products of magnetic fields, being well-suited for omission, do occur at almost all 
places of the gauge field subsystem, for instance in connection with the equations flIII.23bj) - 
( 1III.23cl) . ( 1III.24bl) -( lIII.24c!) or ( 1III.31bl) . Therefore, in order to not pervert the essentially 



magnetic character of the gauge field configurations, we keep here all the gauge field 
equations formally unchanged and do insert therein only the non-relativistic limits for 
the densities ^jo, j a , h , h (see ref. p3]). Clearly, such a procedure implies that all the 
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Poisson identities and exchange identities do formally apply also for the present non- 
relativistic limit with merely the constraints N, n being replaced by their non-relativistic 
versions N,n. 

Thus the last constraint which must be be considered here refers to the wave function 
normalization ( 1III.49I) . The non-relativistic version hereof consists simply in using the 
non-relativistic approximation of the Dirac densities ^ko^r). The relativistic density 
reads [T3] 

Mjfe (f) = (a V+(r) (a V+(r) + (< V_(r) (a V~(r) , (IV.46) 

and since the non-relativistic limit consists in disregarding the "negative" Pauli spinors 
W(p_(r) one concludes for the non-relativistic approximation ( Nd, say) 

N D N D = / d 3 r (a V+(r) ( °V+(r) -1 = 0. (IV.47) 



N on- Relativistic Energy Functional 

But now that the non-relativistic (however gauge-invariant) forms of all the constraints 
have been revealed, one can easily write down the non-relativistic approximation E[t] of 
the original relativistic energy functional Em ( 1IV.20I) : 

Em = + £ A p(a) • N D(a) + Ag> ■ (N« - *#) + a£> ■ (N<f> - ftW) 

a=1 3 (IV.48) 

+ Ag } ■ (4 g) - 4 g) ) = + £ A P (a) ■ N D( a) + A G ■ (N G - n G ) . 

0=1 

Here the three non-relativistic gauge field constraints with Lagrangean multipliers Xq\ Aq \ A 
may be collected again into the single constraint Nq with Lagrangean multiplier Ag = —2 

N G = N§ } - Nff + Njf 5 (IV.49a) 
n G = fig) - ng + ng g) , (IV.49b) 
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cf. the corresponding relativistic arrang ement fliy.34aMIV.34hj) . Clearly, the principal 
structure of this nonrelativistic functional ( ]IV.48j) is the same as for the relativistic origi- 
nal E[t], namely the sum of a physical part (i.e. Ep^) and of the constraints. The physical 
part is of course the collection of all those physical contributions which are already present 
also in the original relativistic version fli\T23l) . however now in their non-relativistic 
approximation, i.e. 

E[™ } = E kin + 4 C) + 4 m) - 4 h) + E® + hcj d*r (H x • ^S x + H u • ™§n + H m • ^„i) ■ 

(IV.50) 

It is true, this is the same form as in the preceding paper, see equation ( 1IV.21I) of 
ref. [13]; however, the invariance of any individual term with respect to the magnetic 
gauge transformations is now guaranteed! Namely, the point here is that the four gauge 
field contributions E^ , E^ , Eq , E^ are formally the same as their relativistic originals 
( liV.28l) -( ]IV.31l) if expressed in terms of the corresponding field strengths E a , H a , X ,Y; but 
observe here the change of the magnetic contribution in the non-relativistic functional Em 
(IIV.48I) which no longer contains the magnetic constraint — tiq of the relativistic 
case ( 1IV.20I) . Indeed this magnetic constraint is eliminated in favour of the positive sign 
of E^ and of the emergence of the dipole interaction energy (~ H • S). Moreover, the 
gauge invariance of the magnetic dipole energy is here manifest since the non-relativistic 
Dirac dipole densities ^ m 'S a are defined by [13] 



(m) 5i = ^- (1 V^ ( V (IV.Sla) 

{m) S2 = ^ c i2 W% + (iv.sib) 

H ^^2i^ (3 H^ ( V- (IV.51c) 

And finally, the total kinetic energy E^ n in flIV.50p of course inherits its sum structure 
from its relativistic predecessor ( 1IV.25j) . i.e. 

3 

^kin = ^^kin(a) , (IV. 52) 
a=l 
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with the manifestly gauge-invariant one-particle contributions 

h 2 r 



^ki„(l) = -r^r- / rf 3 f (V+ (V + ill)' ( V (IV.53a) 
^ ki „(2) = | d 3 f <V+ (V + il n ) 2 (2 V+ (IV.53b) 

^kin(3) = -^fj ( V + (V + iA„i) 2 {3 V+ • (IV.53c) 



But now that the non-relativistic functional Em is precisely fixed in the form (1IV.48P it 
is a standard matter to deduce the non-relativistic RST field equations through extremal- 
ization of that functional: The corresponding variational procedure with respect to the 
Pauli spinors ( a ^(p±(r) (or ( a ^(p±(r), resp.) yields the system of gauge- covariant eigenvalue 
equations flIV.40ap - pV.40cp with the Lagrangean multipliers 



A p(o ) = -Ep(a) (IV. 54a) 

X§ = -\^= -A[, h) = Ag } = -2 , (IV.54b) 

see Appendix B. And similarly, the extremalization of Em with respect to the gauge 
fields ^Ao(r), A a , Bq, B (or their complex conjugates, resp.) lets emerge the electro- 
magnetic Poisson equations and their (gauge-covariant) exchange counterparts (see Ap- 
pendix C). Thus the commutativity of the logical arrangement on p. EH is safely validated 
which in turn supports the general logical consistency of RST. (The variational deduction 
of the non-relativistic eigenvalue equations ( irV.40aP -( irV.40cl ) from the non-relativistic 
principle of minimal energy, SKm = 0, is a very instructive exercise in dealing with 
the gauge-covariant structures, see Appendix B. See also Appendix C for the variational 
deduction of the gauge field equations.) 
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V Ortho- and Para-Positronium 



As a practical demonstration for the usefulness of the constructed energy functional 
Em, or its non-relativistic approximation Em, resp., one may consider now a two-particle 
system (N = 2) which is simple enough in order to avoid unnecessary complications but, 
on the other hand, is sufficiently non-trivial in order to display the physical correctness of 
the theoretical construction. Thus for the sake of demonstration, we turn to positronium 
whose constituents (i. e. electron and positron) have identical rest masses (M e = M p = M) 
but different sign of charge so that both constituents must count as non-identical particles. 
Consequently, these particles cannot feel the exchange interactions, which entails the 
vanishing of all the exchange objects (B^ = 0, G^ v = 0, etc.). But this does not mean that 
the demonstration becomes trivial because positronium is known to possess a sufficiently 
rich level structure [llj. This is mainly due to the well-known ortho/para dichotomy 
which occurs in a similar way also for muonium and hydrogen, albeit with different rest 
masses [10] . Thus in the present context, positronium provides us with a test case for our 
energy functional E^} (or E[t], resp.) in a two-fold respect: (i) one can reveal the specific 
way in which the ortho/para dichotomy emerges in RST (in contrast to the conventional 
theory), and (ii) one can (at least qualitatively) discuss the energy difference between 
certain ortho- and para-levels in order to compare the corresponding RST predictions 
with those of the conventional theory and with the observations. 

Naturally, such a double program is to be carried through in two steps: In this section 
(V.), we first study the emergence of the ortho/para dichotomy within the RST formal- 
ism which is based upon the Whitney sum of one-particle bundles, in contrast to the 
conventional formalism, which deals with tensor products of one-particle Hilbert spaces. 
And in the next section (VI.) we apply the constructed RST energy functional Em (or 
Em, resp.) in order to discuss the origin of the energy difference between the ortho- and 
para-levels. This energy difference appears to be caused not only through different mag- 
netic interactions but also by the different anisotropy of the electric interaction potentials, 
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which itself is due to the different pattern of angular momentum inherent in the ortho- 
and para-configurations. 



1. Conventional Multiplet Structure 



It is true, our model demonstration below (i. e. positronium) does not consist of 
identical fermions; but for the sake of comparison of the conventional and RST situations 
it is very instructive to first consider a system of identical fermions (e. g. the two helium 
electrons) from the conventional viewpoint. 

According to the conventional spin- statistics theorem, i. e. ll Pauli principle" in pop- 
ular terms [TU], any stationary bound system of identical fermions (e.g. spin-^ particles) 
occupies a skew- symmetric quantum state ^^f(fi, si; r 2 , s 2 ): 



Analogously, the bosons (i. e. integer-spin particles) can occur only in symmetric states 



For a classification of such states one resorts to the assumption that their number is 
preserved when one thinks of the spin-spin interactions as being switched-off, with simul- 
taneous neglection of the relativistic effects. Namely, the point with this hypothesis is that 
the wave functions ( IV. II) (or analogously also ( 1V.2I) ) can then be factorized approximately 
into an orbital part ^(r*!, r 2 ) and a spin part x(^i, s 2 ): 



But since the Pauli principle ( IV. II) demands that the total wave function must be an- 
tisymmetric under particle permutation, the factorization (IV. 3D admits exclusively the 
following combinations of symmetry properties for the orbital and spin parts 








(V.3) 



X(*2,si) 



(V.4a) 



ip(r 2 , ft) x(si, s 2 ) = +x(s 2 , st) . 



(V.4b) 
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Here a closer inspection of the configuration space of a two-spin system shows that 
the antisymmetric spin configuration ( lV.4al) is actually a singlet state X[o,o] 

1 

7! 



X[o,o] 



Co""® CI' 5 -c 



1 _1 

2 ' 2 



1 1 

2 ■ 2 



Co 2 ' 



(V.5) 



This antisymmetric decomposition of the spin part xi^it^z) refers to the standard spin 



— i — 

basis {Co 2 2 } of the two-dimensional unitary space, i.e. 



sis 



1) 



h 



S>0 



2 2 



^ 1-1-1 



1+1 
± 2 Co 

0, 



(V.6a) 

(V.6b) 
(V.6c) 



where {s 2 , s z , 1 2 } are the corresponding one-particle operators for spin (s = and orbital 
angular momentum (I = 0, 1, 2, 3, ...), resp. Thus the antisymmetric singlet state X[o,o] has 
zero values for both the total spin S (— s*i © s 2 ) and its ^-component S z (= s^m © s z (2)), 
i. e. 

J ?2 X [o,o] = ^5(5+l)x[o,o] = (V.7a) 
4x[o,o] = . (V.7b) 

In contrast to this singlet structure, the symmetric spin parts (1V.4bjl do form a triplet 
system 

Xm = [C<H ® Co 1 '"" + Co 5 '" 5 ® Co 5 ' 5 ] (V.8a) 

X[i,i] = Co 1 ' 5 ® Co 5 ' 5 (V.8b) 
=d'" 5 ®Co 5 '" 5 (V.8c) 

such that 



S 2 X[i,o] = h 2 S(S + 1) X[i,o) = 2ft 2 X[i,o] 
S z X[i,o] = 



(V.9a) 
(V.9b) 
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and 

?2_ . _ fe2 CVO , 1 \ . _ Qfc2 



S"X[i,±i] = V S(S + 1) X[i,±i] = 2ft 2 X[i,±i] (V.lOa) 
&X[i,±i] = ±ftX[i,±i] • (V.lOb) 



Observe here that the singlet state X\o,o] HV.5|) and the triplet state X[i,o] (TV-8ap are 



degenerate with respect to the total spin component S z , cf. ( lV.7bl) and ( 1V.9bl) . but this 
does not imply that the singlet state ty-p with symmetric orbital and antisymmetric spin 
part 

%? P = ^ + {f u f 2 ) ®X[o,o] 

(V.ll) 

(■0+(fi,f 2 ) = V+(r 2 ,fi)) 

has the same energy eigenvalue as the corresponding triplet state ^o- The reason is here 
that the orbital part ip + (ri,r2) of the para-state typ (IV. lift is symmetric, whereas the 
orbital part ^-(r^f-j) of as a member of the ortho-system is antisymmetric: 

*o = ip-(r h r 2 ) <S> X[i,o] 

(V.12) 

(^-(fi,f 2 ) = -ip-(r 2 ,ri)) . 

This symmetry arrangement implies namely that for the para- state y v (jVTTTll the two 
identical fermions tend to approach one another, whereas for the ortho-state ( IV. 121) 
they tend to retreat from one another. Therefore one expects, since the electrostatic in- 
teraction of both identical fermions is repulsive, that the para-state (i. e. the singlet typ 
(jV.lip ) has higher energy than the ortho-states (such as ()V".12j) ) which themselves 
should constitute a degenerate triplet. Possibly the latter (electrostatic) degeneration will 
be eliminated as soon as one takes into account also the magnetic interactions and the rel- 
ativistic corrections. However, the typical energy scale, being associated with these latter 
types of interactions (~ a s 2 ), is distinctly smaller than that of the typical electrostatic 
interactions (~ some few electron volts [10]); and therefore the energy difference between 
triplet states ( 3 Si) and singlet states ( 1 So) in an atom is mainly of electrostatic origin! 
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For instance, the energy difference between the singlet state Is 2s """So and the triplet state 
ls2s 3 Si in a helium atom amounts to (roughly) 0,8eV [10J. Nevertheless, the helium 
groundstate Is 2 1 So is of the para-type ty-p (IV. lip since it is only for this configuration 
that both electrons can have the lowest possible principal quantum number (n = 1). 

However, this conventional picture of the multiplet structure and the associated en- 
ergy eigenvalues does not apply to positronium, since the latter sytem does of course 
not consist of identical fermions. This implies that the Pauli principle is not valid here; 
and consequently that combination ( 1V.4aP - fV.4bp of permutation symmetries, concern- 
ing the orbital and spin parts of the positronium wave function, does not apply here. 
Nevertheless, there also arises a conventional ortho/para-dichotomy for positronium, but 
with the associated energy difference between the singlet ( 1 S ) and triplet system ( 3 Si) 
being caused now by the magnetic interactions! More concretely, the non-relativistic 
positronium groundstate will be of the form 



*(ri,si;f 2 ,s 2 ) = ^(ri,r 2 ) <g> < 



X[o,o] 

X[i,i) (V13) 

X[i,o] 
, X[i,-i] 

and thus is four-fold degenerate in this lowest order of approximation, since there is 
no spin-induced symmetry restriction for the orbital factor ip(fi, r 2 ). But this four-fold 
degeneracy of the positronium groundstate becomes eliminated in reality, namely mainly 
as a consequence of the spin-spin interactions (~ s\ • s 2 ) of both non-identical fermions, 
see ref. [10J. (For the classification of the positronium states see ref. |llj.) However, 
the energy difference of the ortho-positronium state 1 3 Si and the para-positronium state 
1 1 So amounts to only 8 • 10~ 4 [eV], see ref. [10]; and this is thousand times smaller than 
the ortho/para splitting of the electronic helium states quoted above. 

Summarizing, the conventional theory seems able to predict the ortho/para splittings 
of both the electric and magnetic type in excellent agreement with the observations; and 
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thus the question arises now: in which way can there emerge a comparable multiplet 
structure in RST and what is the precise physical nature of the corresponding ortho/para 
energy splitting? 



2. Ortho/Para Dichotomy in RST 



Since the conventional theory and RST are based upon very different mathematical 
structures (i. e. tensor products vs. Whitney sums) it should not come as a surprise that 
the RST multiplets do arise in a totally other way. Recall here that in a conventional 
entangled iV-particle state ^(fi, r 2 , f^) the individual particles cannot be linked to 
single-particle states but rather do participate in that common iV-particle state $ in a 
completely equivalent way. In contrast to this democratic participation, any individual 
particle of an iV-particle state \P(r) does occupy in RST a well-defined one-particle state 
^ B (r} (a = 1, N) in agreement with the Whitney sum construction 



Anyone of these individual RST particles emits a unique electromagnetic potential 
A a ll (a = 1, N) to be determined from the Maxwell equations, e. g. flH.12aD -f lH.12cD for 
the case of a three-particle system. Additionally, any individual pair of identical particles 
does cooperatively admit also an exchange potential (B^) which is to be determined from 
the exchange part of the Maxwell equations, e. g. ( 1II.12dl) for the considered three-particle 
system. Subsequently we will inspect the positronium system (N = 2) which is built up by 
two non-identical particles (a = 1, 2) so that the exchange potential is zero (B^ = 0) and 
the Maxwell equations flII.12aD - flII.12cD become simplified to the following form (a = 1, 2) 



V(r) = ^i(r) © if) 2 (r) © ... © ip N (r) . 



(V.14) 



F 



<a 



(V.15a) 



y/x pa 



47ra s j 



■ a 



v 



(V.15b) 
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For the stationary bound systems, the latter equation becomes simplified to the ordinary 
Poisson equation for the electric potentials ( a l4 (f) 

A^A (r) = -47ra s Wjo(r) , (V.16) 

and analogously also for the magnetic three- vector potentials A a (f) 

AA a (r) = -Ana s j a (f) . (V.17) 

Clearly this is a special case of the general situation described by the equations (10.3j> - 
f laib]) and flUTjl - dOSbl of the appendix. 

Concerning now the corresponding matter equations, it seems reasonable to adopt 
the hypothesis that both particles do occupy physically equivalent one-particle quantum 
states since their rest masses are identical putting M e = M p = M) and their charges 
differ merely by sign. But obviously these equivalent states can occur in two forms; namely 
either the magnetic fields are parallel 

ortho-positronium: Ai(r) = A 2 (f) = A b (r) (V.18a) 

Hi(r) = H 2 {r) = H b {r) (V.18b) 

fci(r) = -k 2 (f) = k b (f) (V.18c) 

^(f) = j 2 (r) = j b (r) = k b (r) ; (V.18d) 

or they are anti-parallel 

para-positronium: Ai(f) = — A 2 (f) = A p (r) (V.19a) 

H x (f) ee -H 2 (f) = H p (r) (V. 19b) 

hir) ee k 2 {r) = k p {r) (V.19c) 

j^^-j^r)^]^ . ( V .19d) 

But in contrast to the magnetic objects, the corresponding electric fields must always 
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differ in sign because both particles carry opposite charges: 

(1 k (f) = - ( %(f) = {b/p) A (f) (V.20a) 

E^f) = -E 2 (f) = E b/p (f) = -V^Uo(r) (V.20b) 

(%(r) = ®k Q (r) = Wp) k (r) (V.20c) 

Uj (r) s -( 2 >j (r) = (fe/p) j (r) . (V.20d) 

Therefore the electric Poisson equation (IV. 16ft for the common electric potential ( b ' p >Ao 
reads for both positronium configurations 

A {b/p) A (r) = -4vra s {b/p) k (r) , (V.21) 



and analogously the magnetic Poisson equation ( IV. 171) becomes for the two different 
configurations 

AA b/p = -4vra s k b/p (r) . (V.22) 

But here the sources {^ b ^ko(r); k b / p (r)} of the electromagnetic potentials 
{( 6 /p)A (f); A b / p (r)} must now be built up by the matter fields ^ a '(p±(r) (1III.2 1[) . First, 



consider the Dirac density w/cq (f) = ipalo^a (IH-161) . which reads in terms of the Pauli 
spinors ( a '(p±(r) as it is shown by equation ( 1IV.46I) . Consequently, the common form 



of the charge density ^ b ^ p ^ko(r) ( 1V.20cl) yields the following common requirement for the 
Pauli spinors: 

%l(f)% + (f) + (V_(r) (1 V-(r) = {2 ¥+(r) {2) f+(r) + { \l{r) (2 \-{r) . (V.23) 

However, since the analogous requirements ( lV.18cl) and (1V.19cl) for the Dirac currents 
differ (or not) in sign, one is led to two different requirements for the ortho- and para- 
case: 

{%l a (V + ( V- B (1 V+} = T{ {2 ¥ + $ ( V + ( V- (2 W} (V-24) 

where the upper/lower sign refers to the ortho/para-case, resp. The solutions for both 
requirements (jVT23]) - (jV33j) are [H 
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ortho-positronmm: (V.zoa) 

W^_(fO = i (2 V-W 4= (6 V-(r) 



para-positronium: - (V. 25 b) 

$± i) . 

1 1 k p 1 1 

The most striking difference between the ortho- and para-states will surely refer to 
the energy E? being carried by those states. For a first tentative estimate of this energy 
difference of both types of states one may tend to the hypothesis that the magnetic forces 
are much smaller than their electric counterparts so that they can be treated as small 
perturbations. Thus one could perhaps first solve the RST eigenvalue problem under 
the exclusive action of the electric forces alone (electrostatic approximation) which would 
yield from equations ( 1IV.20I) and ( 1IV.23I) the corresponding total energy ^Em of the 
positronium system as 

E [T] => (% T] = ( C )4 IV) = (Z 2 1} + Zf 2) ) ■ Mc 2 + 2T kin + 4 e) , (V.26) 

provided all the constraints are validated. Observe here that (for different particles) all 
constraint terms referring to the exchange subsystem must vanish automatically because 
of 5 M = 0: 

N§ ] = = ng } = = nff = = , (V.27) 
while (1IV.5P is left unchanged and for (1IV.2h the last term is zero. Furthermore, 



the rest masses M p , M e of both particles are identical (M p = M e = M) and are denoted 
simply by M. 

Now it is of course very tempting to think that, when all the magnetic interactions are 
neglected, the ortho-states (b) and the para-states (p) will have the same energy 

(e)£(IV) 
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(fV26l) 



(e) E (IV)| = (c) E (IV)| 



>T I," -T | p , (V.28) 

and therefore it is only the magnetic energy term (IIV.29I) which in lowest order 

of approximation is to be computed by means of the electrostatic solution and is then 
responsible for the energy difference between the ortho- and para-states 
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(V.29) 



Observe here that this kind of reasoning is adopted also in the conventional approach 
( IV. 131) where the orbital factor ipi^i, ^2) is determined by the electric forces alone and 
is thought to be the same for all four spin configurations (otherwise the groundstate 
degeneracy would not be four-fold). If those hypotheses like (IV.28I) - (1V.29() were true also 
for the RST case, the ortho/para energy difference in lowest order would appear by means 
of E^ ] f TIV39|) as 
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(V.30) 



and thus by means of flV.18bj) and flV.19bl) 

he 



47ra s 



d f 
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< . 



(V.31) 



This result says that the RST ortho-positronium (conventionally: 1 So) configuration has 
lower energy than the RST para-configuration (conventionally: 3 Si) which is qualitatively 
in agreement with the observations (see the remarks below equation (1V.13j) ). 

Further information about the magnetic energy difference A^-Er (1"V-31j) is gained by 
observing that, in the absence of the exchange interactions, the magnetic fields Hb/ P are 
simply the curls of the corresponding three- vector potentials A b / p , cf. (1III.24al) -( lIIL24bl) 



H, 



b/ P 



V x A 



b/ P 



(V.32) 
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and furthermore the curls of the magnetic field strengths are related to the Maxwellian 
three-currents j b / p by 

V x H b/p = -AA b/p = 4vra s j b/p . (V.33) 

The combination of the latter two equations yields the standard solution of the vector 
potentials in terms of the Dirac currents k b / p as 

A b/p = a s jd^^l (V.34) 



This result may now be used in combination with the magnetic Poisson identity ( 1IV.4I) 
(1IV.5P in order to recast the magnetic energy into the following form: 



4-) = | / d¥ 



(V.35) 



k 1 'A 2 -k 2 ' A-y 

i. e. concretely for both types of positronium configurations 

4 m) | b/p = ±he J d 3 f (k b/p • A b/P ) . ( V.36) 

Thus it is possible to express the magnetic energy completely in terms of the Dirac 
currents as follows 

and consequently the magnetic energy difference A^-Et (IV.29I) reads in terms of the 
Dirac currents 

A^E T = -e*JJ d*rd*r> ^^1±M^ . (V .38) 

This result presents now a very critical point for the intended estimate of the or- 
t ho/para energy difference. Namely, at first glance it seems very tempting to think that, 
because of the smallness of the magnetic interactions as compared to their electric coun- 
terparts, the neglection of the magnetic forces will induce only a negligibly small change 
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of the Dirac currents k a (f) (a = 1, 2) as they are obtained by the electrostatic approxima- 
tion, so that one should be allowed to equate approximately the ortho-and para-currents, 
i. e. 

k b (r) ~ k p (r) = K(r) . (V.39) 

Furthermore, this common current k*(f) could be eventually taken from the solution 
of the electrostatic (and therefore truncated) RST eigenvalue problem. Indeed, such a 
philosophy was put forward in two preceding papers [TSJ [TB] ; and by such presumptions 
the ortho/para energy difference A^-Et (lV-38j) would be finally found in the following 
form 

A^E T -2e 2 [ [ d 3 fd 3 f' k *(r)'K(p ^ (v 4Q) 



see equation (5.113) of ref. [TB] . (Unfortunately, in ref. [TO] the notations for the Dirac 
currents kb and k p are interchanged relative to the present definitions ( 1V.18cl) and ( lV.19cl) ). 

However, serious objections against such kind of reasoning seem adequate because it 
is highly questionable whether the Dirac currents k a (a = 1, 2) are really (or almost) 
identical for the ortho- and para-system, even in the electrostatic approximation! The 
crucial point in this context is namely that two strongly interacting fermions do form a 
bosonic system when being considered as a whole, and in RST this bosonic character of the 
compound system becomes transferred to the individual particles which as a consequence 
do then occupy bosonic one-particle states. In general, such exotic RST states carry then 
an integer eigenvalue of total angular momentum J z = © jj ^ (the z-axis is adopted 
to be the symmetry axis of the two-particle system): 

J z ^ a (f) = ( Ji+) < a V + (r)) © (Ji- )(a V-(r)) = 3Aa{f) t 

(V.41) 

j, = 0,±l,±2,... 

But since the one-particle wave functions ipair) are Dirac four-spinors, this integrity 
condition (IV.41|) can hold only at the price of non- uniqueness of the spinors ip a {^), see 



the discussion of this effect in ref.s [T3l [T6] . Thus for each one of the two particles of 
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the positronium groundstate one expects j z = (-w singlet) for the para-state and j z = 
0, ±1 (~-> triplet) for the ortho-system. Now it is well-known in (conventional) quantum 
mechanics that wave functions ip a with different eigenvalues of angular momentum can 
display a rather different spatial pattern and this difference will then be transferred also 
to the currents k a which are (more or less uniquely) determined by the wave functions ip a , 
see the Dirac currents (III. 161) . Therefore one concludes that, in contrast to the equality of 
the densities ^k (r) (lV.20c[) . the presumed equality (lV.39j) of ortho- and para-currents 
kb and k p can not be correct in general, not even in the electrostatic approximation, cf. 
fig. 1 vs. fig. 6; and as a consequence the result ( 1V.38I) is to be prefered in place of ( IV. 461) . 
Clearly, this conclusion forces us now to inspect more rigorously the spatial pattern of the 
ortho/para wave functions if} a (r) and of the associated curents k a (f). 



3. Positronium Eigenvalue Problem 

The postulated positronium constraints (lV.18al) -( lV.20dl) for both gauge field modes 
(a = 1, 2) imply the corresponding constraints (1V.23I) -( TV.24I) for the matter modes ^ip±{r) 
(a = 1,2) so that both kinds of constraints are consistent with each other and thus admit 
to deal with only one single matter field (p± for the considered two-particle system: ^<p±(r) 
( 1V.25al) for the ortho-configuration and ^(p±(r) ( 1V.25bl) for the para-configuration. But 



with this arrangement there arises now a serious problem; namely the question whether 
each of these two links (1V.25ap - (1V.25bj) is additionally consistent also with the two-particle 



eigenvalue equations (1III.22al) - (1III.22bl) ? The latter consistency requirement says that any 
of the two equations ( lIII.22al) and (1III.22bl) must collapse to the same equation for the 
Pauli spinors ^ a V±(^*) ; namely either for ^tp±(r) in the ortho-case or for < - p V±(r) m the 
para-case. Clearly, the situation with ortho-positronium ( lV.25al) is simpler and an explicit 
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calculation yields for the one residual Pauli pair of equations [H]: 

i ff • (V + il b ) % ± {f) ~ {b) A ■ = c ■ W M0 

« (V.42) 

(ortho-positronium) . 

Here both mass eigenvalues M a are assumed to differ merely in sign and the rest masses 
of electron (M e ) and positron (M p ) are equal 

M x = -Mi = -M* (V.43a) 
M e = M p = M . (V.43b) 

However, the situation with para-positronium ( 1V.25bj) is somewhat more intricate but 
nevertheless admits also both Pauli pairs of equations (lIII.22ap - (1III.22bl) to collapse onto 
one single pair: 

iff • (V - il p ) (P V±(^) - {p) Ao ■ {p) if T (r) = U ' U c • {p) <p T (r) 

h (V.44) 

(para-positronium) . 

The deduction of this para-equation is not only more difficult than its ortho-counterpart 
(IV. 421) but it yields also an additional constraint upon the para-spinors ^(p±(r) (]V.25bj) . 
i. e. 

j( ±)(p V±(rO = , (V.45) 
which itself is a special case of equation (IV.41I) 

J" 2 ^a(r)=0. (V.46) 



For both cases (IV.42I) and (IV.44j) it was assumed here that both the three-vector potentials 



A a and the Dirac currents k a are of azimuthal character, i. e. 



A a = ■ e (V.47a) 
4 = (a) fc -e . (V.47b) 



72 



But observe that the ortho-case ( 1V.42I) must not be associated with such a restriction 
as ffVTIH]) ! Moreover, it may seem that the definitions flV.lSah - flV.lQdj) for the ortho- 
and para-configurations are rather arbitrary; but actually we do assume here that such 
configurations are the only localized solutions of the two-particle system (lIII.22aP - (lIII.22bl) 
if both rest masses are identical (i. e. M e = M p = M). Clearly, the proof of this hypothesis 
needs an extra discussion. If the hypothesis turned out to be wrong, the RST predictions 
would come in conflict with the observational data which yield a one-particle spectrum for 
positronium [llj. Observe also that, when the magnetic interactions are neglected (Ab = 
A p = 0), the ortho-equation ( IV. 42j) becomes identical to the para-equation ( IV. 441) ; but 



since the angular-momentum restriction (IV.45I) does hold only for the para-case (also in 
the electrostatic approximation) one expects that ortho-positronium has a richer spectrum 
(see below). 

The vanishing eigenvalue of total angular momentum J z (1V.46I) is a very unorthodox 
feature of a Dirac particle's quantum state because one usually expects half-integer eigen- 
values. However, as mentioned above, the bosonic character of a two-fermion system is 
transferred in RST to the individual constituents so that they can individually acquire 
bosonic properties. This amazing effect may be elaborated in RST by an adequate choice 
of a spinor basis relative to which all the Pauli spinors ^(p±(r) are then to be decom- 
posed. Our choice of such a basis does not directly refer to the spinors £/' m (1V.6al) - (1V.6cl) 
but rather to a certain modification of them HEI 



J+) = (V.48a) 

= e** • CT 1 (V.48b) 
J+) = e -iH . g>h (V.48c) 
J-) = . Q-l ; (V.48d) 

where the boson number o is a real number with range — | < 6 < |. This number measures 
the bosonic character of a Dirac particle which becomes more obvious from the action of 
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Ji ± ' ) upon those basis spinors 



j| + M +) = -(6-7,)M +) (V-49a) 

J^J-^ib-^h-J^ (V.49b) 

J^uJ^ = -(b- ^)h-u[ +) (V.49c) 

J^J-^ib-Uh-J^. (V.49d) 



Whenever the boson number 6 is zero, the new basis {4?} flV.48ajl - (IV.48djl coincides 



with the original basis {Cq 2 i 2 } ( lV.6ap - (lV.6cP and one deals then with a purely fermionic 



basis, since the eigenvalues j z of Jz^ are half-integer. However, for S = ±| the eigenvalues 
of J z are integer and one has a purely bosonic basis (the case — \ < 6 < \ with t> ^ is 
something between a fermionic and bosonic basis and needs an extra discussion). 

But once a spinor basis has been selected in such a way, one can decompose the Pauli 
spinors ( b / p )y?±(r^) for the ortho- and para-system in the usual way 

(^V + (f) = ^>JZ + • 4 +) + • w<-> (V.50a) 

^tp_{f) = -i \ (b/p) R- ■ 4 +) + ^<S_ ■ J- ] ) . (V.50b) 



Strictly speaking, such a spinor coordinatization should first refer to the original single- 
particle states ^Lp±{f) (a = 1,2); i.e. one puts 

(°V+(r) = {a) Tl + ■ 4 +) + {a) S + ■ J ~ ] (V.51a) 
W^_(f) = -i |( a >^_ • ^ +) + (a) S_ • wj -3 } , (V.51b) 

but on account of the ortho- and para-identifications ( 1V.25al) -( 1V.25bl) the wave amplitudes 
( - a ^JZ±, ( a >S± are also to be identified in a certain way, namely for ortho-positronium (!V.25aP 
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{l) TZ + = -i {2) TZ + = {b >Jl + (V.52a) 

(%+ = -* (2) S+ = (V.52b) 

(1 >K_ = z (2) K_ = (6 >K_ (V.52c) 

«<S_ = * (2) <S_ = (V.52d) 



and similarly in the case of para-positronium ( lV.25bl) for boson number 6 = | 



2' 



= (2) S + = (V.53a) 

= _( 2 >^ + = (p)5 + (V.53b) 

= _(2) 5 _ # (p>^_ ( V .53c) 

(1) <S_ ee ( 2 >^_ = ^<S_ . (V.53d) 

Observe here for the latter case of para-positronium that the requirement of vanishing 
angular momentum (1V.46I) can be attained in the simplest way by (i) choosing for the 
boson number S the special value 8 = ~, which puts the eigenvalues j z of J z in ( 1V.41I) to 
zero (j z =>- 0), and (ii) adopting the wave amplitudes Wft±, ^S± as being independent 
of the azimuthal angle 0, i. e. 

{p} Jl±(r } 0, 0) (p) i?±(r, #) (V.54a) 
ip) S± (r, 0, 0) ^S± (r, i?) . (V. 54b) 

The parametrizations (lV.51al) - (1V.51b|) of the Pauli spinors ( 6//p )</} ± (f) in terms of the 
(complex- valued) wave amplitudes ^ p >Jl±, ^S± can now be inserted in the two-particle 
Pauli eigenvalue equations for ortho-positronium flV.42j) and para-positronium (1V.44I) in 



order to yield the corresponding eigenvalue equations for these wave amplitudes; i. e. for 
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ortho-positronium 



+ e 



2i{ 



or r o<p r 

r ov r od> 



M + M, 



(V.55a) 



^ _ 1 . ^± + I . (% + _ ■ (%_ + sin <% ■ 

ar r am r 



_ -2i{%-\)4> 



1 a (b) fc + cot^ ^ (fe> ^ . 

r <9$ r + cM 



-cosi?^-^ 



M + M, 



(V.55b) 



dr 



r a© r 



+ cos ( fo U, • ^ *> = c ■ {by R+ (V.55c) 



^ + - • ^ + — ■ + ( 6 U ■ % + - sintf(% ■ 
ar r a<p r 

2i(h i u il a( 6 >K_ cottf r . , ftVo .0 

+ e" 2 ^"^ . ^ — - + 6 ■ ( 6 >K_ + i 

r aw r a<p 

- cos ■ 1 = c ■ . (V.55d) 



Here, the value of the boson number S has been left unspecified for the sake of generality; 
but it is clear from equations flV.49aP - flV.49dl) that the right value for ortho-positronium 
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will be 8 = — |. 

Quite analogously, for para-positronium one will adopt the value 8 = |, which puts 
the eigenvalues of J z to zero, cf. (IV. 45ft ; and additionally one adopts the wave amplitudes 
( a fo±, ^S± (IV.54ap - flV.54bp to be real-valued functions ( {a) R±, {a) S±) over three-space. 
For this arrangement, the two eigenvalue equations (1V.44P read then in terms of those 
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(real-valued) wave amplitudes {p) R±, ip) S± for para-pos itronium 

° ll " n - 1 ■ {P) R + - (p U ■ WiL. + ( % [sintf • - costf • + 1 ^"'^ 



<9r 2r r d'd 

2r + /I 

(V.56a) 

^ + I ■ - V>Ao ■ ^S. - H4, [costf • Wfl^ + sintf • - 
or 2r r ov 



cot*. W „ = M + M, c . Ws 



2r ! ft 

(V.56b) 



Or 2r r ov 



(v) S __ M - M * c (V) R 



2r h 

(V.56c) 



0>I " S ~ + A ■ + ip) A ■ (P)S + + ( % [cost? • WiL. + sintf ■ W&_] ' ^ ^ 



dr 2r L J r ch9 



2r h 

(V.56d) 

Both eigenvalue systems, i. e. the ortho-system ( lV.55al) - (1V.55dl) and the para-system 
(1V.56al) - (1V.56dl) provide now the basis for a thorough inspection of the RST multi- 
plet and degeneracy structure. Observe here that the main feature of this structure 
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is the ortho/para dichotomy for which RST supplies a very natural basis, namely the 
(anti-)parallelity of the magnetic fields. In contrast to this rather dynamical criterion, 
the conventional theory puts this ortho/para dichotomy on the basis of the composition 
rule for the angular momenta of both particles which is rather a kinematical criterion 
than a dynamical one! In this context it is worthwile to reconsider the kinematical con- 
sequences of the dynamical requirement concerning the (anti-)parallelity of the magnetic 
fields. Clearly, the parallelity or antiparallelity of the magnetic field strengths H a (r) 
(a = 1,2) entails the antiparallelity or parallelity of the corresponding Dirac currents 
k a (r) (or Maxwell currents j a , resp.) because the latter are the curls of the magnetic 
fields, cf. (1V.33p . But on the other hand, these currents are themselves generated by 
the wave functions ip a {r), cf. (III.17a[) - flII.17e[) ; and therefore the currents can transcribe 
those (anti-)parallelity requirements for the field strengths H a to the wave functions, i. e. 
ultimately to the wave amplitudes (°vl±, ( a >S±. More concretely, when the Dirac cur- 
rents k a are decomposed with respect to the basis triad {e r ,e,j,e^} of the spherical polar 
coordinates {r, 0} 

K = {a) k r ■ e r + Wfo • e„ + ( % • e^ , (V.57) 
the current components ^k r , ^k$, ^'k^ are linked to the wave amplitudes through 

Wfc p = JL . ( a >R_ + ( a £ + . W&_ - ^yjl_ • ( a >JZ + - • (V.58a) 

ia) h = ~ {e*M>* ■ C [a) - e~^~^ ■ C (a) } (V.58b) 

(c ( «) = {a} k + ■ (a) 5_ + (a) 7L ■ {a) s + ) 

|e 2 ^~^-C (a) +e- 2 ^-^^C (a) } . (V.58c) 



47T 

COS?? 



47T 

Now one can test here whether the ortho-antiparallelity ( 1V.18cl) and the para-parallelity 
(IV. 19cp of the Dirac currents k a (f) are really validated by both the ortho-identifications 
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( 1V.52al) -( 1V.52dl) and the para-identifications ( lV.53al) -( 1V.53d[) . For the latter para-case 
one finds 



Wk,. = -^k r (V.59a) 
Wfo = -( 2 )fo (V.59b) 
<% = + ( % . (V.59c) 



But obviously it is only the azimuthal component (1V.59cP which is consistent with the 



para-identification (IV.19C0 . but not the radial and longitudinal components (IV.59ap 



(1V.59bl) . If one therefore insists on the three-vector identification (1V.19cjl . one is forced 
to let vanish the radial and longitudinal components of the Dirac para-current k p : 

Mkr = {2) K = {p) k r = (V.60a) 
Wfo = {2) h = {p) h = . (V.60b) 

As a result, the Dirac para-current k p must encircle the z-axis in the planes z = const 

^, = W V^> (V.61) 



cf. the anticipation flV.47bj) . Moreover it is easy to see that both requirements of vanishing 



current components ( 1V.60al) -( 1V.60bl) can be satisfied for t> = | by simply taking all the 
wave amplitudes {ay Jl±, {a) S± to be real- valued ( (a >ft ± => (a) R±; (a) <S± => (a) S±), as it was 
already anticipated for the para-system (lV.56aP - (1V.56dj) . 



However, concerning the ortho-system, one finds that the current components ( 1V.57I) 
for the ortho-value of o (i.e. S = — |), together with the ortho-identifications (1V.52ajl - 
( lV.52dl) . appear in the following form 



W kr = - {2) k r (V.62a) 
Mfo = -( 2 )fc^ (V.62b) 
= . ( V .62c) 
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This means that the ortho-antiparallelity of the Dirac currents ( 1V.18cl) can always be 
attained by the present ortho-parametrization of the wave functions ipa{r), with no ad- 
ditional restriction as in the preceding para-situation. Fortunately, this admits to evoke 
here the toroidal configuration {k^ = 0, k# ^ 0, k r ^ 0) of ref. [T7j which may be adopted 
here in order to provide a representation for that ortho-state with vanishing j z (= 0). 
However we will subsequently restict ourselves to that rotational symmetry around the 
z-axis as given by equations flV.47ap - flV.47bj) . for the discussion of both the para- and the 
ortho-system. In this sense, we specialize down the general Poisson equation ( IV. 221) for 
the vector potentials Aj,/ p to their 5*0(2) symmetric form 

(typ)4 

A(6/ % - -r^h = " 4 ™s ■ {h/p) h (V-63) 

r z sin 17 

which together with its electric counterpart ( IV. 2 II) completes the present eigenvalue equa- 
tions of the ortho- and para-type to a closed system. 
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VI Energy Difference of Ortho- and Para-States 

Once the emergence of the ortho/para dichotomy in RST has been clarified, there 
naturally arises the question of degeneracy of the corresponding ortho- and para-states. 
But here it should be clear that, corresponding to the present state of the art, this question 
can be discussed only from a rather qualitative point of view because exact (analytic or 
numerical) solutions of the RST two-particle eigenvalue problem are presently not known. 
Nevertheless, the knowledge of the exact form of the RST energy functional E[t] ( 1IV.20I) 
admits to draw same general conclusions concerning this question of degeneracy. 

First, it will be necessary to define some order of approximation which corresponds 
to the conventional situation (jV.13p . where the singlet and triplet states carry the same 



energy. As a matter of course, this first approximation step will be essentially the same 
in RST as in the conventional theory, namely the neglection of the magnetic interactions 
(-w "electrostatic approximation") . This neglection of the magnetic forces, however, does 
not yet entail the mathematical equivalence of the eigenvalue equations for the ortho- and 
para- configurations; indeed, there persists a crucial kinematical difference: whereas the 
current of the para-states u para- current" k p ) must always encircle the axis of rotational 
symmetry (i.e. the z-axis) and thus realizes a cylindrical structure, see equation flV.6ip . 



the "ortho-current" kb can obey this cylindrical geometry only for non- vanishing angular 
momentum (j z = ±1), but for vanishing angular momentum (j z = 0) it is required to 
undergo a toroidal geometry (previously discussed in ref. |17j). In the present paper, we 
restrict ourselves to the cylindrical structure (i.e. j z = for the para-states and j z = ±1 
for the ortho-states); but even within this subclass of the cylindrical configurations the 
link between the charge densities ( 6//p )/c and the wave amplitudes W p ^Jl±, ^S ± looks 
rather different: for the para-configurations, the essential angular variability of ^ko(r) 
goes like (sin??) -1 , cf. ( 1VI.17I) below; but for the ortho-densities ^ko (f) there occurs 



a certain subdichotomy, see figure 5, namely in such a way that both forms of ortho- 
matter with density (f) are pressed off the rotational axis upon which they become 
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zero ( ( ^oU,. = 0). 

Clearly, this specific distribution of ortho- and para-matter is responsible for the differ- 
ence in angular momentum (i.e. j z = for the para-distribution being concentrated close 
to the rotational axis (~ (sin??) -1 ); and j z = ±1 for the ortho-distribution concentrated 
off the z-axis (~ sin??)). However such a redistribution of the charge densities does im- 
ply also a corresponding alteration of the electrostatic interaction energy E^' so that the 
ortho- and para-configurations would always carry different energy and the ortho/para 
degeneracy would be broken. Thus one is forced to conclude that, in order to encounter 
the degeneracy phenomenon in RST, it is not sufficient to neglect the magnetic interac- 
tions as is the case with the conventional theory, see the discussion of equation f]V.13|) . 
Additionally, one has to resort to a further approximation, i. e. the spherically symmetric 
approximation, which consists in assuming the electrostatic interaction potential ( 6 / p l4o(f) 
(not the wave amplitudes!) to be spherically symmetric: ( b / p 'Ao(r) => ( b / p 'Ao(r) (r = \r\). 
And indeed, it is just this assumption which yields the complete mathematical equivalence 
of the eigenvalue problems for the ortho- and para-configurations, cf. the para-eigenvalue 
equations (lVI.23al) -( 1yL23bl) to their ortho-counterparts (I VI. 1 15aH - H VI. 1 1 5bjl . The (ap- 
proximately) common energy spectrum is sketched in fig. 2. The RST energy levels of 
this common spectrum do agree with their conventional counterparts (1VI.88I) with devi- 
ations of magnitude of some 10% (see the table on p. 1 1Q7D . This is a remarkable result 
because it is based upon a very simple trial function $(r) ( 1VI.55I) with only two variational 
parameters! 



1. Para-Positronium 

As a matter of course, we will first turn to the simplest field configuration which 
is surely the para-system with its singlet character (5 = |). The corresponding energy 
eigenvalue problem is defined by the mass eigenvalue equations (lV.56al) - (1V.56dl) together 
with the Poisson equations (IV. 211) and (IV.63I) for the electric potential ^A (r, $) and 
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the azimuthal component ^A^(r, $) of the three- vector potential A p {f) } cf. equation 
( 1V.22I) . Here it should be self-evident that it is very hard (if not impossible) to find exact 
solutions for the para-system so that one is forced to resort to an approximative treatment. 
Especially for the para-groundstate one expects the highest possible symmetry of the 
solutions; but this is surely not the 5*0(3) symmetry with respect to the origin (r = 0) 
because this would imply hedgehog-like three- vector fields E p (r) and H p (r), k p (f). This 
could be accepted for the electric field strength E p (r) (~ r), but not for its magnetic 
counterpart H p {f) because the latter would imply the existence of magnetic monopoles 
(to be excluded here). 

Electrostatic Approximation 

Hence for a first rough estimate of the para-levels one may resort to the electrostatic 
approximation where all objects of the magnetic type are to be neglected. Here the hy- 
pothesis is that the magnetic forces are much weaker (than their electric counterparts) and 
therefore they can hardly influence the charge densities ^jo(r) whose spatial distribution 
will be dominated by the electric forces. Since the corresponding electrostatic interaction 
energy (E^) is itself also much greater than its magnetic counterpart (E^), the influ- 
ence of the magnetic interactions on the total energy E^ appears therefore to be negligibly 
small. Clearly, such a picture of the magnitude of the intra-atomic gauge forces seems 
to provide the legitimation for neglecting the magnetic interactions altogether. However, 
there may exist certain regions in three-space (here the z-axis) where the magnetic forces 
become (almost) infinite and thus enforce a distinct redistribution of the electric charge 
densities ^jo(r) so that the electrostatic energy E^ may become modified considerably. 
Evidently, for such a situation the electrostatic approximation is not adequate; but we 
will assume here that this magnetically induced shift of the electric energy E$ will not 
spoil too much the results. However, it may well be possible that through the omission 
of such an energy shift the lower bound of the original energy functional E^ becomes 
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underestimated by the electrostatic approximation (see the discussion of the positronium 
groundstate in ref. p3j). 

In any case, by this neglection of magnetism the original para-system (!V.56aP -( 1V.56dl) 
reduces to 

^± + I^ft - »A, . WJL = ^ c • MJL (VI.la) 
or r ov n 



^± _ I^|± _ (PU„ . <»s_ = (VI.lb) 
or r ov n 



^ + I . WJL _ i ^ + . « = ^ c . WAf (vi.lo) 
or r r ov n 



^ + I . WS . + i^f + <>U • <*>S + = ^ c ■ «S + . (VLld) 
or r r ov a 

Here the modified wave amplitudes &R ± , ^S± are defined in terms of the original real- 
valued amplitudes {p) ll± {p) R±) and (p) <S± (=>■ (p) S±) (TV.54ajl - flV.54bjl through 



^R± = Vr sin d ■ {p) R± (VI.2a) 



{p) S ± = Vr sintf • ( p ^ ± . (VI.2b) 

Being satisfied for the moment with that electrostatic approximation, the relativistic 
energy functional Em (1IV.20P becomes simplified for the two non-identical particles which, 
however, must be assumed to occupy two physically equivalent quantum states (1V.25bl) . 
This implies that both mass renormalization factors Z/ a -\ become identical 



Zfi) — 2(2) -j- %v (VI.3) 



with 



(d 2 r = r dr d"d) . 
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(VI.4) 



Furthermore, the adopted physical equivalence of both particle states implies also that 
the kinetic energies T kin(a ) ( 1IV.26I) -( 1IV.27I) must turn out to be identical 

T kin(1) = T kin(2) = Wfwn . (VI.5) 

Here the common one-particle kinetic energy T kin splits up into a radial part (T r ) and a 
longitudinal part (T#) 

^T kin = + (P>f w (VI.6) 

with the radial part being given by 

°9 M = h 4 let ■ ^ - ■ ^ - ir>A+ - " >k - + «5_ • 



4 ./ 1 dr dr r dr 



dr r 

(VI.7) 



and the longitudinal part by 



4 ./ r I " Or) r <9tf <9tf 

(VI.8) 



And finally, the electric constraint Nq^ ( 1IV. lh -( |IV.2j) is simplified to Nq 



= e£ ] - M^c 2 = (VI.9) 



with the electrostatic interaction energy being given by 



4 e ) = / d 3 r\\V^A \\ 2 (VI.10) 



47TOJ 

and its mass equivalent by 

M(e) c 2 = _^ J dV-^U {(Wi? + ) 2 +( (p) S + ) 2 +(^_) 2 +( (p) 5_) 2 } . (VI.11) 

Thus, putting now together all energy contributions yields for the relativistic energy 
functional (1IV.20P with (1IV.23h in the electrostatic approximation 

E [T] E {T} = 2Mc 2 ■ Z 2 V + 4^>f kin + E^ + 2A D ■ A> D + A§ } ■ • ( VL12 ) 
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Here the Lagrangean multiplier Ad is given by the sum of the individual multipliers 

A D = A D( i) + A D{ 2) = (Mi - M 2 )c 2 = -2M*c 2 , (VI.13) 



cf. 0IV.22aP ; and both normalization constraints A^D( a ) (HH.49I) collapse to a single one 
just on account of the physical equivalence of both quantum states 



(2) ^ 



(VI.14) 



with 



iVi 



b = |y d 2 r {( (p) i? + ) 2 + { ip) S + ) 2 + («iL) 2 + (^5_) 2 } -1 = 0. (VI.15) 

Indeed, after all it is now a very satisfying consistency check to deduce the relativis- 
tic para-system (IVI.laft - ljVI.ldl) from the principle of minimal energy (1IV.21I) with the 
functional Em being given by its present electrostatic approximation -E{t} (IVI- 1 2D . But 
for the subsequent construction of (approximate) solutions one should keep in mind those 
boundary conditions which are needed for the variational deduction: 



{p) R+ 




= lim 

r- >0 




= 


(VI. 16a) 


ip) R + 


• 


19 = 7T 


ip) R + ■ {p) S_ 


i9=0 


= 


(VI. 16b) 


(p)r_ 


■ {p) S + 


t9 = 7T 


(P)R_ . (P)§ + 


i9=0 


= . 


(VI. 16c) 



Finally, in order to close the truncated eigenvalue system (IVI. la[) - (1VI. ld|) . one has to add 
only the electric Poisson equation (IV. 211) for the potential ^Aq 



A^A = -4vra s {p) k 



r sin^ 



(VI.17) 



since the azimuthal magnetic potential ^'A^ is to be neglected in the electrostatic ap- 
proximation. Of course in the true spirit of the principle of minimal energy, this Poisson 
equation (IVI. 171) must be deducible also from the approximative energy functional -E{t> 
( IVI. 121) . namely through its extremalization with respect to the electric potential ^Aq\ 
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Spherically Symmetric Approximation 

Despite the simplifying effect of the electrostatic approximation, the resulting eigen- 
value system ( 1VI.lal) -( 1VI.f dl) plus the Poisson equation ( 1 VI. 1711 is still a too complicated 
system of coupled equations in order to find the exact solutions thereof. Consequently, 
it becomes necessary to think of a further approximation step which should refer to the 
gauge potential (p U (r,$) solution of the electrostatic Poisson equation (IVT.17p . 

Indeed, what obviously does complicate things with that equation is the fact that its so- 
lutions will depend in general on both variables r, § even if all the wave amplitudes 
^R±, ^S± would be assumed to be spherically symmetric (i.e. not being dependent on 
the spherical polar angle $). Therefore it suggests itself to adopt also the electric potential 
( p l4o(r, i?) to depend exclusively on the radial variable r but not on the angle i.e. the 
approximation is the following (" spherically symmetric approximation") 

(p U (r,tf) [p U (r) . (VI. 18) 

But observe here that our present assumption of spherical symmetry ( 1VI.18I) does 
concern only the gauge potential ^A , not the wave amplitudes ^R±, ^S± themselves! 
Indeed, we admit these amplitudes to depend on the angle in the following way: 

w &tM) = f R±(r) • /r(0) (VI.19a) 
<P>£t(r,0)= S±(r) • f s {0) . (VI. 19b) 

The angular functions /r(i?) and fs{$) of this product ansatz must now be determined 
in such a way that the four partial differential equations (rVI.lall - flVI.ldjl become reduced 
to ordinary differential equations with respect to the radial variable r. It is easy to see 
that this goal can be attained by demanding that both angular functions and fs($) 

obey the following coupled system 

«tf) 



t v ■ f s (#) (VI.20a) 
-l v ■ f R (#) , (VI.20b) 
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where i-p is some real- valued constant. 

Actually, by this demand the original system ( 1VI.lal) -( jVI.ldl) is simplified to the fol- 
lowing coupled set of ordinary differential equations: 

- ■ 27 + (r) - ^U (r) ■ ft_(r) = M + M * c • ^(r) (VI.21a) 
dr r n 

d -^[± - . ft +(r ) - bU (r) ■ 2?_(r) = M+ fc M> c • S_(r) (VI.21b) 
dr r n 

d L^z!A + ±.> R M + ^. > S _(r) + ^Uo(r) • ft+(r) = M ~ M * c • ft + (r) (VI.21c) 

+ - • S_(r) + ^ • fc_(r) + [p U (r) • 'S + (r) = ^L^L c . > S+ ( r ) . (VI.21d) 
dr r r n 

Clearly, such a system of ordinary differential equations needs further discussion [13] since 
the same eigenvalue M* emerges in each equation! In view of this fact one will suppose 
that the four solutions 'R±(r), 'S±(r) can not be too different and thus one tries the 
identifications 

'R±(r) = 'S±(r) = $±(r) . (VI.22) 

And indeed, by this assumption the four eigenvalue equations (IVI.21af )- (lVI.21dj) become 
contracted without any contradiction to only two equations 

^ - * ■ $ + (r) - W Mr) . $_ (r) = ^+Ml c . $_ (r) (VL23a) 
dr r a 

^ + <*±± ■ *_(r) + M Mr) . $ +(r) = c . $ + (r) . ( V I.23b) 

These two equations, together with the spherically symmetric approximation of the elec- 
tric Poisson equation (jVI.17p . constitutes now the eigenvalue problem of para-positronium 
in the spherically symmetric and electrostatic approximation. 

The spherically symmetric form of the Poisson equation will readily be discussed, 
but first one should think about the possible values of the constant ip appearing in the 
angular system (1VI.20ap - (lVI.20bj) . For this purpose, one converts this first-order system 
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( lVI.20al) -f jVI.20bl) into a second-order system by differentiating once more 



d 2 f R ($) 
d 2 fs(#) 



-Up ■ fniti) (VI.24a) 
-4 ' fs(#) • (VI.24b) 



d$ 2 

The solutions of this standard problem are well known and are written here as 

f R (#) = cos (£ v . <&) (VI.25a) 

f s (ti) = -an(l r -0), (VI.25b) 

or alternatively 

f R (#) = sm(£ v ■ 0) (VI.26a) 

f g (&) = cos(£ v ■ 0) . (VI.26b) 

But with the general shape of the wave functions being known, i. e. 

(p '*R±(r,-&)=$ ± (r)-f R (0) (VI.27a) 
Vs ± {r,#)=$ ± {r)-fs{&) , (VI.27b) 

one can reconsider the former boundary conditions ( lVI.16bl) -( lVI.16c l ) which yields the 
requirement 

sm(2n£ v ) = , (VI.28) 
and this can be obeyed by fixing the value of £p by the (half-) integers 

*^ = °- (I) • 1 - (i) ■ 2 - (i) ■ 3 - - ■ (VL29) 

In this context it is also interesting to remark that, through this disposal (1 VI. 291) of £p, 
the para-current k p (IV. 611) remains finite on the z-axis (i.e. = and i? = 7r). Indeed, 
substituting the present product form flVI.27ap - flVI.27bp of the wave amplitudes back into 
the general form flV.58cl) of the Dirac current, one finds 

r \ $ + (r)-$ (r) sin \(2£ v + 1) • ■&] 

r, = • q J • VI.30 

2irr sm 
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(The positive/negative sign refers here to the cases (VI.25/26)). This result says that the 
para-current k p remains finite but nevertheless is singular on the z-axis; however, as we 
will readily see, such kind of singularity does not spoil the energy functional, even if the 
(presently neglected) magnetic interactions are included. The alternating flow directions 
of the para-current k p around the z-axis are sketched in fig. 1. 

The energy functional -E{t} ( j VI- 121) is needed now in the context of the spherically 
symmetric approximation because it is not obvious how the corresponding approximate 
form of the electric Poisson equation ( IVI. 171) should look like. In order to deduce this 
from the principle of minimal energy, one substitutes those spherically symmetric ansatze 
for the potential ^A (1VI.18P and for the wave functions (IVI. 19a|) - (1VI. 19bj) into the elec- 



trostatic approximation E^y fIVI. 121) of the original energy functional Em and thereby 
finds the following results for its constituents: 

(i) the mass renormalization factor Zj, (1V"I.4j) emerges as 



POO . - 

Zl^Z\ = -^ drr{($ + (r)) 2 -(<MO) 2 } (VI.31] 



(ii) next, the radial and longitudinal kinetic energies (p) T (r) (EED and Wf w flYL8} 
become simplified to 

M ~ rol - 7T _ f°° , f ~ . , d<f>+(r) ~ . , d§Jr) $+(r) ■ $Jr) 
{p) T ir) lp] T (r) = ^ ftc jf drr U(r) • —)±± - $ + (r) • - -^U_ 

(VI.32a) 

poo 

(P)f m Wf w = -- hd v j dr 4>+(r) • 4>_(r) (VI.32b) 



2 



o 



(iii) the electric Poisson constraint Nq (IVI.9I) is left formally unchanged, i. e. 

N§ = 4 C] - M [c] c 2 = 0, (VI.33) 
with the electrostatic interaction energy (IVI. 101) appearing now as 

91 



and analogously its mass equivalent ( 1VI.11I) 

M {e) c 2 M [c] c 2 = -^hc c/rr [p U (r) • |($+(r)) 2 + ($_(r)) 2 | (VI.35) 

(iv) and finally the normalization constraint (1VI.14P reads now 

N D ^N^ = ^J^ drr j($+(r)) 2 + ($_(r)) 2 j - 1 = . (VI.36) 

Putting now together all four contributions (i)-(iv) yields the desired spherically sym- 
metric approximation say) of the electrostatic approximation -E{t} ( 1VI.12I) as 

E m = 2 Mc 2 ■ Z% + 4 [pl T kin + E§ + 2A D • + A§ } • N§ . (VI.37) 

Here it is again a nice exercise to convince oneself of the fact that the spherically symmetric 
eigenvalue equations (lVI.23ap - (1VI.23bj) actually do arise as the extremal equations due 



to the present functional flVI.370 ! But furthermore, the extremalization of i?r$i with 
respect to the electric potential H4o yields the wanted spherically symmetric form of the 
electric Poisson equation (IV"I.17[) : 

d^Ao + 2 ^Uo = _ E Qs (^(r)) 2 +ft-(r)) 2 ^ 
dr 2 r dr 2 s r 

with the standard solution 



2 



[p] Mr) = yJ— ~ W^Wji • ( VL39 ) 

As a consistency check, one easily deduces hereof the asymptotic Coulomb form (1III.60P 
of the electric potential H4o(r) 

lim [p U (r) = — (VI.40) 

1 — >oo T 

by simply regarding the normalization condition in its spherically symmetric form (1VI.36I) . 
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Summarizing, the combination of the electrostatic and spherically symmetric approxi- 
mation consists of the mass eigenvalue equations ( lVI.23al) -( 1VI.23b[) plus the electric Pois- 
son equation f]VI.38|) . If the (exact or approximate) solutions of this eigenvalue problem 



could be attained, one would substitute these solutions into the associated energy func- 
tional E[$] ( 1VI.37I) and would thus obtain the corresponding energy spectrum (i. e. the 
singlet spectrum) of para-positronium. But although the corresponding energy levels 
would be still of relativistic nature, the corresponding relativistic effects would not be 
of much help! The reason is that the neglected magnetic and anisotropic effects are in 
general of the same order of magnitude as the relativistic effects, so that the numerical 
usefulness would be spoiled in any case. Therefore it is more meaningful to finally do the 
last step of the approximations, namely by turning to the non-relativistic domain. 

N on- Relativistic Approximation 

Though the two preceding approximations (i. e. the electrostatic and the spherically 
symmetric one) did considerably simplify the positronium eigenvalue problem, a further 
rather technical simplification (i. e. the transition to the non-relativistic domain) is highly 
welcome when one is satisfied with a first rough estimate of the positronium spectrum. 
The point of departure for such a purpose is the eigenvalue system ( lVI.23ap - (lVI.23bl) in 
combination with the Poisson equation (1VI.38I) for the electric interaction potential ^Aq 
which has already been assumed to be spherically symmetric. Now in order to deduce 
the non-relativistic approximation of those eigenvalue equations flVI.23aP - flVI.23bp one 



resolves (approximately) the first equation (IVI.23al) with respect to the "negative" wave 
amplitude $_(r) in terms of the "positive" amplitude $+(r), i.e. one puts 

and then one substitutes this into the second equation (1VI.23bp . This then yields the 
following Schrodinger-like equation for the non-relativistic approximation ($(r), say) of 
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the "positive" amplitude $ + (r): 

- w + ^} + m% •<-> - ^ ' *m - * ' *« <™ 2 > 

where the Schrodinger eigenvalue £7* emerges as the difference of the mass eigenvalue M* 
and the rest mass M 

E* = (M* - M)c 2 . (VI.43) 

Finally, in order to close this eigenvalue equation (1VI.42I) . one has to add the non- 
relativistic approximation of the electric Poisson equation (1VI.38[) . This, however, is a 
simple thing because the general form of the Poisson equation survives the transition to 
the non-relativistic domain so that one merely has to resort to the non-relativistic approx- 
imation for the charge density ^ko(r). But since the latter approximation merely consists 
in neglecting the "negative" wave amplitude $_(r), the non-relativistic approximation of 
the electric Poisson equation (IVI.38P is immediately written down as 

2 



d 2 ^A 2#U _ 7i {®( r 



dr 2 ^ r dr 2 s r ^ 

Thus the non-relativistic eigenvalue problem for para-positronium ultimately consists of 
both equations (1VI.42D and (1VI.44I) . 

However, notice here that the "Schrodinger eigenvalue" E* ( 1 VI . 42 1) -( j VI. 43 1) is a one- 



particle quantity and therefore cannot be identified with the (non-relativistic) RST energy 
(E$, say) of our two-particle problem. Rather, in order to get this two-particle energy 
E$ due to the solutions of the present one-particle problem ( 1VI.42I) -( IVI.44I) . one has to 
substitute these solutions into the non-relativistic descendant (E[$] , say) of its relativistic 
predecessor Em flVI.371) . Thus it becomes now necessary to deduce that non-relativistic 
functional E[$] which then simultaneously provides us also with the possibility to check 
whether the non-relativistic system of eigenvalue and Poisson equations ( 1VI.42[) -( 1VI.44I) 
does really represent the extremal equations of the wanted non-relativistic functional Er$] . 
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But of course, the desired non-relativistic version E[$] of the relativistic E[$] ( 1VI.37I) 
will simply consist of the non-relativistic constituents of Em as they are collected by 
equations (1VI.31j) - (1"V"I.36D . Therefore it is sufficient to simply list up those non-relativistic 
versions which come about by eliminating the "negative" wave amplitude $_(r) by means 
of the equation ( 1VI.41I) : 



(i) mass renormalization QVI.31]) : 

ft 2 f°° 

2 Mc 2 Zl => 2 Mc 2 - 2 (^E kin + ^E kin ) +2tt — £ v J dr $( 



(ii) radial and longitudinal kinetic energy < f VI. 32a\) - < f FX 511) 

7i H 2 f°° 7 ( d${r) \ 2 

'■(r) 



2 2M 



o 



W^=.M^-?^rdr*W ' /<I,(r) 



2 2M 
vr h 2 „ 



dr 



kin 



2 2M 



(iii) electric Poisson identity jlVI33\) - fVI35\) : 



N, 



G 



N 



M®c 2 



4 C] - M^c 2 = 
frc^ drr [p] A ■ ($( 



(iv) normalization condition $ VI. 361) : 



N* 



71 



" 2 



drr ( $(r 



1 = 0. 



d$(r) 
dr 



(VI.45) 



(VI.46a) 
(VI.46b) 

(VI.46c) 



(VI.47) 
(VI.48) 

(VI.49) 



Collecting now all these partial energies (with omission of the rest mass 2 Mc 2 ) into the 
desired functional Km as the non-relativistic limit of Em one arrives at the following 
result 



E m ti m = 2(«£ kin + ^E kin ) + 4 C] + 2A S • N* + A§ } - Ng , 



(VI.50) 
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where the Lagrangean multiplier As is related to the one-particle Schrodinger eigenvalue 
E* ( TVT431) through 

A s = -E* . (VI.51) 

Indeed this is a very satisfying result because it provides an instructive exercise to verify 
that the non-relativistic eigenvalue system ( ]VI.42l) -( 1VI.44j) actually consists of just the 
extremal equations due to the functional E[$] ( 1 VI. 501) ! Observe also that for l-p = 1 the 
present functional E[$] ( ]VI.50|) becomes identical to the functional Ep5 of ref. [T3] (see 
there the equations (V.109)-(V.113), which served to treat the first excited state 2 1 P 1 of 
positronium). 

This non-relativistic functional E™ owns now some interesting properties which are 
worth while to be considered in some detail. First, notice that both the mass renormal- 
ization factor Z\ (1VT.45j) and the longitudinal kinetic energy Mf w (lVI.46bl) contain the 
same boundary term, albeit only up to sign 

2>K^-t v r dr^r)-^^ = -n-^-e v -(^(0)) 2 . (VI.52) 
2M 1 J y ' dr 2M 1 K y " y ' 

It is true, this term cancels for the construction of the non-relatvistic functional E™ 
(jVI.50P ; but if one wishes to insist on both non-relativistic contributions fjVI.451) and 
(jVI.46b|) being individually independent of boundary terms, one has to demand that the 
wave amplitude $(r) must vanish at the origin (i. e. $(0) = 0) whenever the quantum 
number of angular momentum l-p is non-zero, cf. (1VI.29I) . Or the other way round: it is 
exclusively for tj> = that one can admit $(0) ^ 0; this circumstance must be regarded 
also for the construction of variational solutions (see below). 

Next, observe that the mass renormalization term (1VI.45I) contains minus twice the 
one-particle kinetic energy E kin (= ^'E ^ + ^E kin ), but since the one-particle kinetic 
energy ^T\^ n enters the relativistic functional EWi flVI.37p with a pre-factor of four, the 
combination of both effects yields the double one-particle energy Ey m (= ^-Ekin + ^-ffirin) 
for the non-relativistic functional Er$i (1VI.50P . as it must naturally be expected for a two- 
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particle system. This effect has been observed already in some of the preceding treatises 
[15] . Summarizing, the non-relativistic functional E[$] ( 1 VI. 501) appears physically plausible 
and logically consistent so that its practical usefulness may now be tested in a concrete 
example. 



Non-Relativistic Para- Spectrum 

The non-relativistic functional E[$] ( 1VI.50I) may now be taken as the point of departure 
for building up an approximation formalism for the non-relativistic energy level system of 
para-positronium. For this purpose, one will select some trial function $(r), obeying the 
normalization condition ( 1VI.49j) ; and furthermore one will then solve the electric Poisson 
equation ( 1VI.44I) for the potential ' p l4o(r) which is due to the selected trial function $(r). 
As a consequence, the electric Poisson identity (1VI.47P will thus be obeyed automatically 



and therefore the functional E[$] (1VI.50P becomes reduced to its physical part E^j flIV.231) 
which reads for the present situation 

= 2 (MEW + W E km ) + Eg 1 . (VI.53) 

This yields essentially an ordinary function E^(/3, u k ) of the variational parameters 
{/3, Vk} being contained in the selected trial function $(r); and thus the extremal values 
of this energy function E^ Iv ^(/3, which are to be determined by 



dp 



= (VI.54a) 
, (VI.54b) 



will then represent approximately the desired energy levels of para-positronium. 

The general structure of this level system is then obviously the following: since the 
energy function E^ Iv ^(/3, u k ) depends also on the quantum number t-p ( lVI.46cl) of angular 
momentum, one expects for any fixed value of t'p — 0, 1, 2, 3, ... a finite number of station- 
ary points flVI.54aP - flVI.54bl) of the function E^ Iv ^(/3, u k ). This latter set of excited states 
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for any fixed l-p may be enumerated by the integers n-p with n-p = ip + 1, lp> + 2,ip + 3, ... 
so that any para-level is specified by the two quantum numbers {np,£p}, see fig. 2. Nat- 
urally, for such a kind of level system there must arise the question of degeneracy, namely 
whether or not the energy of the levels due to the same rip but different ip turns out to 
be the same? Clearly, this question could be answered with sureness only if the exact 
solutions of the eigenvalue problem (1VI.42I) -( 1VI.44I) could be found which of course is not 
possible; the subsequent approximative treatment of the para-level system with only two 
variational parameters can yield only the "groundstates" (rip = ip + 1) for any i-p. 

As a simple example of the proposed approximation scheme we select a trial function 
$(r) which contains only two (real-valued) variational parameters (3 and u, i.e. we try 

<J(r) = <f> #r "e-P r . (VI.55) 

Here, the normalization constant must be fixed by the constraint of wave function 
normalization (lV"I.49j) which yields 



^ V+2 , (VI.56) 
* 7i T(2z/ + 2) ' 1 ; 

so that our trial function $(r) (1 VI. 5 51) may be recast to a dimensionless form $(y) 

$(r) = J-- . 2(3 %) (VI.57) 
y J V 7T ^(2^ + 2) V ; V ; 

with 

<%) = y» e -l . (VI.58) 
Here the dimensionless variable y is defined through 

y = 2pr , (VI.59) 

and the dimensionless form of the normalization condition f1VI.49|) appears now as 



F(2is 



— - / dy y^ +1 e-y = 1 . (VI.60) 
' + l ) Jo 
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The next step must now refer to the calculation of the first term of the functional 
flV1.53p . i.e. the kinetic energy E^in (= ^'Ek m + ^E^). The radial part hereof is defined 
quite generally by equation ( 1VI.46al) . which adopts here for the chosen trial function $(r) 
m.57ll -f lVL60l) the following form 

{r) E km = ^ (2(3) 2 ■ {r) eU") , (VL61) 
with the (dimensionless) kinetic function ^£kin(^) for the radial motion being given by 



r(2z/ + 2) J * \ dy 4(2z/ + 1) 



Analogously, the longitudinal part ^-Ekin fjVI-46c|) of the kinetic energy is found to look 
like 

[ " ] E km = ^ (2P)% ■ (VL63) 
where the kinetic function ^£ kin (z/) for the longitudinal motion turns out as 

Thus the total kinetic energy Ekm is found as 

£ kin = (r) £ ki „ + m E kiQ = ^ (2/3) 2 {^ ki » + 4 • 

(2a B /?) 2 -e k inH (VL65) 



2a B 
ft 2 

(ae = ^ 2 ... Bohr radius) 

with the kinetic function £kin( zy ) for the total motion being given simply by the sum of 
its radial and longitudinal parts, i. e. 

eU") = (r) eU") + 4 " = (l + i) • ( VL66 ) 

Let us remark here that this result ( 1VI.65I) - (1VI.66I) for the kinetic energy Eki n has been 
discussed extensively for the special cases i-p — 0, 1 in two preceding papers: putting first 
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t-p = and then v = 0, one has from ( 1VI.66I) £kin =>• 7 and thus from ( 1VI.65I) 



E kin ^^—(a B (3) 2 (VI.67) 
za B 



which is just the equation (IV. 22) of ref. p2] (hint: rescale here a-Q (3 =>■ /3 and put there 
Tn =>- for N = 0). Furthermore, putting here £p = 1 lets the present equations ( 1 VI. 6511 - 



(IVI.66P collapse to the equations (V.127)-(V.128) of ref. [13J. This provides us with the 
possibility to compare our subsequent results to those of the two precedent papers. 

The present calculations of the value of the kinetic energy functional E-^n = ^E^ + 
^i?kin upon the selected class of two-parameter trial functions $(r) (1VI.55|) represents a 
somewhat fortunate situation because the dependence of the kinetic energy on the two 
variational parameters j3 and v could be expressed by very simple elementary functions 
(e. g. .Ekin ~ P 2 , see equation ( 1 VI. 6511 ). Such a happy chance does not apply to the electro- 
static interaction energy E^ (TVI.34I) which is the second constituent of the energy func- 
tional ErlT ( 1 VI. 531) under consideration. Indeed, a comparable simplification occurs here 



only for (half-)integer values of the variational parameter v (i.e. 2v = 0, 1, 2, 3, ...);but 
for the extremalization one wishes to deal with a continuous variable u, and this suggests 
to inspect the electric Poisson identity somewhat more thoroughly, see appendix D. 
Recall here that, just on account of this identity, the interaction energy E^ ( 1VI.34I) 

*-^f*(r.^) , -^ W )f^.(^)*(V I «) 
(putting [p U (r) = 2f3a s ■ a v (y)) 

must equal its mass equivalent M^c 2 (1VI.48D . i.e. 

/OO -| poo 

dr r H4 (r) • <f>(r) 2 = -e 2 (2(3) T[2v + 2) J dyy~a u (y) ■ §{y) 2 . 

(VI.69) 

Consequently, the electric Poisson identity reads in dimensionless form 
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which may be verified also directly through integration by parts with use of the dimen- 
sionless form of the Poisson equation ( 1VI.44I) 

d 2 a u (y) 2 da v {y) _ _ y y 2 ^ 1 (VI 71) 

dy 2 y dy T(2u + 2) ' 1 ' 

The problem is now to determine the potential function e vo t{v) by one or the other of 
the two possibilities (1VI.70P . not only for (half-) integer values of the variational parameter 
v (as it has been done in the preceding paper [13]) but by treating v as a truly continuous 
parameter. Clearly it is only for this case that the extremalizing condition (1VI.54bl) can 
be meaningful, with the energy function E^ Iv ^(/3, u) being given by 

^ w \(3, v) = 2E kin + 4 C] = — {(2a B (3) 2 ■ £ ki » - (2a B 0) ■ e pot (u)} (VI.72) 

where the energy scale is determined by the well-known atomic unit (a.u.) 

e 2 

1 a .u. = = 13.605... [eV] . (VI.73) 

The specific form of the present energy function E( IV )(/?,z/) NL72§ has been deduced 
for the first time in the preceding papers (see equation (V.138) of ref. |13j); and it has 
already been pointed out that the well-known virial theorem [10J is implied by that type 
of energy function (I VI. 720 . Indeed, this can easily be verified by simply evoking the first 
equilibrium condition ( lVI.54al) which yields for the present situation (1 VI. 721) 



. = Ae 2 | (2a B P) ■ e kin (v) ~ ^e pot (u) | = . (VI.74) 

This says that the energy-minimizing value of (3 is given (for arbitrary v) by 

2a B /5, = |^T ( VL75 ) 
2e k in(^) 

which then yields for the kinetic energy flVI.651) 

E km (P*, v) =► ^kin(^) = ^— ■ ^4^r (VI.76) 

2a B 4e kin (^) 
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and analogously for the interaction energy ( 1VI.68I) or ( 1VI.69I) . resp., 



4 ei = pc 2 =► E pot (^) = • ^44 • ( vL7r ) 

Consequently, the kinetic energy 2£kin(z / ) (1VI.76I) at the equilibrium value of (3 equals 
always (minus) half of the interaction energy E^, i.e. 

2 • E kin (v) = -1 4 C] = -1 M^c 2 , (VI.78) 

which in the conventional theory holds for the Coulomb-like interactions and is generally 
known as the "virial theorem" [TOl [T7j. But the amazing point is here that this virial 
theorem (IVI.78P does hold also in RST quite generally for the electrostatic approximation 



where the interaction potential ^Ao(r) 

[p %(r) = ^ f — ,,51,, ^ * , (VI.79) 



jyW | | ry ry ' | | ry 

as the solution of the Poisson equation (1VI.44I) . adopts the Coulomb form only in the 
asymptotic region (r — > oo), but not for finite distance r! Notice also that the RST virial 
theorem holds for arbitrary values of the variational parameter u, not only for that value 
(z/*, say) which further minimalizes the energy function E< IV ) (/3, v) flVL72|) beyond the 
minimalization through the first variable f3. 

This final step of energy minimalization is described now by the second equilibrium 
condition (1VI.54bj) which then reads for our present energy function E^ IV ^(/?, v) (1 VI. 721) 

= £. { (20B p ,y . «W!0 _ (2ob;3>) . )± o (vi.80) 

and thus would evidently require also the explicit knowledge of the potential function 
£ P ot( z/ ) flVI-70p . analogously to the kinetic function e^^u) (1VI.66|) . However, the deter- 



mination of the analytic form of the potential function e pot {v) is a difficult problem so 
that we have to be satisfied here with a power series representation (see appendix D). 

If the interaction energy of both positronium constituents is interpreted as the energy 
E^ (IVI.68P being concentrated in the gauge field ^A (r), then the potential function 
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£ P ot(^) is found as 



-pot 



V) 



2 f da v (y)\ 2 1 1 T(Av + 3 + m + n) 

dyy • -*H =^t5 



% 7 2 4 ^+ 3 ^ 2 m + n r(2z/ + 3 + m) -r(2z/ + 3 + n) ' 

' m,n=0 ' y 

(VI.81) 

cf. equation flD.8p . On the other hand, if one prefers to interprete the physical interaction 
energy as the mass equivalent M^c 2 ( 1VI.48I) then the potential function £pot(^) is found 
to be of the form 

1 f°° 

£ P otO) = - ^ J dyyh u (y)-<$>(y) 2 



2u 



1^ 1 y» r(4z, + 2 + n) 1 

+ l\ 2 4 -+ 2 ^2" T(2z/ + 2) -r(2z/ + 2 + n) J ' 1 J 



cf. equation (ID.180 . Consequently, the electric Poisson identity implies the numerical 
equality of both potential functions e pot {u) (jVI.8ip and f1VI.82h which may be checked 
explicitly by selecting some integers for the variational parameter 2u, see the table in 
appendix D and fig. 3. 

The desired minimal value (Ei IV \ say) of the energy function E^ Iv ^(/3, v) (j VI. 72ft is 
found by resolving both equilibrium conditions ( 1VI.74D and ( 1VI.80[) for the equilibrium 
values /3* and v* and substituting this into the energy function such that 

EM=EW(J3,u) ^ = — {(2a B f3*) 2 -e kin (v*)-(2a B {3*)-e pot (is*)} • (VI.83) 

However, in contrast to the first equilibrium condition (1 VI. 74ft . the second equilibrium 
condition (lVI.80ft is hard (if not impossible) to solve exactly for the equilibrium values 
v*\ and therefore we have to resort to an approximative solution of the extremalization 
problem defined by those equations ( 1 VI. 741) and ( 1VI.80I) . To this end, one substitutes both 
the kinetic energy E kin (f3*, v) (lVI.76j) and the potential energy = E pot (z/) fl VI. 77ft 
into the energy function WP^(j3,j/) flV"I.72[) which yields a resulting function Ex(V) (== 
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E (IV) (/^)) 



Et(i/) = 2£ kin (z/) +E pot (z/) 



e 2 f 4ot(^) eJotWl e 2 £ 2 ot (^ 



(VI.84) 



2a B t2£ki n (i/) e k h» J 4a B £kin(^j 

Recalling here the fact that the conventional non-relativistic groundstate energy (E conv ) 
of positronium is given by 

E coav = = -6.802 . . . [eV] (VI.85) 

4a B 



the present RST energy function E T (z/) (1VI.84j) does appear as 



Er(u) = E conv ■ S v (u) (VI.86) 
with the spectral function S-p being defined by 

S v {p) = -&)-L . (VI.87) 

£kin(^) 

Now according to the principle of minimal energy, the extremal points of the spectral 
function S-p{v) determine an approximate solution of the original eigenvalue problem 
defined by the former equations flVI.42p - flVI.44p where the extremal values /3* and v* 
do then fix the corresponding wave function $(r) flVI.55D . potential H4 (r) II VI. 79ft and 



associated energc E* ( 1VI.83H . But on behalf of the simplicity of the chosen trial function 
$(r) (1VI.55P one will not expect to obtain by this procedure the totality of para-states; 
but rather one will find the state of lowest energy for any admitted value of the quantum 
number l-p of angular momentum, see fig. 2. By the former convention, this state has 
principal quantum number n-p — i-p + 1. By comparison with the conventional theory 
(see ref. [12]), the conventional spectrum 

p 2 1 

E(n) = — _ _ (VI 88) 

conv 4a B n 2 { } 

(n = 1,2,3,4,...) 
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would obviously be reproduced by the present RST approximation method if the minimal 
values of the spectral function Sp ( 1 VI. 8 71) turned out for any prescribed t-p as 

M =^=w^- (VL89) 

mm fl.p \tp ■+* 1 ) 

In a preceding paper [T5], the conventional groundstate energy Ecorw'* ( 1VI.88I) has 
exactly been reproduced for lp = by choosing v = and extremalizing the energy 
function merely with respect to the remaining variational parameter (3. But clearly, since 
we are dealing here with an approximation method, the true RST groundstate prediction 
(in the electrostatic and spherically symmetric approximation) must be somewhat lower 
than the conventional value (jVL85]), see the discussion of this in ref. |12j . In a subsequent 
paper [13], the first excited positronium state (£p = 1) has been treated by means of 
the present approximation method with two variational parameters (3 and v\ and the 
corresponding RST prediction did approach the conventional value Ecoi^ ( 1VI.88I) up 
to some few percents. However, this result came about by approximating the spectral 
function Sp{y) flVI.871) by its values upon the (half-) integers v — 0, |, 1, |, ... together with 
the help of an interpolating polynomial. But clearly, the reliability of the RST predictions 
can be improved now by numerically calculating the spectral function Spiv) (1VI.87P for 
continuous parameter v and thus looking for its minimal values for any ip (= 0, 1, 2, 3, . . .), 
see fig. 4. 

Perhaps the most remarkable result of the present calculations refers to the positron- 
ium groundstate (i.e. n-p = <^ i-p = 1). Here the conventional value (IVI.85P is exactly 
reproduced also with the present two-parameter trial function $(r) f)VI.55[) . This is a 



somewhat amazing result because the same prediction could be obtained in a preceding 
paper [15] with only a one-parameter trial function which emerges from the present two- 
parameter trial function (1VI.55P by putting one of the two variational parameters to zero 
(i.e. v — > 0). Indeed, both coinciding predictions are consistent because, in the present 
two-parameter approach, the groundstate energy emerges as a boundary minimum, namely 
just for v = (see fig. 4). 
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Concerning the energy of the excited states (n-p = 2, 3, 4, 5, 6) the RST predictions are 
found to differ from the conventional values E<£l y ( 1VI.88I) by roughly 10% (see the table 
below). This difference may be attributed partly to the use of the spherically symmetric 
approximation and partly to the use of an imperfect trial function $(r) (1VI.55P with 
only two variational parameters f3 and v. A better trial function <£>(r) would surely shift 
the RST predictions towards its conventional counterparts -Eionv since the RST energy 
functional Er$i flVI.50p is bounded from below for any chosen value l-p of orbital angular 
momentum (strictly speaking, the quantum number l-p emerges only in the context of 
the product ansatz ( 1VI.19al) -( 1VI.20bl) ; but it is assumed here that the hypothetical exact 
solutions of the RST eigenvalue problem can also be enumerated by some quantum number 
of orbital angular momentum which then is expected to have the same classifying effect 
upon the set of solutions as is the case with the present i-p). On the other hand, for the 
purpose of obtaining more exact solutions to the RST eigenvalue problem one should go 
beyond the spherically symmetric approximation of the interaction potential ^l4 (r, ■#), cf. 
( 1VI.18j) ; and whether or not this also shifts the RST predictions towards their conventional 
counterparts appears to be unclear at the moment. 
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n v (= i v + 1) 


E£l [eV], fjVI.88) 


E T (^) [eV], (|VI.84|) 




TC 1 <-,, \ Z7'( n ) 
^TXy*) — -C/conv ro/1 

r \ \M 






\E (n) 

\ - 1 — ' conv 


1 


-6,802. . . 


-6,802. . . 








2 


-1,700. . . 


-1,550... 


1, 794 .. . 


8,8 


3 


-0, 755 .. . 


-0, 669 .. . 


3, 752 .. . 


11,4 


4 


-0,425. . . 


-0, 373 .. . 


5, 874 .. . 


12,2 

* 


5 


-0,272. . . 


-0, 239 .. . 


8, 130. . . 


12,1 


6 


-0,188... 


-0, 166. . . 


10,504. . . 


11,7 


10 


-0, 06802 . . . 


-0,06068... 


20,953... 


10,7 


20 


-0,01700... 


-0,01554... 


51,919... 


8,5 


50 


-0, 00272 . . . 


-0,00256... 


172,161... 


5,7 



Observe here that the relative deviation of the conventional and RST predictions (last 
column) decreases for increasing principal quantum number n-p (e.g. 3,7% for n-p = 100)! 

2. Ortho-Positronium 

The discussion of ortho-positronium will be carried through along a certain path of 
arguments very similar to those which were presented for para-positronium in the pre- 
ceding subsection. Especially, we will be satisfied again with the electrostatic, spherically 
symmetric, and non-relativistic approximation. The consideration of the magnetic in- 
teractions will be deferred to the next subsection when both the para- and ortho-level 
system will be sufficiently worked out in the above mentioned approximation which then 
will serve as a first rough estimate of the energy difference of the ortho- and para-versions 
of the positronium states. 
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Eigenvalue Problem for the Ortho-States 



The point of departure is of course the ortho-system flV.55ap - flV.55dj) for which the 



boson number 5 must be put to minus the value of the para-configuration, i. e. 5 =>- — |. 
Furthermore, the electrostatic approximation requires to neglect the magnetic potential 
^A^. With these two specializing arguments the ortho-system (lV.55a[) - (1V.55d[) becomes 
simplified to the following system: 

i d®H A 



M + M 'c<^_ (VLDOa) 




h 



dr 



i_ 9% 
r <9</> 



- • % 



2i</> 



1 ggjgj. cot?? l (b , 
r dtf r L 2 H 



M + M, 

h 



c ■ (6) <S_ (VI. 90b) 



TT— + — • ( ^ + ( ^0 • ^ 

ar r a© 2r 



-2i<)> 



cot-d r 

+ — [ 



2 



1 f)sh i 



M-M* 

h 



c • {b) TZ + (VI.90c) 



+ + tt • 5 - + u ° • 5 + 

ar r a<p 2r 

+ e« .{i.^^-^^}.^..^ . ( VI,0 d) 
This eigenvalue system is to be complemented again by the electric Poisson equation 

A( fe U = -a. ■ • + ■ + «jL • ^ + ■ . (VIM) 
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The crucial point with the present ortho-system ( 1VI.90a[) -( lVI.90d[) refers of course to 
the emergence of its triplet character, i. e. one wishes to see in what way there do arise 
from the present ortho-system three different types of solutions (with j z = 1, 0, —1) which 
have in common the same energy albeit only for the special case of neglection of the 
magnetic interactions. Furthermore with the latter presumptions, the common energy of 
these ortho-states should be the same as the energy due to the corresponding para-states, 
see the discussion of the analogous situation in the conventional theory (below equation 
(jYH). 

Vanishing Angular Momentum j z = 

Naturally, it is very tempting (but not correct) to think that the RST counterparts 
of the conventional triplet states with S — 1, S z — are included in the set of solutions 
of the ortho-eigenvalue system ( lVI.90al) -( lVI.90dl) . namely in such a way that both kinds 
of wave amplitudes ^TZ± and ^S± do occur in a symmetric way, just in order to have a 
vanishing ^-component (j z = 0) of total angular momentum J z . In this sense, one tends 
to try the following ansatz for the complex- valued amplitudes ( - b ^JZ±(r) and ( b 'S±(r): 

(b) R ± {r, 0, 0) = e"^ • {b) R±{r, <&) (VI.92a) 
{b) S±{r, 0, <P) = ■ {b) S±{r, 0) , (VI.92b) 

where the reduced wave amplitudes ^i2±(r, "&) and ( b 'S±(r, i?) are adopted to be real- 
valued and independent of the azimuthal angle <fi. Inserting this ansatz into the present 
ortho-system flVI.90aP - flVI.90dl) yields the following system for the real- valued amplitudes 
R±i S±, which for the sake of convenience may be transformed to ^R±, ^S± accord- 
ing to 

®R ± = Vrsintf • {b) R ± (VI. 93a) 

^S± = Vrsintf • ^S± . (VI.93b) 
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This lets reappear the eigenvalue system ( 1VI.90al) -( 1VI.90dl) in the following form: 

£^± + I . _ ( % . m n _ = c . W 5. (v ,. 9 4a) 

Or r ov n 



or r ov n 



d ^L + I . « ft. _ I . £^|z + < % . ® A, = " c . Wfi+ (VI.94C) 
Or r r ov a 



^ + I . WS _ + I . ^ + ( % . M S+ = ^ c . WS+ . (V I.94d) 
ar r r ov a 



However, this special subcase of the ortho-system (1VI.90aP - (lVI.90dj) is exactly identical 



to the para-system (1VI.lal) -( lVI.ldl) ; and therefore the present ortho-system (lVI.94aj) - 
(IVI.94dj) does actually not describe a new quantum state of the ortho-type! 



The reason for this "failure" of the ortho-ansatz (jVI.92a|) -f lVI.92b|) is the following: 
If that ansatz is resubstituted into the Pauli spinors ^f±{r) ( 1V.50al) -( 1V.50bl) with the 



basis spinors being given by ( lV.48al) -( 1V.48dj) for boson number 6 = — |, then one 
will ultimately end up with the following decomposition of the proposed Pauli spinors for 
the ortho-system 

<V = {b) R + ■ 4 +) + {b) S+ ■ J - ] (VI.95a) 



{®K..ui +) + ®S..-ui- ) } (VI.95b) 



where now the basis spinors cjgi are those for 8 = +|, i.e. one recovers merely the 
para-case being described above the para-system ( 1V.56al) -( lV.56dl) . In other words, the 
ortho-ansatz 0VI.92aft -( IVL92bl for the wave amplitudes rescinds the transition from the 
para-basis (6 = ~) to the ortho-basis (8 = — ~). It is very instructive to reformulate 
this effect also in terms of angular momentum: The action of J z on the ortho-spinors 
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(lV.50ap - dV30bll reads explicitly 

JW <V = (L^>JZ + ) ■ w « + (L.C*lS + ) ■ W H + (^ + . JM U W + . 

(VI.96a) 

-i { (4 (6) ^-) ■ ^S +) + {L z [b) S-) ■ + 0>>R_ ■ J<M +) + {b) S- ■ JtM' } ] , 

(VI.96b) 



where the ortho-form (S = — |) of the basis spinors cJq 1 flV.48ap - flV.48dp is to be used 
But here one concludes for the ortho-value (i. e. 8 = 
eigenvalue equations of angular momentum ( TV.49al) -f rV.49bl) 



|) of the boson number 6 from the 



J i+) J +) 



h ■ uj, 



(+) 



(VI.97a) 
(VI.97b) 
(VI.97c) 
(VI.97d) 



and similarly for the wave amplitudes {b " r R ±1 ^S± f lVI.92al) -f rVL92bl) 

L z {b) Tl ± = -h- {by R± (VI.98a) 
L z {b) S ± = h ■ (fe) <S± . (VI.98b) 

Obviously, both the spinor basis and the wave amplitudes (b >K ± , {b) S± carry the 
expexted quantum of angular momentum, but if this result is inserted into the relations 
for the Pauli spinors flVI.96ap - flVI.96bp one recovers the para- situation 

j(+) ( b ) ip+ = o (VI.99a) 



o 



(VI.99b) 



Thus the result is that the orbital angular momentum of the wave amplitudes compensates 
the (total) angular momentum of the spinor basis so that the Dirac spinor field i/jq 



(VI.100) 
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would carry here vanishing angular momentum: 

J z ^ = j(+) © j(-) ( V = . (VI.101) 

In contrast to this, for the para-state ipp = © ^V- neither the spinor basis nor the 
wave amplitudes do carry any angular momentum; but the net effect (j z = 0) is the same. 

Summarizing, the triplet states with j z = cannot be described by those field equa- 
tions where the associated current j(r) (or k(r), resp.) encircles the z-axis, as shown by 
equation (IV. 47b j) . Rather, one must evoke here the toroidal geometry which has already 



been studied in a previous paper [IT] . This, however, requires an extra discussion not 
to be presented in this paper. For the subsequent discussion of the ortho-states we will 
restrict ourselves to the RST counterparts of those conventional states with j z = ±1, cf. 
flV.lOajl - flV.lObjl . 



Ortho-States with j z = — 1 

The states with definite orientation of their total angular momentum J z along the 
(positive or negative) z-axis can intuitively be associated with an ortho-current kb which 
encircles the z-axis in the usual way 

^ = ( %e^, (VI. 102) 

cf. the analogous para-case fjV.6ip . For analyzing this ortho-situation, one has to restart 
from the original eigenvalue system flVI.90ap - flVI.90dl) for the ortho-states, but now avoid- 
ing such an ansatz for the wave amplitudes like (lVI.92aP -f lVI.92bl) which turned out to 
be adequate for the para-states. Clearly, the wanted ansatz for those ortho-states with 
definite orientation along the z-axis must respect the fact that either the positive (z > 0) 
or the negative (z < 0) direction is singled out; and therefore both wave amplitudes 
( b '7Z± and w«S± cannot enter the description of a triplet state in such a symmetric way 
as it was the case for the sing let states, cf. fjV.54ap - fjV.54bp . Thus we try now for the 
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treatment of ortho-positronium the following ansatz for the wave amplitudes {by Jl±, ®S± 
flV.50aj> - flV.50bj> : 



-2i(/> 

(fc) K ± (r, <&, 0) = = • {b) R±(r, <&) (VI. 103a) 

sin^vr sin-i? 

(6) <S± (r, <&,</>) = • W S± (r, 0) , (VI. 103b) 

where the new amplitudes ( b ' ) R± and ^S± do not depend upon the azimuthal angle (p. 
Obviously, the present ortho- ansatz (1VI. 103al) - (1 VI. 103bl) is the analogue of the para-case 
flVOa] - flVObl) . 

Following now further the line of those arguments which were put forward for the para- 
configurations in order to reformulate the eigenvalue equations in terms of the new am- 
plitudes WR±, {p) S± (cf. (TVI.lall - flVI.ldjl ). one substitutes the present ansatz flVI.103ajl - 



(jVI.103bj) into the original ortho-system of eigenvalue equations (jVI.90ap - (1VI.90dl) which 
then reappears in the following form: 

^± + I . <«fl + + «4 . = ™ c • <»>/L (VI.104a) 

Or r r ov n 



§M±-L. m§+ L_ . _ < % . w §_ = « + c . (V i.io4b) 

or r rsin v ov n 



9^ _ sin^ ^ + = M-M, c 

Or r ov n 



^ + 1 . ms_ + i . ^ + < % . « S+ = ^ c • ™S + . fVLlOM) 
ar r rsnri? ai? ft 

Indeed, this is just the ortho-counterpart of the former para-system AVI. laj) - (jVI.ldh : and 



therefore it suggests itself to look here also for some product ansatz similar to the para- 
case ( jVI.19aj) -( jVI.19b|) ; i.e. we try here again the spherically symmetric approximation 
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(in its ortho-form) 

(6) .R ± (r,<i?) = "R±{r) ■ g R {$) (VI. 105a) 

®S±(r, &) = "S±(r) ■ g s (#) . (VI.105b) 

For the angular functions and gs{"&) we postulate the coupled first-order equations 

S mH-^^- = i -g R ^) (VI.106a) 

1 d9Rm -io-9s(#) (VI.106b) 



sin 2 

(^c^o = const.) , 

since by this choice the original ortho-system (lVI.104ap - (lVI.104dl) becomes neatly reduced 
to a system of ordinary differential equations 



d"R + (r) 
dr 


+ - ■ "R + 
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(r) + 
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"S+[ 


r) - [b] A ( 
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"R. 


(r) 


M + M* 
= ft C " 




(r) 


(VI. 107a) 


d"S + {r) 
dr 


- - ■ "S + 
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"R + 


(r) - 1% 


(r) 




(r) 


M + M* 

= ft c " 




(r) 


(VI. 107b) 




d"R_ 
dr 
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io 
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r) + [b U ( 


r) • 
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M — M* 
= ft C " 
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(VI. 107c) 


d"S-(r) 
dr 


+ - • 
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(r) + 


£° 
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"R. 


(r) + [6 k 


(r) 


"5+ 


(r) 


M-M* 
= ft C " 


"5+ 


(r) 





(VI.107d) 

Evidently, this system is again the ortho-counterpart of the para-case (lVI.21a[) - (1VI.21dh 
where the role of the para-number E-p is now played by the two ortho- numbers lo and Iq- 
In order to clarify their meaning, we defer for a moment the discussion of the present ortho- 
system flV1.107ap - flVI.107dD and first analyze the angular system flV1.106ap - flVI.106bD 
which, of course, is the ortho-counterpart of the former para-system (IVI.20al) -( rVL20bl 
but is considerably more complicated. 

In order to see more clearly the type of solutions admitted by that angular system 
( jVI.106a)) - (lVI.106b|) one differentiates once more and thus deduces a decoupled second- 
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order system (as the ortho-analogue of the para-case ( lVI.24aft -l jVI.24bj) ): 

^01 + 2coti) . d jgl = {e J o) . gsW ( vi.l08a) 

*2@. - 2cottf • d -igl = (W' o) • ShW . CVLlOBb) 

Sometimes it may be rather instructive to first consider an unphysical solution of the 
problem since this lets one truly appreciate the elaboration of the physically correct case. 
Such an unphysical solution can easily be guessed and is given by 

9s{#) (VI.109a) 
sin 17 

g R (>&) coad (VI. 109b) 

where the constants to, to must both be fixed by 

t Q -1 (Vl.llOa) 

to^-1. (VI. 110b) 



However, this solution must be rejected because the wave amplitudes (jVI.103ap - (jVI.103bl) 
would adopt the following form 

p -2i</> „„„ „Q . p —2i<j> 

VR+Mfi = . -== .g R {<t>) ■ "R ± (r) =► . , — ■ — — • "Mr) (VLllla) 
sin v v f sm 1/ sin 17 V ^ sin 17 

^ ± (r,i9,0) = J**±.g 8 {p) ■ "S ± (r) => -J= • "S ± (r) . (Vl.lllb) 

V T yjr Sin 17 

Indeed, since the wave amplitude ® r R.± (jVI.lllaj) diverges like ~ (sini?)~2 on the z-axis, 
the normalization condition ( 1III.49I) cannot be satisfied! And clearly it makes no sense 
to deal with non-normalizable wave functions. The normalizable solutions of the angular 
system flVl.106ajl - flVI.106b!) . or flV1.108al) - (IVI.108bD . resp., are worked out in detail in 



appendix E. 

For a better understanding of the ortho-system it is very instructive to elucidate the 
meaning of those constants £@, to emerging in the angular systems (jVI.106ap - (jVI.106bl) 
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or f lVI.108a[ )-f lVI.108bl) . resp. To this end, one passes over to the non-relativistic limit 
of the one- dimensional eigenvalue system ( lVI.107al) -( 1VI.107dl) . Here the usual procedure 
consists in first identifying both "positive" and both "negative" amplitudes with each 
other, i. e. one puts 

"R±(r)= "S±{r)=n±{r), (VI.112) 

so that the radial eigenvalue system becomes reduced to only two equations for Q±(r). 
However, this reduction can be performed in two ways. Here the first version is the 
following: 



dQ + {r) 
dr 



+ — -Sl + (r)-">%(r)-SUr) 



^LM_^.n_ (r) + ^ o(r) .n +(r) 

dr r 



M + M* 

h 

M-M* 



c-Q(r) 



h 



(VI.113a) 
(VI. 113b) 



which is valid when the introduced ortho-number £q of orbital angular momentum is 
linked to the angular constants £q and £q by 



9 o = ^ ~ 1 

!o = -(£o + l). 



(VI. 114a) 
(VI. 114b) 



The second version looks slightly different 

±L± - . n + ( r ) - [b] A (r) ■ fi_(r) = 



dr 



dfljr) 



+ 



to 



+ 1 



■fi_(r) + [ fe U (r)-fi + (r) 



h 

M-M, 



c-Qjr) 



-c-ttjr) 



dr r ' h 

where the link of the ortho- number Iq and the angular constants £q reads now 



(VI. 11 5a) 
(VL115b) 



?o = -(^o + l) 



(VI. 116a) 
(VI. 116b) 



Observe here that both links flV1.114ap - flVI.114bp and fl V1.116aj) - flVI.l 16bh entail a com- 
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mon form of link, namely 



£ .£ = l-£ 2 (VI. 117a) 

£o+io = -2, (VI. 117b) 

and this is closely related to the posible values of £q 2, 4, 6, 8, ...), see appendix E. 

Clearly, the present ortho-number £q is the ortho-counterpart of the para-number l-p 
defined through (1VI.20al) - (rVI.20bl) ; and this ortho/para analogy refers also to the pair of 
radial functions $±(r) (1VI.22|) and f2±(r) flVI.112p . However, both corresponding eigen- 



value systems, i. e. the former para-case m.23al) -f lVL23bl) and the present ortho -cases 
(IVI. 1 13a]) - (IVI. 1 13bl) plus ( IVI. 115a]) - (IVI. 115bl) . own a rather different degree of complex- 
ity so that the treatment of the ortho-case requires an extra discussion (see appendix E) . 
On the other hand, the relative intricacy of the ortho-case disappears with respect to 
the eigenvalue equations if one writes down their non-relativistic approximation which is 
shared by both versions HVI. 1 13ap ~ (IVI. 1 13bj) and (IVI.115aD - (IVI.115bD : 



Clearly, since the ortho-number £q appears here in quadratic form ( ~ Pq) both relativistic 
versions of the eigenvalue equations must have the same non-relativistic limit ( IVI. 1181) ; 
namely observe that the change of sign Iq — > — io leaves invariant the non-relativistic 
equation (IVI. 1181) but interchanges both relativistic versions (IVI. 1 13aj) / (lVI.113blW 



(IVI. 115a]) / (jVI.115b|) . Furthermore, the interpretation of the ortho- number Iq as a quan- 



tum number of orbital angular momentum seems legitimated by its appearance in con- 
nection with the centrifugal term ( ~ ^o) in the Schrodinger-like equation (IVI. 1 18l) . cf. 
the para-counterpart (IVI.42[) . 

Finally, it is also instructive to convince oneself of the fact that our improved ansatz 
( IVI.lllaD - flVI.lllbl) (with general g R {-d) and g s ($)) for the wave amplitudes (by JZ ± , {b) S± 
actually equips now the Dirac spinor ipo (IVI. 100|) with the desired one quantum of angular 

117 



momentum, i. e. 

Mo = -Ho , (VI. 119) 

in contrast to the spoiled result (jVI.lOip ! The present nice result for the ortho-wave 



function ipo comes about by observing its Whitney sum structure (cf. (jVl.lOOj) ) with 

= -h-Mp+ir) . (VI.120) 

Here it is especially interesting to observe that the general wave amplitudes ^R,± (IVLlllaft 
carry two quanta of orbital angular momentum, i. e. 

L z {b) K ± = -2h {by R± , (VI.121) 



which however is partially compensated by the one quantum being carried by the spinor 
basis Uq 

J«4 ±} = ±fr-4 ±} • (VI. 122) 

On the other hand, the amplitude field ( b ]S± (IVI.lllbl) carries no angular momentum at 
all 

L z ib) S± = 0, (VI. 123) 

whereas the spinor basis contributes one quantum 

J^u^ = ±h-u{ ±) , (VI. 124) 

cf. ( 1V.49al) -( 1V.49dl) . Thus the negative Pauli spinor carries also one quantum of total 
angular momentum 

jH < b V-(r) = -/i (6 V-(r) , (VI. 125) 

so that both results flVI.1201) and (IVI.125I) actually add up to the claimed result flVI.HQh 
for the Whitney sum ip ( 1VI.119I) . 
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Energy Functional for the Ortho-States 

The ortho-eigenvalue system (either in its relativistic form ( 1VI.115aT) -( lVI.115b|) or 
in its non- relativistic approximation ( 1VI.118I) ) appears to be very similar to its para- 
counterpart (lVI.23ap - (1VI.23bj) . or (1VI.42I) . resp.; but this should not lead one astray to 
think that such a formal similarity is sufficient in order to ensure the numerical equality 
of the corresponding energy spectra (due to the ortho- and para-types). In fact, for such 
a numerical equality (i.e. the ortho /para degeneracy) it is necessary that two additional 
conditions must be satisfied; namely the numerical equality must also refer to (i) the 
interaction potentials (i.e. ^A H4 ) and (ii) to the energy functionals (E[$j <^ Ep). 
This means that, for the sake of comparison, we have to set up now the ortho-functionals 
E[n], E[n] as the ortho-counterparts of the para- functionals (1VI.37[) and E[$] (1VI.50P . 
resp. 

Turning here first to the relativistic form (-Bp), one restarts with the proposed ortho- 
ansatz flVI.lllal )- flVI.lllbj) . with general angular functions gn($), gs($) an d under ob- 



servation of the ortho-identification (IVI.112P ; and one substitutes this ansatz now into 



the original energy functional E\t] ( 1IV.20I) with vanishing magnetic and exchange objects. 
Thus one obtains the desired ortho-counterpart E[a\ of the para-functional E^ ( 1VI.37I) 
in the following form 

E M = 2Mc 2 ■ Zl + 4Wf kin + 4 el + 2A D • N n + A§> ■ ivg 1 . (VI.126) 

The problem is now to adapt all the individual contributions to the peculiarities of the 
ortho-configurations. First, consider the mass renormalization factor Zq 



Zl = J d 3 rMr) <Mf) = J d 3 r | (6 Vf +) (r) (6 V( +) (r) - (6 VJ_)(r) (6 V(-)(r) 

(VI. 127) 

which reads in terms of the wave amplitudes ( 6 fc±, ^S± (lV.50al) -( IV30bl) 

Zl = j |(^ + . (% + _ . (6)^_ + (6)£ + . (b) s+ _ . ( b)(S _| (VL12g) 
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But substituting here our general ortho-ansatz ( 1VI.lllaT) -( lVI.llTbl) lets reappear the mass 
renormalization factor Zq in the following factorized form 

(dtt = smddddcf)) . (VI. 129) 

This should be compared to the wave function normalization and thus suggests to establish 
separate normalization conditions for the angular and for the radial parts of the wave 
amplitude, cf. (1111.491) and flTV46l) 



1 1 / d 3 r ( %(r) = / d 3 f | (b V+(f) {b) Mr) + ( V-(0 {b) ^-if) 



^hlW +sM . glw ).m±m, (v ,,30) 



4tt \sin 3 $ b J r 

i. e. we satisfy this normalization condition by putting separately 



f 9 Q 2 Jr)+Q 2 (r) 
Nq = drr 2 +V ; - K J -1 = 0. (VI.131b) 



With this convention, the mass renormalization factor Zq flVI.1291) becomes reduced to 



Z 2 = drr 2 =^ . (VI. 132) 



Clearly, this is again the ortho-counterpart of the para-case ( 1 VI. 311) . 



Concerning now the kinetic energy ^T^n, i. e. the second contribution to the energy 
functional Em\ (1VI. 126j) . one is forced to plunge into some technical details which are 



presented in appendix F. The result is, cf. equation (1F.15I) 
MT kin = -- [drr 2 p+^-W fi-(r)rffi + (r) | (1t2V n+(r) • fL(r) 



2 J \ r dr r dr r 2 

(VI. 133) 
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Furthermore, the electrostatic interaction energy in the spherically symmetric ap- 
proximation is the same as for the para-configurations, cf. (IVI.34() 



,l He f°° , / d^A (r) \ 
4 ] = -- / dr (r-—^) , (VI.134) 



a s J V dr 



2 



whereas the ortho-counterpart of the electrostatic Poisson identity ( 1VI.35I) requires to 
restart from the general form of the mass equivalent M^c 2 

M (e) c 2 = -He J d 3 f {b) k (r) ■ {b) A (r) . (VI.135) 



But by observation of the ortho-density ^k (r) (jVI.130p this object factorizes again to 



M N c 2 = -He 



dVL 
Air 



+ sin^ • gt;($) 



n 2 Jr) + Q 2 (r) 



drr 2 ( 6 U (r) ■ +V ; ~ V ; , (VI.136) 



sin 3 "d 

i. e. by reference to the angular normalization ( !VI.131aT) 

M^c 2 = -He J drr 2 ®A (r) ■ ^±±^1 . (VI.137) 

Thus the electric Poisson identity reads also in the ortho-configuration in the spherically 
symmetric approximation 

N§ = E§ - M [c] c 2 = , (VI. 138) 

where E^ is given by fl VI. 134ft and M^c 2 by fIVI. 137ft . The important thing here is that, 
apart from the different normalizations of the wave functions (cf. the para-normalization 
(IVI.36P vs. the present ortho-normalization flVI.131bj) ). both mass equivalents flVI.35ft 
and ( 1VI.137I) are formally identical as expected; and this similarity holds then also for the 
corresponding Poisson identities, cf. (1VI.33|) vs. (1VI.138j) . However, it must be stressed 
that these formal similarities between the ortho- and para-configurations are valid only for 
the spherically symmetric approximation because here the different angular dependencies 
become hidden behind the separate normalization to unity of just the angular parts of 
the wave amplitudes! 
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With all the individual contributions to the ortho-energy E\qi ( 1VI.126I) being now 
specified, it is again very satisfying to find the former relativistic ortho-system ( 1VI.113aT) - 



( 1VI.116bl) being identical to just the extremal equations of the present ortho-functional 
(i.e. 8E\q] = principle of minimal energy). Especially, the solutions of the electric 
Poisson equation for the ortho-configurations 

^U (r) 2 S b %(r) Ql(r) + Ql(r) 

1 1 = ~ a B (V1.139) 

ar z r ar r 

extremalize the energy functional Em] (IVI. 126|) with respect to ' 6 l4o(r) within the class of 
spherically symmetric trial functions ' b i4o(r). Therefore this equation ( IVI. 139ft is then to 



be considered as the spherically symmetric approximation to the original Poisson equa- 
tion (1VI.91P for the ortho-configurations, i. e. the ortho-counterpart of the para-equation 



( jVOSl) . 

Of course, the similarity of the present ortho-equation (IVI. 1391) with its para-counterpart 
( 1VI.38P raises now the general question of the extent to which the former para-eigenvalue 
problem and the present ortho-eigenvalue problem are representing essentially the same 
mathematical problem? Indeed, this ortho/para equivalence (within the framework of 
the spherically symmetric approximation) holds not only for the above-mentioned elec- 
tric Poisson equations but equally well for both eigenvalue systems (cf. the para-case 
( 1VI.23ap - flVI.23bp to its ortho-analogy f1VI.115ah - (1VI.115bp ): and it holds also for both 
energy functionals E^j (1VI.37h and Em] ( IVI. 126ft . Obviously, both eigenvalue problems 
are completely equivalent from the mathematical viewpoint which in physical terms is 
nothing else than the expected ortho/para degeneracy. And therefore it must be clear 
that both eigenvalue problems do generate, albeit only in the spherically symmetric ap- 
proximation, the same arrangement of energy levels (fig. 2)! This is to be conceived as the 
RST counterpart of the conventional energy degeneracy of ortho- and para-positronium 
when the magnetic interactions are neglected, see subsection V.l (Conventional Multiplet 
Structure). 

But the present RST picture of this ortho/para-degeneracy is somewhat more intricate 

122 



than its conventional competitor, namely because it elucidates also the fact that the 
ortho/para-degeneracy becomes broken not only if the magnetic interactions are taken 
into account but also if the spherically symmetric approximation is dropped. In the latter 
case the reduced wave amplitudes 'R±, 'S± (lVI.19ap - (1VI.19bl) for the para-configurations 
and their ortho-counterparts "R±, "S± ( lVI.llla]) - (1VI.lllbl) will depend also on the polar 
angles d and 0; and additionally the electric potential ( b ' p 'Ao(r) will become a function 
of r and i?, even if the magnetic interactions are kept inactive. Here it should be evident 
that such an angle-dependent charge density ^ko(r), as it is observed for the ortho- 
configurations (fig. 5), will entail a rather different interaction energy (of the electric type) 
in comparison to the "almost" spherically symmetric charge densities ^^(r) (~ of 
the para-configuration, cf. (lVI.17p . Thus there remains in RST the problem to clarify the 
mutual numerical relationship of the anisotropic and the magnetic corrections which are 
to be added to the predictions of the spherically symmetric approximation. 



Anisotropic and Magnetic Effects 

A first preliminary treatment of the anisotropic effects can be carried out most easily 
in the non-relativistic approximation. Recall here the remarkable fact that both rela- 
tivistic eigenvalue systems (IVI.l 13a|) - (lVI.113b|) and (lVI.115aj) - (1VI.115bj) have the same 



non-relativistic limit ( IVI.l 181) . However, observe also that this common limit form of 
the eigenvalue equations does not imply the numerical equality of the corresponding en- 
ergy levels if the anisotropic or magnetic (or simultaneously both) effects are taken into 
account. Actually, it is only in the electrostatic and spherically symmetric approxima- 
tion that both coinciding limit forms ( 1VI.118I) of the relativistic eigenvalue equations do 
imply the same energy levels. The reason is here that both relativistic versions (lVI.113a|) - 



(lVI.113bl) and (IVI.l 15a|) - (1VI.115b[) are to be linked to different angular parts of the charge 



density ^ b 'ko(r), even in the spherically symmetric approximation. 
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In order to present a simple example of this charge-density dichotomy of the ortho- 
configurations, one may consider the smallest possible value for the ortho-number i. e. 
£o = 2 (see the first line of the table given in appendix E). For this case the angular 
functions g^fi) (lVI.105ap and <7s(i?) (1VI.105bj) are given by 



g R (ti) = p 3 ■ sin 3 $ 

gs($) — &i • cosi? . 



(VI. 140a) 
(VI. 140b) 



Furthermore, the Dirac density ' b 'ko(r) is seen from equation (1VI.130I) to factorize into 
an angular and a radial part according to 



where the radial part is given by 

[b] k (r) 

and similarly for the angular part 



nl(r) + n 2 _(r) 



g_m 

sin 3 "d 



+ sin$ • g 2 s (d) 



(VI.141) 



(VI. 142) 



(VI.143) 



This latter part is the marked difference between both types of ortho-configurations and 
becomes with reference to the chosen simplest situation ( 1 VI. 140 a]) - ( I VI. 140bl 

'2 



W *b(*) =► ^ 
47T 



sin 3 •& H — §■ sin 1? cos 2 
Pi 



(VI. 144) 



The claimed density ambiguity of the ortho-configurations arises now through the 
circumstance that the ratio — can adopt two different values, namely for to = 2 either 



Ps 



-1 



(VI.145) 



for the first system (lVI.113al) -( 1VI.113bl) ; or alternatively for the second system ( 1VI.115al) - 
(IVI.llBbD : 

(VI. 146) 



^=3 

P3 
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see equations ( 1E.16jl -l lE.17j) and the first line of the table in appendix E. Accordingly, one 
obtains two different angular factors ^ko($) ( 1VI.143I) which both can be required to obey 
the separate normalization condition (!VI.131aP and do then appear as follows: 

1 (VI. 147a) 



"M^(.„ ....... 



-(l + 8cos 2 i?) , (VI. 147b) 

see fig. 5. Of course, such an angle- dependent density would generate also an angle- 
dependent potential ^A (r, $) but this becomes truncated in the spherically symmetric 
approximation to a symmetric potential (i.e. ^Ao(r, $) =>- ^Ao(r)) so that the angular 
discrimination between both types of ortho-configurations becomes ineffectual. 

But, strictly speaking, both states (due to the same value of but different ratio — ) 
must be expected to carry different energy E^, even under the neglection of magnetism! 
The reason is that this quantity Et reacts sensitively to any redistribution of the electric 
charge density ^ko(r), namely via the electrostatic interaction energy E$ ( IVLTOl) which 
may be rewritten in terms of the Dirac density ( b ' p >ko(r) as 

4^-//^¥-' W ^^l'y° (f '). (VI.148) 

Nevertheless, both charge distributions for given quantum number £q do generate the 
same electrostatic interaction energy (1VI.134h provided that one resorts to the spher- 
ically symmetric approximation] The reason is that the approximate Poisson equation 
( 1VI.139I) implies the validity of the associated Poisson identity ( 1VI.138I) with the electro- 
static interaction energy E^ being specified by equation (IVI.134I) and its mass equivalent 
M^c 2 by (IVI.137I) . However, neither of the latter two quantities E^, M^c 2 is influenced 
by the angular parts gn{i9) and gs{$) of the wave amplitudes as long as one is satisfied 
with the spherically symmetric approximation! From this reason one finally arrives at a 
multiple degeneracy phenomenon: namely (i) the mathematical equivalence of the para- 
eigenvalue problem { ( I VI. 23al) - (1VI.23b"l) . ( 1VI.37I) . ( 1VI.38I) } and the ortho-eigenvalue problem 
{ f lVI. 1 15aD - ffVI. 115bl) . rtvTTM . (IVI. 1391) } says that the solutions of both eigenvalue prob- 
lems must be the same and therefore must carry also the same energy Ey, and (ii) the 
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solutions of both ortho- sub systems flVI.113aj) - (lVI.113bj) and (lVI.115al) - dVI.115bl) do also 
carry the same energy (albeit only in the spherically symmetric and non-relativistic ap- 
proximation) . 

Furthermore, it should be obvious that one has to expect also a certain interplay 
between these anisotropic and the magnetic effects; for a more distinct anisotropy of 
the wave amplitudes will imply also a more marked anisotropy of the azimuthal current 
(IV.58c|) . For instance, the latter quantity becomes quite generally by means of our 
ortho-ansatz fIVI.lllajl - flVI.lllbjl 

% - {§g - si ° 2 ^w - 2c ° s ^i ' 9sW } • ™ 

and if we are satisfied for the moment with the simplest possible ortho-configuration (i. e. 
the first line of the table in appendix E with to = 2) one is led to the following result 
(see also fig. 6) 

~ , ~ , x . . 2 f 1 , — = -1 (VI. 150a) 

2r / -(l-16cos 2 tf), ^ = 3. (VI.150b) 

The key features of this ortho-current become obvious through a comparison to its 
para- counterpart (jV1.30p . see also fig. 1. The first observation is here that the 

para-current is finite (and therefore singular) on the z-axis ($ = 0,7r) whereas the 
ortho-current ^'k^ flVI.150a]) - flVI.150bl) vanishes for •§ = 0, it. This is physically plausible 
because the angular momentum J z of the para-current is zero, cf. equation (IV.46I) . which 
may be associated with a very close revolution of the quantum fluid around the z-axis 
(fig. 1). On the other hand, the ortho-current ( % (jVI. 150ajl - ( jVl. 1 50bj> is zero on the pr- 
axis and is maximal off the z-axis which intuitively is to be associated with a larger angular 
momentum; in fact, the eigenvalue of the present ortho-configuration ipQ corresponds to 
one quantum of angular momentum, cf. equation ( 1VI.119I) . Clearly, such a difference of 
the angular distribution of the currents will entail also a different magnetic interaction 
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energy of the para- and of the ortho- (sub) configurations, beyond the different relative sign 
of the individual magnetic fields Hi and H 2 due to H p and iff,, cf. the discussion of the 
magnetic energy difference A^E T (jV30jl - (jV3II ). 

But clearly it remains to be clarified whether those density ambiguities as for ^k 
flVI.1471) or for flVI.1501) are artefacts due to the product ansatz (IVI.lllajl - flVI.lllbh 
in connection with the spherically symmetric approximation, or whether that ambiguity 
is a real physical effect which would persist also for the exact solutions of the eigenvalue 
problem? 

Non-Relativistic Energy Functional 

For a first rough estimate of the energy content of the ortho-configurations it may 
again be sufficient to consider merely the non-relativistic approximation. Its deduction 
follows essentially those paths having been described already for the para-configurations, 
cf. flVI.50P ; however, there arise here some few points which are somewhat more subtle 
and therefore need a closer inspection. The main difference is in the present context that 
the range of the ortho-number Iq is restricted to Iq > 2, whereas the para-number i-p can 
adopt any non-negative integer; and furthermore the passage from formally £q =^ is 
associated with a distinctly different ortho-density ^fco(^) ? cf. (1VI.147al) vs. (1VI.147bl) . 
whereas the para-density (1VI.17I) does not react at all to the substitution i-p =>- 
—i-p, cf. the angular dependence of the para-amplitudes &R ± , V>S± flVl.25ap - flVI.27bp . 



These differences must now be elaborated in greater detail in order to better understand 
the ortho/para degeneracy and its breaking by the magnetic interactions and by the 
anisotropy of the electric interaction potential. 

It is true, both relativistic versions flV1.113ap - flVI.113bp and flV1.115al) - flVI.115bl) of the 



eigenvalue equations for the ortho-configurations are related to each other through the 
formal substitution Iq =>■ — io which then would correspond to the analogous replacement 
(i-p =>• —i-p) f° r the para-system flVI.23ap - flVI.23bp . The same change of sign would 
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then also apply for both kinetic energies lp] f kin = ^f [r) + Wf w flVI.32ajl - (IVI.32bjl and 
(6) T kin (1VI.133I) . but the case with the charge densities ^ b ^ p ^ko(f) is somewhat different. 
Namely, while the replacement i-p —t-p leaves invariant the charge density ^k as 
the source of the electric potential (p U , cf. flV1.17|> . the passage from the first ortho- 
case (IVI.113al) -( lVI.113"bl) to the second ortho-case (1VI.115al )-( 1VI.115bl is accompanied 
by an essential change of the ortho-density ^ko(r): its ang ular part Mfc o (0) (j VI 143ft 
receives a considerable change of the magnitude of both angular functions gn($) and 

gs{$)- This relative magnitude may be measured by the ratio of the coefficients — , 

Ps 

cf. equations ( 1E.16ft -( lE.17l) of appendix E. A special example of this effect is given by 
the equations (1 VI. 145[) - (1 VI. 147bj) . It is only on account of the spherically symmetric 
approximation that this change of the charge density ^k (r) leaves us without any effect 
since its angle-dependent part is integrated off for the spherically symmetric form of the 
energy functional, see for example the deduction of the mass renormalization factor 
(IVI. 132ft or also the deduction of the kinetic energies ^T r , ^T#, in appendix F. 

On the other hand, both relativistic ortho-systems ( lVI.113aT) -( 1VI.113bl) and ( IVI.llSaT) - 
( lVI.115bl) have the same non-relativistic approximation ( 1VI.118I) . albeit only within the 
framework of the spherically symmetric approximation. This non-relativistic limit form 
(TVI. 118ft of both ortho-versions is the same as the corresponding limit form of the para- 
versions, cf. ( 1VI.42I) ; and also the relativistic energy functional E\g>] (1VI.50I) and E[n] 
(TVI. 126ft are the immediate counterparts of each other. Finally, the electric Poisson equa- 
tions are also formally the same since the non-relativistic limit of flVI. 139|) 

^Uo(r) 2 #U (r) fi(r) 2 

— 1 j = -« s (VI.151) 

ar z r ar r 

is again the immediate ortho-counterpart of its para-predecessor fTVI.441) . 

Thus there remains the last problem, namely of setting up the non-relativistic version 
E[n] of the present relativistic ortho-functional E[a\ ( 1VI.126I) . Clearly, one expects here 
that this will turn out again as the immediate ortho-counterpart of its para-predecessor 
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E[$] ( 1VI.50[) : i.e. the expectation is the following 

E [n] = 2E kin + 4 C] + 2A S • No + A§ } • N§ ] (VI. 152) 

so that we are left with the problem to validate this assertion ( 1VI.152I) by deducing the 
non-relativistic forms of all the relativistic contributions to the functional Em ( 1VI.126I) . 



To this end, one will of course follow again the path traced out already by the analogous 
para-calculations, see equations (1VI.45I) - (1VI.49I) : 



(i) mass renormalization /[VI. 130) : 



2 o } 

2-3-2 _ o 1/ 2 ^ /" J„«. J / \ Iq ~ 2 I H fid / - 2 



(VI.153) 



(ii) jfcifieftc energy M T fcin ( fVT.j,^) : 

ft 2 . f . d r~, . \ 2 



M T kin => £ kin ± ^ i oJdr-^- (fi(r)) 2 (VI.154a) 




(iii) electric Poisson identity /iVI.138\) 



Ng =>• N£ J = ££ J - M [< V = . (VI. 155) 

Here the interaction energy is formally the same for both the relativistic and the 
non-relativistic situation, cf. (IVI. 134ft ; but for the non-relativistic form of the mass 
equivalent M^c 2 one neglects again the "negative" wave amplitude f2_(r) against 
its "positive" counterpart fi+(r), i.e. one puts 

M [e] c 2 M [e] c 2 = -He J drr [b] A (r) ■ (ft(r)) 2 , (VI.156) 
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(iv) normalization condition /IVI.lSlb]) : 



N n ^N n = J drr f2(r) 2 -1 = 0. (VI.157) 



Thus the non-relativistic ortho-functional Er^i 0VI.152P is uniquely fixed; and it is again a 



rather satisfying aspect that the boundary term (~ Eo) cancels through the cooperation 
of mass renormalization and kinetic energy: 

2Mc 2 Zl + 4 ^f kin =s> 2E kiQ , (VI.158) 



with the one-particle kinetic energy E^ n being specified by equation (lVI.154bj) and the 
rest mass 2Mc 2 being omitted. This cancellation is an essential fact for the non-relativistic 
limit because it just ensures its uniqueness both relativistic ortho-versions ( 1VI.126I) 
with ±£q in the kinetic energy ( 1VI.133P generate the one non-relativistic form 



( 1VI. 152ft with quadratic dependence upon £ , cf. ( 1VI.154bj) ). 

Finally, it is also a very satisfying exercise to identify both the former non-relativistic 
eigenvalue equation (1VI.118I) and the present Poisson equation (1VI.151I) as the extremal 
equations (5Ep = 0) due to the claimed energy functional Ep flV"I.152[) . This proves 



the mathematical equivalence of the non-relativistic eigenvalue problems for the ortho- 
and para-configurations, albeit only in the electrostatic and spherically symmetric ap- 
proximation. The physical consequence hereof is the ortho/para degeneracy as the RST 
counterpart of the conventional degeneracy (IV. 13|) . For the breaking of this ortho/para- 



degeneracy one first has to work out more exact solutions of the present RST eigenvalue 
problems, especially with regard to the anisotropy of the electric interaction potential 

(b/p)A (r). 
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Appendix A: 
Manifest Gauge Invariance of the Energy 

Functional E[ T ] 



As it stands, the proposed energy functional E[t] ( 1IV.20I) with its physical part E^ V ' 



[T] 

( 1IV.23I) is surely not manifestly invariant with respect to the magnetic gauge transfor- 
mations (IIII.20ah - f1III.20cj) . because there are two kinds of contributions which lack that 
desired invariance: these are (i) the kinetic energies Tki n ( a ) (a = 1,2,3), as they are de- 
fined by equations (1IV.26I) -( 1IV.27I) . and (ii) the magnetic constraint term (1IV.5I) . 



Perhaps a somewhat critical term (concerning the gauge invariance) is Uq (IIV.6P which 



contains the integral of the scalar product of 21 and I and transforms as follows: 

J d 3 f{^'l) ^ Jd 3 r(W-l) = J d 3 r(K-T) - J d 3 fV(a 2 -a 3 )-t. (A.l) 
Thus the gauge invariance demands the vanishing of the last integral 

J d 3 fV(a 2 -a 3 )'f = - J d 3 f (a 2 -a 3 )V«f = , (A.2) 



and this would ensure the gauge invariance of tiq (IIV.6P . This requirement (IA.2p , how- 



ever, is satisfied for the stationary bound systems. Namely, these systems are to b( 
expected to have source-free currents, cf. ( 1II.18b[ )-( 1II.18c[) 

V^ aM = - = , (A.3) 



and this assumption lets vanish the source of the entanglement current (IIII.47P for the 
stationary situation 

VZ M = V-f=0 (A.4) 

which obviously validates the invariance requirement (1A.2I) . 

Now it should be self-evident that the claim of invariance of Em can be true only if 
the sum of the two non-invariant terms (i) and(ii) is itself manifestly invariant. This, 
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however, can easily be proven by first decomposing that magnetic constraint term Nq ^ 
(1IV.5P into a gauge-dependent part Nq^ and a manifestly gauge- independent part Nq , 
i. e. we put 



N, 



(m) 



G 



with 



N, 



(m) 



G 



•y / dV ( .V./i +^j,vy 2 + .4, n vy, 



A, 



(m) 



d 3 f 



(; 4vra 

+ f V x A 3 ) • ( V x Ax 



VxAi • V x A 



Vxi 2 * Vxi 3 



(A.5) 

(A.6a) 
(A.6b) 



But here the Maxwellian three-currents j a f (a = 1, 2, 3) in (IA.6aj) are identical (up to sign 
to the Dirac currents k a r ( lII.17al) -( 1II.17cl) so that one has (e.g. for the second particle) 

"(V+^V- + ( V-<? (2 V+ 



drfAa' 



(A.7) 



On the other hand, the kinetic energy T kin ^ 2 ) ( 1IV.27I) reads in terms of the Pauli spinors 



(2) 



tp±r 



T, 



kin(2) 



-i — I d 3 r 
2 



(A.f 



(V^-v (2 V+ + (2 V^*v ( V 

so that the sum of both non-invariant terms for the second particle (a = 2) actually 
appears now in a manifestly gauge invariant form 7k in ( 2 ) : 

I d J fAu']2 

(A.9) 



%\n{2) = T kill{2 ) - — J d fA U 'j 2 



-i — I tf'r 
2 



Wplff. (v + a„)( 2 V+ + {2 V^- (v + ^n) { V 

When the first (a = 1) and third (a = 3) particle are treated in the same way, the 
energy functional emerges finally in its manifestly gauge-invariant form: 



E m = Z( 1} ■ M p c 2 + Z( 2) ■ M e c 2 + Z( 3) ■ M e c 2 + 2T kin 



+ £ A D(a) • iV D(a ) + Ag } ■ (N$ - 4 e) ) + X^- (N™ - 4 m) ) + A£ } • (N^ - n^) (A.10) 



(m) / ," r (m) 



,(m)v 



0=1 
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Here, the first line constitutes now again the physical part of the RST energy Et (which, 
however, is now manifestly gauge invariant) whereas the second and third line represent 
the constraints with being replaced by its invariant part (lA.6bl) . But as a 

consequence of this rearrangement, the residual magnetic constraint term — 
contributes now non-trivially to E^]. 
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Appendix B: 



Variational Deduction 
of the Non-Relativistic Eigenvalue Equations 



The detailed variational deduction of the RST field equations is an instructive exercise 
since it mediates a better feeling of how to handle the magnetic gauge structure. In 
order to present a brief example of this, consider the deduction of the Pauli eigenvalue 
equations ( lIV.40al) - (1IV.40cl ) . First recall here the fact that all the gauge field energies 
(IIV.28p - flIV.3ip do not depend at all upon the matter fields ^^p±(r); and therefore the 
corresponding variations (<5 mat , say) of the energy functional Em 0IV.48P with respect to 
the conjugate Pauli spinors ^tp^ becomes simply 



^mat^[T] — #mat 



E kin + hcj d 3 r [Hj • H£ + e u • ^S 2 + E IU • H£ 3 } 



+ ^A P(a) -N D{a) + Ag ) -Ng ) 



a=l 



(b.i; 



+ 5. 



mat 



A^-f^ + AifNg) 



Here, carrying out the variational procedure for the first part proceeds along standard 
rules and yields 

d 3 f{(5( 1 V+) t [W(7) + A P( i)l 



Sl 



mat 



+(* (2 V + ) 



(e) . pj(c) 
t 



+ Ap(2) 



+ Ap(3) 



(3) 



(B.2) 



where the one-particle Hamiltonians (7i) have already been defined through equations 
( lIV.42aP -f lIV.42cp . and the previous value ( lIV.22bp for the multiplier Aq^ has also been 
used. 

However, the variational formalism for the second part (i. e. exchange part in IB. 1[) 
requires a more thorough discussion. It is true, the "electric" part (Nq"*) hereof presents 
no problem because the matter fields ^ip±{r) enter the exchange density ho(r) in the 
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form [13] 

h = ^ °V + ( V- (3 V- (B.3) 

and thus its non-relativistic approximation h.o(r) is obtained by simply neglecting the 
"negative" Pauli spinors ( a )(/?_: 

h =► ho = (2 Vl (3 V+ • (B.4) 

Therefore the variation (t^at) °f the non-relativistic constraint term Nq ( 1IV.12I) with 
respect to the spinors ^(p + becomes 



(B.5) 

d¥ 5 ■ (5 (2 V + ) t(3 V + + £o • (5 (3 V + ) t(2 Vn 



Tic 

~ ~ y „ 

But the variation of the magnetic constraint term Nq (in the second line of (IB.ip ) 
represents a more subtle problem. Indeed, here it is necessary to first split up the non- 
relativistic approximation lh of the exchange current h into a convection part b and a 
polarization part is, i. e. 

h^h = h + z , (B.6) 

see the preceding paper [13] for a more detailed discussion of this splitting. The convection 
current b is given by 
ih 

h = 

2M e c 

and the polarization current % emerges as the curl of the magnetic polarization density Z 



{ [(v + zln) (2 V + ] 1 (3 V+ - ( (2 V + ) t (V + ^„i) (3 V + } (B.7) 



i = VxZ (B.8a) 

2 *^bWV+. (B.8b) 

Clearly as a consequence of this splitting, the exchange mass equivalent M^c 2 of the 
magnetic type (see ref. [13] ) 

M (g) c 2 = — [ d 3 r <J B-h + B-h )■ , (B.9) 



2 
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which for the present intention is the relevant part of the constraint term Nq (1IV. 1 5j) . 
must also be split up into two components 

M (g) c 2 = M^ nv c 2 + M^c 2 . (B. 10) 

Evidently, the convection part M 

conv C cipp GclXS here as 

M^c 2 = ^ J d i f{B-h + B>h\ , (B.ir 
and similarly for the polarization part M^c 2 



T(g)^2 _j_ 



he 



c z = — d'ri B'% + B'v\ . (B.12 



P oi- • 2 

The reason for such a close inspection of the mass equivalent M^c 2 is of course that 
the matter variation of the constraint term N^p (IIV. 15[) is identical with that of the mass 
equivalent: 

= -iat (M (g) c 2 ) = -4 at {MiHc 2 ) - 5l at (M^c 2 ) . (B.13) 

Turning here first to the convection part McolwC 2 , one recasts this term through integration 
by parts into the following gauge-invariant form: 

(B.14) 

V'-B + 2B-(V + iA u ) 



ih2 f d*f(6^v + y 



AM, 



(2 Vh 



Next, consider the non-relativistic polarization part M^c 2 ( IB. 121) . Since the exchange 
polarization current % ( lB.8al) is the curl of the exchange polarization density Z, one 
can convert the polarization part (through integration by parts) to the following gauge- 
invariant form 

4at(Mgc 2 ) =^f d B f{(^)W<*V + + (*< 3 V + ) W< 2 V + } • (B.15) 
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Finally, one puts together the variation of the second part of ^ at E[T], in the second 
line on the right-hand side of equation ( 1B.1I) , by means of the partial results (|B~5l) . (lETTIl) 
and (1B.15P and thus one finds 



Si 



A "mat 



d 3 f (^ 2 V + ) T W 3 V+ + (* (3 V + ) T 



(B.16) 



Here the exchange Hamiltonians (()) have already been defined through equations (1IV.43al) - 
(lIV.43bj) and the Lagrangean multipliers Aq \ Xq are taken from equation (IIV.22cl) . But 
now that both essential parts ( IB. 21) and ( IB. 161) are explicitly known, one arrives at 



<W E [T] 



d 3 f{ (5 (1 V+) t W(J) + A P(1) (1 V++(^ (2 V + ) T (7i:(//) + Ap (2) )( 2 V + 



+ ^ (// ) (3 V + l + (5 (3 V+) + + Ap (3 )) (3 V + + W 2 W 



(B.17) 



If this variation is now put to zero according to the RST principle of minimal energy 



4at%] = 



(B.18) 



one just ends up with the claimed eigenvalue equations (lIV.40ap - (1IV.40cp . provided one 
identifies the non-relativistic Lagrangean multipliers Ap( ) with the Pauli energy eigenval- 
ues E P{a) (ITVTO -fllVTO 

Ap(a) = —Ep( a ) . (B.19) 

In this way, the logical consistency requirement of commutativity on p. EH is actually 
realized. 
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Appendix C: 

Variational Deduction of the Gauge Field Equations 



The set of non-relativistic gauge field equations is chosen to be formally the same as in 
the relativistic case, with the exception that for the charge and current densities one has 
to substitute their non-relativistic forms [TJ] . Therefore one can deduce the desired gauge 
field equations from the relativistic energy functional E[t] ( 1IV.20I) rather than from its non- 
relativistic approximation Em ( 1IV.48I) . Here it is advantageous to consider the electric, 
magnetic and exchange subsystems separately because any subsystem does require for its 
treatment only a few terms of the whole functional Em. 

Thus turning first to the deduction of the electric source equations (1III.29ah -f 1III.29cl) . 
one denotes the variation of Em with respect to ( a l4 (r) by ^5q ( a — 1, 2, 3) and then 
finds for the first gauge field mode ^A 

( %(£ [T] ) = { %(4 C) ) + Ag } - ( %(4 e) ) = , (C.I) 
whereas for the second and third modes ^Ao(r) (a = 2, 3)one has 

{ %(E [T] ) = < a fc (4 e) )- (a ^^ = • 

(C.2) 

Obviously the first mode (a = 1) contributes to the variation of the functional Em only 
two terms, cf. (IC.ip . whereas any of the identical particles (a = 2,3) contributes five 
terms, cf. ( 1C.2I) . The resaon for this is that the second and third potentials ^ 2 Aq and ( 3 ko 
(but not the first potential ^A ) do influence the electric exchange energy 
via the exchange field strength X fllll. 25a^ ; and furthermore these potentials enter also 



both constraint terms (1IV.12I) and (1IV.15D via 2l and 2t M - As a consequence, 
the (generalized) Poisson equation for the modified mode ^Aq ( 1III.8a[) is found from ( IC.lf l 
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to look relatively simple, i. e. 

(C.3) 

( (I) Jo = (2) Jo + (3) Jo) , 

whereas the occurence of the exchange effect in the variational equations ( 1C.2I) lets appear 
the corresponding Poisson equations for the other two potentials and ( m ko somewhat 
more intricate: 

A( n k + 2V • W = -4vra s ( (II) j - k) (C.4a) 
A (m k - 2V • W = -4vra s ( (III) j + k) (C.4b) 
( (II) Jo = (1) Jo + (3) Jo, (m) Jo = (1) Jo + (2) Jo) • 

It is true, at first glance these electric Poisson equations may appear somewhat unaccus- 
tomed, but by means of the equations for W ( \lll.52ah and Iq (1III.57H it is easy to show 
that the present Poisson equations for ®Aq, ^ u A , ^ iu A (1C.3[) - f1C.4b|) are actually equiva- 



lent to the original source equations ( 1III.29al) - (1III.29c l ) for the electric field strengths E a 
(a = 1,2,3). 

A similar situation occurs for the case of the magnetic Poisson equations. Denoting 
here the variations of E[t] with respect to the three- vector potentials A a by one finds 
again that the variation with respect to the first mode A\ yields only two terms, i. e. 

% (E m ) = -% (E^) + ( % (Nt ] ) , (C5) 

whereas the variation with respect to the other two potentials A a (a = 2, 3) produces six 
contributions 

(E m ) = -(% (E^) - % (i#>) + (% (Eg) + 1% (N™) 
+ • ( jvg° ) + Ag } • (N® ) . 

This corresponds to the circumstance that the energy functional Em is gauge-invariant 
with respect to the magnetic transformations, namely through the use of the covariant 
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derivative V' ( 1III.26I) which relies on both vector potentials A2 and A3, but not on A\. 



As a consequence, the first magnetic Poisson equation emerging from ( 1C.5I) looks again 
rather simple 

AA/ = -47ra s j/ 

. . . ( C - 7 ) 

(j/ = J2 + js) , 

whereas its second and third counterpart (a = 2, 3) necessarily must appear somewhat 
more complicated 

AA n - 2(V x V) = -4rra s (j n - f) (C.8a) 

A4j + 2(V x y) = -4vra s (j m + f) (C.8b) 

(In 4= ii + J3, jm = Ji + J2) • 

Here it is again a simple exercise to show that the present (generalized) Poisson equations 
( 1C.7p - (lC.8bj) are actually equivalent to the original source equations (1III.30al) - (1III.30cl) 
(hint: use the equations flIII.52bD and (1IIL58H ). 

The situation with the exchange counterparts of the electric and magnetic Poisson 
equations needs also sufficient attention because here the electric and magnetic field 
strengths E a and H a are again influenced by both the electromagnetic and exchange 
potentials; see equations (1III.23bl) - (1III.23cj) and (1III.24bl) - (1III.24cl) . Keeping this in mind, 
one finds the following five contributions to the variation (^£0, say) of Em with respect 
to the electric exchange potential B : 

iB) 5 (E m ) = n (E^)-n (E^) + \^5 (N^) . 

(C.9) 

Carrying here through all five variations yields the electric exchange equation in the 
following form: 

£B + %l l B -2iB*(£-W) = -4ira s h (CIO) 

where the gauge-covariant Laplacean A' is nothing else than the square of the corre- 
sponding gradient operator V (IIII.26I) . i.e. A' = V* V. Clearly, the present exchange 
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equation ( 1C.10I) is again equivalent to the source equation ( 1III.31al) for the exchange field 



strength X. This equation is also well-suited in order to demonstrate that oversimplified 
linearization procedure of ref. [13]. Indeed, truncating the exchange equation (IC.lOp by 
means of the replacements 

A' => A (C.lla) 

*m=>A (C.llb) 
am 2 

B'(i-W) =4>0 (C.llc) 
lets appear that equation ( 1C.10I) as 



AB + -±-B = -Ana s h , (C.12) 
am 2 



and this is just the equation (III. 10a) of ref. 

But clearly, it should be self-evident that the most intricate case refers to the varia- 
tional deduction of the curl equation ( lIII.31bl) for the magnetic exchange field strength Y . 
Indeed, one finds that the variation (^5-[, say) of Em is built up by seven contributions: 

[Aj. 13 ) 

+ A§> • + A&> ■ (N^) + A« ■ (4 S) ) • 

However, carrying properly through all the variations yields finally the expected curl 
equation (lIII.31bj) . provided the Lagrangean multipliers adopt those values being specified 
by equations ( 1IV.22bl) -( 1IV.22cj) . Moreover, the somewhat oversimplified equation (III. 10b) 
of ref. p3] for the magnetic exchange field B can of course be deduced from the present 
correct equation (1III.31bl) by destroying the magnetic gauge invariance in the following 
way: 

2l M X^0, B £ + iBx$=>0 (C.14a) 
V' xy4Vx^Vx(VxB) = —AB + V(V • B) 

-AB + — V5 -AB l —B. (C.14b) 
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Of course, this chain ( lC.14bl) of neglections ( "=^" ) demonstrates clearly in what way 
magnetic gauge invariance becomes spoiled by the various steps of "linearization" . 
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Appendix D: 
Electric Poisson Identity 

The computation of the gauge field energy ( 1VI.681) or its mass equivalent 
(IVI.69P is much more intricate than the analogous case with the kinetic energy E^ in 
( 1VI.65|) and therefore the numerical equality ( 1 VI. 701) . as the immediate consequence of 



the electric Poisson identity, may be used as a welcome check for the correct calculation 
of the interaction energy. 

Considering here first the part of the gauge field energy E$ VI. 680 . the associated 
potential function £ po t(^) ( 1 VI. 701) arises in the form 

which requires to first determine the field strength in dimensionless form. Trying for this 
the following ansatz 

^ = -1 {l My)) e-«}, (D.2) 

the equation for the ansatz function f u (y) is to be deduced from the Poisson equation in 
the dimensionsless form ( 1V"I.71j) and thus is found as 



dl(y) ff , . y 2u+1 m ,x 

In order that the field strength ( 1D.2I) . i.e. the derivative of the potential a v (y), remains 
finite at the origin (y = 0) one demands as an initial condition for f v {y) 

W) = , (D.4) 
and then the solution for f v (y) contains the variational parameter v in the following form: 

oo n 

Uy) = e y -l- y 2v+2 ■ Y V - , (D.5) 

n=0 v ' 
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which simplifies for integer values of 2v to 

y 



n 

71=1 



This solution may now be substituted back into the field strength ( ID. 21) which in the 
general case (ID.5[) then reappears as 

^ = _e-»£_jg_ (D . 7) 

Obviously the electric field strength is finite at the origin (y = 0) for v = but is zero 
for v > 0, see fig. 1 of ref. [19]. Furthermore, by use of the result (ID .71) the potential 
function e po t{v) ( ID. II) becomes now 

= _L_ ^ _L_ r(4z/ + 3 + m + n) 

£ potl^J o4l/+3 2^ om+7i TV 0;/ _|_ Q _|_ rn^\ ■ Pf 9zy -1-3-1- [ J 



2^+3 2 m + n r(2z/ + 3 + m) -r(2i/ + 3 + n) ' 

m,n=0 



For integer values of 2v (i.e. 2v = 0, 1, 2, 3, . . .) this potential function becomes simplified 
again considerably: first, one recasts the original definition (ID.ip by means of some simple 
manipulations into the following form 



2 / x 2 



^ 2 (^r) S (,!l( *vt) -S (hl[< (d.9) 

and if the special function /^(y) ( ID. 61) is substituted herein, one arrives at |13J 

2 ^±,' 1 (,„ + „ _ 2)! ^ 1 / IN , 

m,n=l n=2 V / \ / 

Thus for the lowest values of v one gets the following table (for v — the second sum in 
( ID. 101) is to be omitted): 



V 





1 

2 


1 


3 
2 


2 


5 
2 


3 


7 
2 


£pot(v) 


i 

2 

=0,500 


5 

16 
=0,3125 


11 

48 

=0,2291. . . 


93 
512 

=0,1816. . . 


193 
1280 

=0,1507. . . 


793 
6144 

= 0.1290. . . 


1619 
14336 

=0.1129. . . 


26333 
262144 

= 0.1004. . . 
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But as the dimensionsless Poisson identity ( 1 VI. 7011 says, the positron-electron interac- 
tion energy may also be viewn to be represented by the mass equivalent M^c 2 ( 1VI.69I) 
and then the potential function £ po t(^) appears in the following alternative form 



1 

r(2z/ + 2) J ^ Vtf/ ^ vy/ T(2z7+2) 



-potM = T^TTT^ I dyya v (y) ■ $(y) 2 = ^ ■ \ dyy 2v+1 e~y ■ ~a v {y) . (D.ll) 



If one wishes to calculate the potential function £pot(^) i n this alternative way (in contrast 
to (ID. Ill ), one obviously has first to calculate the electric potential a u (y). The differential 
equation for a v (y) is given by (ID. 21) . with f v (y) given by (ID. 51) in the general case, or by 
(ID.6P for integer values of 2u, resp. However, for elaborating the solution a v (y) of (ID.2I) 
for given f u (y) one has to fix the value of the potential a u (y) at the origin (y = 0) which 
works as an initial condition for that first-order equation ( ID. 21) . But how to fix that initial 
condition a v (0)l 

The solution of this problem comes from the standard potential [ p i4 (r) (1VI.39I) which 
reads in the non-relativistic approximation 

M Mr) = ^f^ (D . 12) 

or, resp., in dimensionless form for the selected trial function $(?/) (1VI.58I) 

aJy) = — / d 3 y' y \ ,, . (D.13) 

vKy) 4vr ■ T(2u + 2) J y \\y-y'\\ 

Here it is obvious that this solution has the right limit form (i.e. Coulomb potential) in 
the asymptotic region (y — > oo), i.e. 

lim 5„(y) = - (D.14) 

y^oo y 

cf. ( 1111.6011 ; and furthermore its value at the origin (y = 0) is found by explicit calculation 
of the integral in ( ID. 131) as 
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But with both boundary conditions ( ID. 141) - ( ID. 151) being thus fixed, the solution of 
equation ( ID. 71) for a v (y) is unique and looks as follows: 



+ 1 I ^ 

L n=0 



Mv) = ~t~\ { 1 - e ~ y ■ > ; wo " . x y 2v+n ) ■ ( D -!6) 



2z/ + l ^T(2z/ + 2 + n 

^ n=0 

For integer values of 2v i.e. (2z/ = 0, 1, 2, 3 . . .) this result is converted to 



n=0 y \ n=0 



(for z/ = the first term is to be omitted). And finally, the potential function e po t(' / ) 
( ID. Ill) emerges by means of the present potential a v (y) ( ID. 161) as 



1 ( 1 
■ + 1 I 2 4iy + 2 ' ^ 2 n 

\ n=0 



r(4i/ + 2 + n) 



(D.18) 



1v + 1 \ 2 4 ^+ 2 2 n r(2f + 2) • T(> + 2 + n) 
or for integer values of 2v, resp. 



1 1 / 1 ^ (2z/+l + n)! _ (2i/ + w)! \ 

£pot( ^ ~ 27TT + (2i/ + 1)! I 27TT ^ n! 2 2 -+ 2 +™ ~ ^ n! 2 2 -+ 1 +™ j ' ( ^ 

v 7 \ n=0 n=0 / 

Clearly, this latter result can equally well be used in order to generate those numbers 
which are displayed in the table below equation (1D.10I) : for a plot of £ po t( z/ ) see fig. 3. 
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Appendix E: 



Angular Part of the Ortho- Wave Amplitudes 

< 6 te ±) < 6 >s ± 

The angular dependency of the wave amplitudes must necessarily transcribe to all the 
physical densities which originate from these amplitudes. Therefore it is worthwile to 
inspect somewhat more closely the angular functions gR{$), gs{^) of the product ansatz 
flVI.lllaT) - (1VI.lllbl) . Here it should be clear that the physical difference between the 
ortho- and para-configurations will originate from just the different angular structure of 
the corresponding densities. For instance, the ortho-analogy ^k^ of the para-current ^k^ 
( 1VI.30j) will surely appear in a rather different way because the angular dependency of 
<7.r($) and gs($), as the solutions of the coupled system ( 1VI. 106aj )-( 1VI. 106bl) . cannot be 
of that simple form as their para-counterparts and fs($) (lVI.26ap - (1VI.26bj) . Quite 



generally speaking, both ortho-functions <7h($) and gs{$) do enter the ortho-current ^k 4 
in the following way: 



n sin^ J r r 

cot * SrWss W • gtWg^HgzMgtM . (E . i; 



2vr 

Indeed, this result is obtained by simply substituting the transformations \^ b 'TZ±, ^ b 'S±j 
{ "R±i "S±] ( 1VI. 1 1 laf )-( 1VI. 1 1 lb|) with general angular functions gR($), gs($) into the gen- 
eral expression (]V.58cp for the ortho-current ^'k^. Consequently, the wanted angular 



dependency of the ortho-current ^k^ is determined essentially by the (yet unknown) 
functions gR($) and gs{$)- This becomes more explicit by resorting to the identifications 
(IVI.112I) which yields the following product structure of the ortho-current w/ty (IE. II) : 

% - • { {^) - m - ^ > • <-) 
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Obviously, this ortho-current ^k^ displays an angular characteristic quite different from 
its para-counterpart ^'k^ ( 1VI.30I) ! 

The determination of the angular functions can be performed most conveniently by 
resorting to the decoupled second-order equations (1VI. 108aP - (1VI. 108bl) . Here it is advan- 
ageous to introduce the auxiliary variables x and z by 

x = sini) (E.3a) 
z = cos$. (E.3b) 

This is to be combined with the hypothesis that <7.r(i?) depends exclusively upon x and 
<7s(i?) exclusively upon z, i. e. we put 

g R (#) = G R (x) (E.4a) 
gs(#) = G s {z) . (E.4b) 



This arrangement recasts the second-order system (1VI.108aT) - (IVI.108bl) to the correspond- 



ing system for the new functions Gr(x) and Gs(z): 

(1 _ z2) _ ^ dGs^ = (ldo) Gs(z) (E gb) 

The general shape of the latter equations suggests now to try some ansatz in terms of 
a power series expansion, i. e. we put 



G R (x) = Y,PnX n (E.6a) 

p=3 

ns 

G s (z) = J2°nZ n (E.6b) 

n=l 

with constant coefficients p n and a n . For these coefficients one deduces by means of 
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standard techniques the following recurrence formulae: 

{to ■ 'to) + n(n - 2) 
Pn+2 = („ + 2)(n-l) ' Pn (E ' 7a) 
(i -*o)+n(n + 2) 

°n+2 = 7 —ZTJ — TT ■ (E.YbJ 

(n + 2)(n + 1) 

Since the proposed power series expansions ( ]E.6al) -( lE.6bl) are adopted to stop at certain 
maximal integers n R or ns, resp., there arises a halt condition for each series, namely 

(to ■ to) = ~n R {n R - 2) (E.8a) 
(i -t ) = -n s (n s + 2) . (E.8b) 

Thus the maximal powers n R and ns are linked to each other by 

n R = n s + 2 (E.9) 
(ns = 1,3,5,7, ...) 

and the allowed values for the product of the constants £q and Iq turn out as 

(to ■ to) = -n R ■ n s => -3, -15, -35, ... (E.10) 

A possible physical meaning of this product of the constants is seen more clearly by 
building the quantum number of orbital angular momentum to'- 



{to ■ to) = (n s + I) 2 => 4, 16, 36, ... (E.li; 



Naturally, it is very temptive to interprete this in the way that ortho-positronium 
owns an orbital angular momentum of magnitude toh (=> 2h, 4h, 6h, ...). But this would 
just meet with the intuitive expectation that each positronium constituent owns the same 
integer number of units (H) of (orbital) angular momentum. However, one should not 
forget here that, in RST (as a fluid-dynamic theory), the physical observables are obtained 
by integration of the corresponding densities, not as the eigenvalues of certain Hermitian 
operators. Therefore such quantities as the above-mentioned discrete numbers £q (IE. Ill) 
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and their para-counterparts l-p ( 1VI.20al) -( 1VI.20bl) . or similarly the energy eigenvalue E* 
in the Schrodinger-like eigenvalue equations ( 1VI.42I) and ( 1VI.118I) . are to be conceived as 
auxiliary objects which help to build up the true physical observables (e. g. the value 
of the energy functional upon some solution of the RST eigenvalue equations). 

At first glance, the spacing of the ortho-numbers £q (=>• 2, 4, 6, ...) is twice the spacing 
of their para-analogies i-p = 0, 1,2,3, ... (fig. 2); but a more thorough scrutiny reveals a 
certain dichotomic character of any value of £q so that the total number of physical 
states is (roughly) the same for the ortho- and para-configurations. In order to realize 
this dichotomy of the ortho-configurations more clearly, reconsider the definition of the 
ortho-number £q (IE. lip and combine this with the equation (lVI.117bl) 

£o+to = -2 (E.12) 

so that the two constants £q and £o, being originally introduced by equations (lVI.106ap - 
(IVI.106bl) , become fixed now in terms of £q in the following dichotomic way, namely either 



£ a = £ a - 1 (E.13a) 
£o = -(£o + l) (E.13b) 

or the other way round (substituting £q =>- — to) 

£o = ~(£o + l) (E.14a) 

e = e -i. (E.i4b) 

Indeed, both equations ( IE. IIP and (1E.12I) are obeyed by each of these two possibilities. 

But once the origin of the claimed dichotomy has become clear, one can look now 
for the physical consequences. From the mathematical viewpoint, one first wishes to 
see the integration constants o\ and p% of the solutions ( 1E.6al) -( lE.6bl) being completely 
fixed. But this is not possible because one integration constant (p3, say) must remain 
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free for the sake of normalization of the angular part of the wave functions (see equa- 
tion ( 1VI.131aT) ). Therefore it is possible to nail down only the ratio This, however, 
can easily be done by substituting the solutions gn($) and gs($), being associated with 
the auxiliary functions Gr(x) and Gs{z) (!E.6aP - (1E.6bl) . back into the original first-order 
system (IVI.106aT) -( lVI.106bl) which then yields 



1± 

P3 



-3 
+3 



i 


o 




3 


p 
c o 




1 


To 


( 


o 




3 


p 




1 


~T 



for ns 
for ns 



1,5,9, 



3,7,11, 



(E.15) 



Furthermore, the pair of auxiliary constants £q, io has already been shown to occur in 
two versions, cf. flE.13aj] - (IE.14bjl . and this of course transcribes then also to the ratio of 
integration constants (1E.15j) . i.e. 



1l 

P3 



3 



for £ = 2,6,10, 



(E.16) 



and 



+ 



+ 



P3 



3 



for i = 4,8,12, 



(E.17) 



'o 



- 1 



The upper case of each of both equations (IE. 16[) - (1E.17[) refers here to the first eigenvalue 
system f rVT113a[) - (lVI.113bl) ; and the lower case to the second system ( IVI.115al) - (1VI.115bj) . 
The subsequent table displays a collection for the three lowest-order cases io = 2, 4, 6. 
Observe here that for the pairs £o, (third column) and for the ratios ^ (last column) 
one obtains always two different possibilities for a fixed value of Iq which entails the 
existence of two different physical states for any admitted io- In the non-relativistic 
limit, both relativistic systems ( 1VI.113aT) -( 1VI.113bl) and ( 1VI. 1 15aT) -( 1VI. 115bl) collapse to 
the same non-relativistic equation (jVI.llSj) . 
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n s 


n R 


to 


io 




G s (z) 


G R {x) 


£1 

P3 


1 


3 


1 


-3 


2 


U\ • z 


P3 ■ X 3 


-1 


-3 


1 


3 


3 


5 


3 


-5 


4 


a\-(z- 2z 3 ) 


p 3 • (x 3 - fa; 5 ) 


3 
5 


-5 


3 


-1 


5 


7 


5 


-7 


6 


ai -(z-fz 3 + fz 5 ) 


p 3 -(x 3 -f X 5 + f X 7 ) 


3 
7 


-7 


5 


3 
5 
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Appendix F: 

Kinetic Energy ^Tkin of the Ortho- Configurations 



The most intricate contribution to the energy functional E\a\ ( 1VI.126I) is surely the 
kinetic energy ^T>i n . First, the kinetic energy ( 1IV.26I) -( 1IV.27I) reads in terms of the Pauli 
spinors ^ b '<p±(r) f° r the ortho-configurations 

he 



{b) T kiQ = i— I d 3 r 



(p+ (r) + {b) ip\ (r) a • V (b) (r ) 



(f.i; 



and when the Pauli spinors are decomposed here as shown by equations (lV.50ap - (lV.50bl) 
the kinetic energy appears as a sum of three contributions 



(6) T kin = ®T r + WT, + {b) T d 



(F.2) 



The radial part ^ b 'T r itself consists of two contributions where the first of them refers to 
the wave amplitudes ( b ^JZ±(r) and the second to the amplitudes ^S±(r): 



®T r = (fc) T r [ft] + {b) T r [S] 



(F.3) 



with 



{b) T r [JZ} 



he 



16tt 



d 3 r I 



+3 



dr ' dr 



dr 



Or 



(FA) 



and analogously for "T r [S] where merely the replacement 1Z S must be applied. Just 
as the radial part ^T r is a collection of all derivative terms with respect to the radial 
variable r, the longitudinal term collects the derivative terms with respect to the 



153 



angle &. 



he f d 3 



16tt I r 



-2i(f> 



(F.5) 



Finally, the azimuthal part represents a collection of all the derivative terms with 
respect to the azimuthal angle and looks as follows: 



7r / r 



+^ f Seattle* 



167r / r 



d(j> 



+e 



-2i</> 



90 



(F.6) 



In the next step, one substitutes herein the spherically symmetric approximations 
(IVI. 105a|) - (IVI.105b[) for the wave amplitudes ^ b ^JZ±, ^S± and thus obtains for both radial 
parts ^T r [R] and ®T T [S\: 



®T r [K]=>®T r [B!'] = ~ [d 3 r 

07T J 



9RW 

V sin 3 $ 



"R+(r) d"R(r) "RJr) d"R+(r) 



+3 



dr r 
"RAr) ■ "RJr 



dr 

(F.7a) 



«T r [5] =► (6) T r [S"] = / ^ (v/sTn^-^W) 2 - 



2 f (r) d "S_ (r) "S_ (r) d "S+ (r) 



r dr r 

„ "£+(r) • "5_(r) 



dr 

(F.7b) 



Now one resorts to the ortho-identification of the wave amplitudes ( 1VI.112j) and thus finds 
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the radial kinetic energy ^T r ( 1F.3|) in the following product form: 

2 



(b) T ^ (b)f 



he f dQ 
~27 4^ 



9R&) 



drr 



V sin 3 $ 



+ ( Vsintf • g s (d) 

Qjr) d&Jr) 



(F.8) 



n+(r) ■ fi_(r) 



dr 



dr 



But this is indeed now a very pleasant result because the product structure allows us to 
apply here the separate normalization condition flVI.131al) for the angular factor which 
then leaves the radial kinetic energy in the final form 



(b)f ^ [b)f 



he 



drr 



n+(r) dQJr) Qjr) dn,(r) 



tt + (r) ■ fL(r) 



dr 



dr 



(F.9) 



In quite a similar way one may calculate the longitudinal kinetic energy flF.5p 
and one finds the result again in form of a product structure, i. e. 



(b) T ^ [b]f = hc 



J 4tt |_Vsin 3 tf dd V ) dd Vvsi 



sin i? 



dr r 



(F.10) 



For the angular part one obtains here by use of the first-order angular system (lVI.106al) - 
flVI.106bl) 

d ( gd$) 



g R {V)-g s {$) sk&) 



dn . 

— < 2 COS 17 ■ 

An I sin -8 



sin d 



V sin 3 $ J . 



(F.ll) 



o 



sin 



+ sintf-^Ki/ 1 ) 



The term containing the quantum number (of orbital angular momentum) would 
present no problem because one could apply again the normalization condition (1VI.131al) 
for the angular part of the charge density ^ b 'ko(r) (1VI.130I) . (The upper sign in front of Iq 
refers to the first possibility (IE.13ap - flE.13bp and the lower sign to the second possibility 
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( ]E.14al) -( 1E.14bl) ). But the other three terms on the right-hand side of ( IF. 1 II) would not 
fit into the right picture of the ortho-configurations. 

But fortunately, there is a third contribution to the kinetic energy ^T kiD (1F.2I) . i.e. 
the azimuthal part ^T^. And indeed, this term is found to be again of the following 
product form: 

wr^Hfi = nc f g (Mffl _ 2cos , £gW -gffl I . f drt , ^ ■ tLw 



4-7T I sin $ 



sin 



(F.12) 



Obviously, it is just this form which is needed in order to bring the transversal kinetic 
energy (^T$ + ^T^) into the expected shape: 



[i±£o)hc 



dQ 

47T 



v 9_m 

sin 3 $ 



+ sin?? • g§($) 



drr 



QJr)-tt(r) 



(F.13) 



or if the angular normalization ( 1VI.131a) is used again 



[b] f# + 



:i± 



o 



)Hcj 



drr' 



njr)-n(r) 



(F.14) 



Observe the difference of this ortho-result with respect to its para-counterpart ( 1VI.32bh ! 

Summarizing, one puts together all three partial results (1F.9P and (IF. 14ft and then 
ultimately finds for the desired kinetic energy ^T kin ( IF. 2ft 



{b) T km ^f km = -^ Idr 



\njr) dftjr) Qjr) dQJr) 



r dr 



r dr 



+ (1 =F 21 



o 



Q + (r)-fL(r) 



(F.15) 



where the upper/lower sign in front of Iq refers to the first /second possibility (1VI.113ap 
f lVI. 1 13bl) / (EmgaD - (IVI. 1 15bl) . 
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Fig. 1: Azimuthal Para- Current ^k^i!)) 

The azimuthal current (VI. 30) is sketched for fixed radial variable r. The current 
strength vanishes on the cones C1/2, C±, C 3 / 2 , • • • where it changes sign, and it becomes 
maximal between the cones, especially on the z-axis (1? = 0, tt). Therefore the magnetic 
field strength H p is singular on the z-axis. For l-p — (i.e. the para-groundstate) the polar 
curve is a circle centered in the origin. The flow direction of the para-current k p = ^k^-e^ 
becomes inverted across the cones ( off-going; Q approaching). 
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Fig. 2: Non-Relativistic Para-Spectrum 

The non-relativistic energy levels of para-positronium can be classified by the quantum 
number i-p of orbital angular momentum (i-p = 0,1,2,3,...). For any (integer) value 
of i-p there exists a subspectrum of excited states which may be enumerated by the 
"principal quantum number" n-p. Degeneracy does occur if states with the same np (but 
different i-p) carry the same energy Et; whether or not this is possible in the electrostatic, 
spherically symmetric and non-relativistic approximation can be clarified not until more 
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exact solutions of the para-eigenvalue problem ( lVI.42l) -( rVT44j) are at hand. For the use of 
approximate variational solutions, see fig.s 4a-4b. Observe also that the discreteness of the 
RST spectrum arises here without any reference to the (anti)-commutation of operators. 
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Fig. 3: Interaction Energy E{j ] = M^c 2 QVI.68\) - fVL7Q) 



The interaction energy flVI.68p - flVI.70p for para-positronium is essentially determined 



by the potential function £pot(^) ( I VI- 70ft and is continuously increasing (up to zero) with 
increasing trial parameter v ( 1VI.55[) . Concerning the values of e vo %{v) for the half- integer 
values of v see the table in Appendix D. Since for l v > the kinetic energy E kin (IVI.65[) - 
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( 1VI.66I) becomes infinite for v — > 0, the total energy E T (» (IVL841) develops a minimum 
(fig. 4) which represents the binding energy of para-positronium. 
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Fig. 4a: Energy Curves Et(V) ( [ VI. 84\ ) 



The energy function Et(V) ( 1VI.84I) has for any quantum number £-p one minimum 
which represents the positronium binding energy. For ip = (n-p = 1) there occurs a 
boundary minimum at v — which then yields exact coincidence of the corresponding 
RST prediction and its conventional counterpart E mw (1VI.85I) . If one is willing to accept 
also weakly singular wave functions (i.e. v < 0) there occurs a local minimum also for l-p = 
0, see fig. 4b. For l-p > the local minima of the energy curves Et(z^) predict the excited 
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positronium energies not better than up to a deviation of some 10% because the chosen 
trial function $(r) ( 1VI.55I) is too simple for the excited states, see the table on p. 11071 
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Fig. 4b: Non- Relativistic Groundstate Energy 

The energy curve Et(V) ( 1VI.84j) for the groundstate [t-p = 0) hits the conventional 



prediction -E con v (1"VI.85j) for vanishing variational parameter {y = 0). This then yields a 
boundary minimum if one wishes to exclude singular solutions which have v < (fig. 4a). 
However if one admits a weak singularity [y < 0) one finds the true (local) minimum of 
the groundstate curve {t-p = 0) at — 7, 23 . . . [eV] for u — — 0, 2049 .... This supports the 
previous prediction of — 7, 68 . . . [eV] which has been obtained by the use of hydrogen-like 
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wave functions, see ref. [12J. This latter method must be considered to yield a more exact 
estimation because it is based upon trial functions with more than the two variational 
parameters (/3, v) due to the present choice $(r) (1VI.55[) . 
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Fig. 5: Angular Factors ^ko^) ( [ Vl.lj.'ty for Ortho-Positronium 



The ortho-configurations do occur in two forms which are characterized by their dif- 
ferent angular distributions [b] k ($) flV1.147al) - flVI.147bl) . The first distribution flV1.147al) 
(broken lines) is characterized by circles with center C at a distance 1/2-k 2 from the z- 
axis. The second distribution ( 1VI. 147b|) (solid curves) equals the first one (lVI.147al) on 
the intersections 5* with the three-cones of vertex angle d(s) = tt/3, 27r/3. 
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For the spherically symmetric approximation both types of charge distributions ' b 'fco( r > i?) 
are replaced by the symmetric density f2(r)/47rr which then generates the spherically sym- 
metric potential ^A (r) according to the Poisson equation (lVI.15ip . By this neglection 
of the anisotropy, there arises the desired ortho/para degeneracy which thus appears as a 
very crude approximation within the framework of RST. 
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Fig. 6: Azimuthal Ortho- Current ^k^tf) ^VI. 15(h ) 

The ortho-current ^k^(r, i?) ( 1VI.150I) factorizes into a radial and a longitudinal part 



%(r,0) = %(r).%(#) , 

quite analogously to its para-counterpart ^/^(r, i?) ( 1VI.30I) . see fig. 1. However, whereas 
the angular part ^^(i?) of the para-current is unique for any value of i-p, the present 
ortho-case fl VI. 1501) for £ = 2 splits up into the two subcases flV1.150ap - flVI.150bD 
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with 



, ' o \ 2 1 , solid line 

| (1 — 16 cos 2 1?) , broken line . 



7T 



Observe here that in the second case (broken line) the ortho-current ^^(r, i?) changes 
sign on the three-cones determined by cos 2 $ = 1/16. 
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